NASA TECHNICAL 
R E PORT 


^e 

■>0 



NASA TR R-466 

A* * 


OAN copy: re^ 
;AFWI- TECHNlCALg 
^>iKIRTLAND AFB^£ 


tin ! 




VIBRATION-TRANSLATION 
ENERGY TRANSFER IN 
VIBRATIONALLY EXCITED 
DIATOMIC MOLECULES 


Robert Laivrence McKenzie 

Ames Research Center 
Moffett Field, Calif. 94035 






■■ i 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


WASHINGTON, D. C • OCTOBER 1976 


TECH LIBRARY KAFB, NM 




TECH LIBRARY KAFB, NM 



0DbB57b 


1. Report No. 

NASA TR R-466 


2. Government Accession No. 


3. Recipient's Catalog No. 


4. Title end Subtitle 

VIBRATION-TRANSLATION ENERGY TRANSFER IN VIBRATIONALLY 
EXCITED DIATOMIC MOLECULES 


5. Report Date 

October 1976 


6. Performing Organization Code 


7. Author(s) 

Robert Lawrence McKenzie 


8. Performing Organization Report No. 
A-6453 


9. Reforming Organization Name and Address 

Ames Research Center 
Moffett Field, Calif. 94035 


10. Work Unit No. 

506-25-61-01-00-21 


11. Contract or Grant No. 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 


13. Type of Report and Period Covered 
Technical Report 


14. Sponsoring Agency Code 


15. Supplementary Notes 

Submitted to York University, Toronto, Ontario, as partial 


fulfillment of Ph.D., April 1976. 


16. Abstract 

A semiclassical collision model is applied to the study of energy transfer rates between a 
vibrationally excited diatomic molecule and a structureless atom. The molecule is modeled as an 
anharmonlc oscillator with a multitude of dynamically coupled vibrational states. Three main 
aspects in the prediction of vibrational energy transfer rates are considered. The applicability 
of the semiclassical model to an anharmonlc oscillator is first evaluated for colllnear encounters. 
Second, the colllnear semiclassical model is applied to obtain numerical predictions of the vibra- 
tional energy transfer rate dependence on the Initial vibrational state quantum number. Thermally 
averaged vibration- translation rate coefficients are predicted and compared with CO-He experimental 
values for both ground and excited Initial states. The numerical model is also used as a basis for 
evaluating several less complete but analytic models. Third, the role of rotational motion in the 
dynamics of vibrational energy transfer is examined, A three-dimensional semiclassical collision 
model is constructed with coupled rotational motion Included. Energy transfer within the molecule 
is shown to be dominated by vibration-rotation transitions with small changes in angular momentum. 

The rates of vibrational energy transfer in molecules with rotational frequencies that are very 
small in comparison to their vibrational frequency are shown to be adequately treated by the preceding| 
colllnear models. 


17. Key Words (Suggested by Author(s)) 

Diatomic molecules 
Vibrational relaxation rates 
Vlbratlon-fotatlon energy transfer 
Semiclassical collision models 


18. Distribution Statement 


Unlimited 

STAR Category - 72 


19. Security Oassif. (of this report) 
Unclassified 


20. Security Classif . (of this page) 
Unclassified 


21. No. of Pages 
242 


22. Price* 
$7.50 


Tor sale by the National Technical Information Service. Springfield. Virginia 22161 





TABLE OF CONTENTS 


Page 

!• Introduction ^ 

1.1 New Aspects in Modern Vibrational Relaxation Processes .... 1 

1.2 Role of Vibration-Translation Energy Transfer 4 

1.3 Purpose and Objectives of This Study g 

1.4 Overview of the Contents and Results 

2. Concepts in Vibrational Energy Transfer • ^5 

2.1 Historical Suiranary and Review Literature Guide 3^5 

2.2 General Considerations 21 

2.2.1 Modes of Energy Transfer in Diatomic Molecules .... 2 I 

2.2.2 Rate Coefficients from the Collision Dynamics 24 

2.2.3 Controlling Variables in Vibrational Energy Transfer • 26 

2.3 Theoretical Methods for Modeling Collision Dynamics 28 

2.3.1 Classical Collision Theories ... 29 

2.3.2 Quantum-Mechanical Theories 3 O 

2.3.3 Semiclassical Collision Model 32 

2.4 Basic Criteria for the Semiclassical Approximation 33 

3. Intramolecular and Collision-Interaction Potentials 37 

3.1 Intramolecular Potential ..... 37 

3.1.1 Vibrational Anharmonicity .......... 37 

3.1.2 Vibration-Rotation Coupling 43 ^ 


iii 



3.2 Collision- Interaction Potential 45 

3.2.1 Spherical Features 47 

3.2.2 Anisotropic Features 51 

3.2.3 Three-Dimensional Interaction Potential Model .... 53 

3.2.4 Collinear Interaction Potential Model 55 

4. A Comparative Evaluation of the Semiclassical Approximation . . ^ 57 

4.1 Semiclassical Model for Collinear Collisions 59 

4.1.1 Incident Particle Motion 61 

4.1.2 Oscillator Motion 63 

4.1.3 Coupling of the Oscillator Motion and the Classical 

Trajectory 65 

4.1.4 First-order Perturbation Solutions ^ 65 

4.1.5 Numerical Solution Methods ^ 67 

4.2 A Comparison with Fully Quantum-Mechanical Solutions . . . ^ 68 

4.2.1 Energy Conservation and the Classical Parameters . . ^ 68 

4.2.2 Influence of Oscillator Response on the Classical 

Motion 74 

4.2.3 Applicability of the First-Order Perturbation Theory . 77 

4.3 Summary 7 g 

5. Vibrational Quantum Number Dependence of Energy Transfer Rates . . 81 

5.1 Collision Model . 82 

5.1.1 Features Influencing the Excited-State Collision 

Dynamics 82 

5.1.2 Aspects of the Semiclassical Numerical Model 85 

iv 



5.1.3 Thermally Averaged Rate Coefficients from a Collinear 

Semiclassical Model 87 

5.2 Analytic Approximations 89 

5.2.1 Keck-Carrier Formula for Anharmonic Oscillators .... 90 

5.2.2 Mies Perturbation Solution for Anharmonic Oscillators . 91 

5.2.3 Kerner Solution for Linearly-Forced Harmonic 

Oscillators 93 

5.3 Comparisons with CO-He Experiments 95 

5.3.1 Effective Interaction Potential Parameters 96 

5.3.2 Comparisons with Excited-State Rate Measurements ... 98 

5.4 An Evaluation of the Analytic Approximations 99 

5.5 Multiple-Quantum Transitions 104 

5.6 Concluding Remarks 107 

6. Effects of Rotational Transitions on Vibrational Energy Transfer. . Ill 

6.1 Vibration-Rotation Collision Model 112 

6.1.1 Interaction Potential 113 

6.1.2 Quantized Molecular Dynamics 114 

6. 1.2.1 Radial Matrix Elements 117 

6.1.2. 2 Spherical Matrix Elements 119 

6. 1.2.3 Complete Factored Matrix Elements 120 

6.1.3 Classical Trajectory 126 

6.2 Representations of the Results 131 

6.3 Numerical Methods of Solution 134 


V 



6.4 Effective Hamiltonian and Other Approximations 136 

6.4.1 Sudden-Rotatlon/Perturbed-Vlbratlon Approximation . . . 136 

6.4.2 Maximum Coupling Approximation 137 

6.4.3 Effective Hamiltonian 138 

6.5 Aspects of Convergence 140 

6.5.1 Evaluation of the Projection-State Decoupling 

Approximations 141 

6.5.2 Convergence Requirements for Vibration-Rotation Energy 

Transfer 144 

6.6 Ambiguities of Unconserved Energy In a Three-Dimensional 

Semiclasslcal Model 151 

6.7 Three-dimensional Inelastic Collisions and Their Relation to 

Colllnear Encounters 154 

6.7.1 H^-He Collisions 156 

6.7.2 N^-He Collisions 161 

6.7.3 CO-He Collisions 171 

6.8 Classification of Rotational Coupling Effects on Vibrational 

Energy Transfer 173 

A Review and Some Considerations for Future Study 177 

7.1 Review of the Newfound Aspects of Vibrational Energy 

Transfer 177 

7.2 Considerations for Future Study 182 



APPENDIX 


A. Notation 185 

B. General Formulation of the Semlclasslcal Collision Theory .... 189 

B.l Classical Trajectory 190 

B.2 Quantum-Mechanical Molecular Motion 192 

B.3 First-Order Perturbation Theory 196 

C. Radial Matrix Element Integral 199 

D. Spherical Matrix Element Integral 209 

E. Analytic Trajectory Equations for Nonzero Impact Parameter .... 215 

F. Symmetries of Vibration-Rotation Matrix Elements and Their 

Compact Computer Storage 219 

F.l Symmetries of 222 

F*l.l Exchange Symmetry 224 

F.l. 2 Projection-State Sign-Reversal S)nnmetry 224 

F.2 Primary Matrix Element Storage 225 

References 229 


vii 



VIBRATION- TRANSLATION ENERGY TRANSFER IN VIBRATIONALLY 


EXCITED DIATOMIC MOLECULES 
Robert Lawrence McKenzie 

V 

Ames Research Center 
CHAPTER I 
INTRODUCTION 

1.1 New Aspects in Modern Vibrational Relaxation Processes 

While the collisional excitation of vibrations in diatomic molecules has 
been a frequently studied topic for decades, an increasing interest in 
processes that depend on the details of energy transfer to specific vibra- 
tional states has put new demands on the analysis of such collisions. The 
following comparison of the early class of relaxation processes with those 
introduced in the past decade demonstrates the new features required in a 
theoretical model of vibrational energy transfer. 

Early studies of vibrational relaxation in gases were concerned mainly 
with the influence of vibrational energy transfer on the bulk thermodynamic 
properties during the relaxation process. Phenomena such as the absorption 
and dispersion of ultrasonic waves^*^ or the vibrational excitation behind 
shock waves^»® were described analytically in terms of an effective **relaxa- 
tion time,” x, that characterized the rate at which the collisional exchange 
of vibrational and translational energies brought the total energy in vibra- 
tions toward equilibrium. Only a single relaxation equation was then 
necessary, in the simple form: 

^Submitted to York University as partial fulfillment of Ph.D., April 1976. 
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( 1 . 1 ) 



dt T 


where e^ is the total specific energy in vibrations and e* is the corre- 
sponding value at thermodynamic equilibrium with the local state of the gas 
(see Appendix A for a list of symbols). The factors controlling these early 
processes are evident in the derivation leading to equation (1-1) (e.g-, 
Vincent! and Kruger,® Ch. 7). The gas is considered to be an admixture of 
harmonic oscillators in a thermal heat bath of structureless inert atoms. The 
harmonic oscillator quantal properties allow only single-quantum transitions 
and lead to the Landau-Teller^ relationship describing the quantum number 
dependence of rate coefficients: 

^v,v-l " ^^1,0 (1*2) 

where denotes the rate coefficient for transitions from oscillator 

state V to v-1 induced by oscillator-atom collisions. A set of simplified 
rate equations describing the detailed kinetics can then be collected into a 
simple equation from which the parameter, x, emerges in the form 

T = [(1 - (1.3) 

where o) is the fundamental oscillator frequency and T, p, and m are the 
kinetic temperature, density, and average molecular mass of the gas mixture, 
respectively. Equations (1-1) and (1-3) combine to illustrate the principal 
common feature of processes described this way; viz, the only collision param- 
eter required is the rate of single-quantum transitions to the ground state. 
The description is independent of the population distributions among higher 
vibrational states and no assumptions regarding their definition is made. In 
a more detailed account of the kinetics that includes the exchange of vibra- 
tional energy between oscillators. Montroll and Shuler® show that a population 
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distribution of harmonic oscillator states rapidly recovers from an arbitrary 
distortion and achieves a Boltzmann distribution described by some nonequilib- 
rium "vibrational temperature." The recovery occurs in a time period small 
compared to t, causing the subsequent relaxation to proceed through a contin- 
uous sequence of Boltzmann distributions. This result further reduced any 
concern for the details of energy transfer to excited vibrational states 
beyond the description given by equation (1.2). 

As experimental studies of vibrational relaxation became more detailed, 
the kinetic models based on harmonic oscillator properties appeared less 
capable of describing the observations. Motivated by some large discrepancies 
between theory and experiment in nonequilibrium supersonic expansions, 

Treanor, Rich, and Rehm^® recently showed that the small anharmonicity of most 
diatomic molecules was sufficient to generate non-Boltzmann distributions 
among upper vibrational states during some relaxation processes. The impor- 
tance of oscillator anharmonicity in the vibrational kinetics has since been 
amplified by the introduction of infrared gas lasers, where the effects 
of anharmonicity are essential to produce the necessary population inversions 
among vibration-rotation states- Some recent proposals using lasers to 
selectively excite specific vibrational states for photochemical or isotope 
separation experiments will also be strongly influenced by the effects of 
anharmonicity. These modern applications of vibrational nonequilibrium con- 
stitute a new class of relaxation processes that depend on the degree of dis- 
tortion from a Boltzmann population distribution. Their analysis requires a 
detailed solution to the set of relaxation equations — one for each contribut- 
ing vibrational state — that describes the change in number density, N^, of 
each state v. A general form of the relaxation equations can be written as 
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dN. 


dt YjYj K' 


(1.4) 


1 V 


where N_* denotes the number density of all species or states i and k , 

V ,v 

is the rate coefficient for transitions from state v' to v. Equation (1.2) 
is not an accurate description of the quantum number dependence of k^i ^ 
when the oscillator is anharmonic, and the selection rule allowing only 
single-quantum transitions is also invalid. Thus the simplifications of 
equation (1.4) leading to equation (1.1) and the concept of an effective 
relaxation time given by equation (1.3) are no longer applicable. Rate coeffi- 
cients for transitions from excited vibrational states are as essential to the 
analysis of such processes as k^ q is and vibrational anharmonicity will have 
a major influence on their values. Since this study is concerned mainly with 
the rate of energy transfer to excited states, oscillator anharmonicity is 
therefore a basic feature to be included. 


1,2 Role of Vibration-Translation Energy Transfer 

In most modern applications involving vibrational nonequilibrium and 
particularly in the analysis of infrared gas lasers, the collisional exchange 
of vibrational and translational energies must usually be considered as one of 
several paths for energy transfer to the molecular state in question. As 
equation (1.2) indicates, the rate of vibration-translation (V-T) energy 
transfer increases with quantum number even in the simplest model of the 
oscillator. Thus, the V-T process can dominate the flow of energy from upper 
vibrational levels even where it may be insignificant to the kinetics of lower 
levels. In some cases, it may provide the principal path for vibrational 
energy loss from the system. An essential feature of the V-T rates is there- 
fore their dependence on the initial state quantum number, particularly in 
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deexcitation processes where excess vibrational energy has been produced. 
Unfortunately, very little quantitative information defining the quantum num- 
ber dependence of V-T rates for even the simplest diatomic molecules is 
presently available from either experiment or theory. 

Experimental ground-state excitation rates have been obtained from mea- 
surements behind shock waves or in fluorescence experiments^^ for many years 
by determining the value of x in equation (1.1) that best fits the observa- 
tions. However, the difficulty of obtaining experimental V-T rates for mole- 
cules in well-defined excited vibrational states is indicated by the sparsity 
of attempts. Numerous experimenters have recently measured the rates of 
vibration-vibration (V-V) energy exchange between pairs of oscillators in 
excited states^®’ because the fast V-V transfer can easily be made a dominant 
mechanism; but to date, only one comprehensive set of upper-level V-T rate 
measurements has been reported.^® Even then, while the experiment was cleverly 
designed and carefully analyzed, the conditions were complex and the measure- 
ments required substantial correction to compensate for extraneous modes of 
energy transfer. 

Theoretical studies addressed to the analysis of initially excited oscil- 
lators have been similarly sparse. Since the relaxation time, x, is deter- 
mined solely by the rate of single-quantum transitions to the ground state, 
kj Q, the usual theoretical approach has centered on a harmonic oscillator 
model of the molecule initially in the ground state. The small-amplitude 
oscillations characteristic of the ground state then allow a linearized inter- 
action potential between the oscillator and an incident particle to be used. 

By assuming further that the particle trajectory is collinear with the molecu- 
lar axis and by adopting a semiclassical approximation, the linearized 
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interaction makes possible an exact and convenient analytical solution for the 
oscillator transition probabilities for any initial state. However, the 
inaccuracy of the harmonic oscillator model has been demonstrated by Mies^^“"^^ 
even for transitions originating from the ground state. Mies found that the 
use of an enharmonic oscillator potential introduces matrix elements associated 
with oscillator transitions that are no longer equal on the diagonal. (A har- 
monic oscillator with an interaction linear in the oscillator coordinate has 
constant diagonal matrix elements.) The nonzero differences in the diagonal 
matrix elements introduce additional phase differences between the time- 
dependent oscillator eigenfunctions during a collision and can lead to large 
corrections to the harmonic oscillator model. Because the origin of these 
corrections resides in the unperturbed oscillator eigenfunctions (from which 
the matrix elements are computed) , their effects are not always reproduced by 
the popular practice of simply inserting oscillator eigenenergies , corrected 
for anharmonicity , into a harmonic oscillator theory. Nevertheless, in the 
absence of better analytic solutions, such theoretical models are frequently 
used to predict the quantum number dependence of V-T rates. Thus, a need 
clearly exists for the development of a suitable analytic solution containing 
anharmonicity as a fundamental feature. 

1.3 Purpose and Objectives of This Study 

While the modern literature is abound with comprehensive and detailed 
studies of the collisional excitation of diatomic molecules in vibration and 
rotation (see almost any recent issue of the Journal of Chemical Physics) , 
vibrational states higher than the second are rarely considered. The objec- 
tives are usually either to examine improved techniques for calculating the 
collision dynamics or to obtain a quantitatively accurate estimate of the 
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inelastic cross sections. The collision models used in the latter case are 
usually fully quantum mechanical and hence are as exact as the form of the 
interaction potential chosen for study. These studies are clearly a necessary 
step in the understanding of molecular collision dynamics since they provide 
the most precise test of our ability to explain the experimental observations. 
Unfortunately, the computational requirements to obtain such precision are 
expensive and tend to limit the scope of such studies. To a pragmatist con- 
cerned with the analysis of a modem macroscopic process, these studies of 
microscopic collision dynamics are seldom able to provide much useful informa- 
tion about the thermally averaged rate coefficients for molecules in excited 
states. Furthermore, even if exact calculations were typically carried far 
enough to produce rate coefficients, a means of numerically reproducing the 
results inexpensively would be required before they could be conveniently 
applied in a solution of the macroscopic rate equations. This study is 
addressed to the pragmatist and to four corresponding objectives. 

The first objective is to examine the inelastic collision dynamics of 
diatomic molecules in an arbitrary initial state struck by a structureless 
atom. The purpose is to explore the qualitative features of such encounters 
and to identify the parameters and physical features contributing most to the 
prediction of the associated energy- transfer rates. By including a complete 
account of the coupling between interacting vibrational states, especially as 
it is amplified by anharmonicity , the results provide a basis for evaluating 
more approximate treatments of the collision process that may have emphasis on 
other aspects, such as those with sufficient simplifications to allow analytic 
solution or those including coupled rotational motion. 
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A second objective of this study is to evaluate several analytic colli- 
sion models in popular use for predicting the quantum number dependence of V-T 
rate coefficients. This objective is motivated by the importance of having an 
inexpensive means of generating values of when solving the detailed 

rate equations typified by equation (1.4). A collision model with sufficient 
generality to be applicable for all conditions of interest in modern applica- 
tions will necessarily require numerical solution, and the first objectives of 
this study are met only with such a model. The solutions are time consuming, 
however, and would be prohibitively expensive in a practical application. 
Consequently, the usual practice is to obtain a simple analytic approximation 
by introducing sufficient assumptions to decouple the interactions between 
oscillator states and to linearize the interaction between the oscillator and 
its collision partner. Several such solutions have been extracted, but they 
all exclude one or more properties of the collision process that remain impor- 
tant in a generalized model. Little definition, if any, of the range of 
applicability of these analytic solutions appears in the literature in other 
than the most general terms. The second objective described here is an effort 
to define more explicitly their range of applicability for predicting excited- 
state rate coefficients. 

A third objective of this study is to evaluate the consequence of 
several simplifying assumptions regarding the equations of motion and the col- 
lision geometry that were necessary to meet the preceding objectives. A fun- 
damental simplification to the equations of motion is achieved here by adopt- 
ing a semiclassical or ^’impact parameter” description of the collision 
dynamics . The path of the incident particle is obtained classically, 
while the oscillator response is treated quantum mechanically. As a result. 
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the second-order quantum mechanical equation of motion Is reduced to two first- 
order differential equations, that are subsequently decoupled. The complexity 
of solution is thereby reduced greatly. However, the semiclassical approxima- 
tion fails to properly account for the quantal interference between colliding 
nuclei while the decoupling of first-order equations obviates the conservation 
of total energy. These shortcomings have been partially compensated in similar 
harmonic oscillator models, but the success of the compensations has not been 
tested for enharmonic oscillators. Part of this third objective is to examine 
and define the limitations of the semiclassical approximation when applied to 
enharmonic oscillators. The results will contribute to a more complete under- 
standing of the associated analytic solutions that are also based on the semi- 
classical approximation. 

A final objective of this study is to evaluate the influence of coupled 
rotational motion on the rate of vibrational energy transfer. The necessity 
of including a multitude of vibrational states with large quantum numbers to 
study their interactions required a reduction elsewhere in the complexity of 
the molecular motion to keep the problem within practical bounds. The obvious 
choice was to eliminate any account of the rotational motion by limiting the 
collision geometry to one-dimensional collinear encounters. This procedure is 
commonly applied throughout the literature for similar reasons and is usually 
based on the presumption that collinear collisions are the most effective for 
inducing vibrational excitation. However, intuitive notions suggest that the 
inelasticity of three-dimensional collisions is partitioned among vibrational 
and rotational degrees of freedom in the molecule, in varying amounts, depend- 
ing on the molecular inertial properties and the initial state. Kelly and 
Wolfsberg^^ have used a fully classical model to demonstrate that collinear 
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collisions are not always the most effective for converting vibrational energy. 
For example, atom collisions with molecules possessing widely spaced rotational 
states can induce vibrational transitions with very little energy converted to 
translational motion. Vibration-rotation energy transfer is then a more cor- 
rect description of the event. On the other hand, a full account of the 
coupled vibration-rotation motion in a three-dimensional collision model must 
include at least all of the energetically accessible rotational states in each 
vibrational level. As a result, collision energies sufficient to induce vibra- 
tional transitions will encompass hundreds of multiply degenerate rotational 
states in most molecules. Since the occupation of each state must be treated 
as a separate variable, an extremely large system of coupled differential 
equations is required whose numerical solution is intractable for all but the 
simplest cases. Fortunately, some methods have recently been introduced that 
average the combined action of degenerate states and reduce the problem to an 
expensive, but tractable, size. Several of these methods are evaluated and 
applied here, in conjunction with a three-dimensional semiclassical collision 
model. From another point of view, the objective is to determine the conse- 
quences and validity of using a collinear one-dimensional model to predict the 
vibrational quantum number dependence of vibrational energy transfer rates. 

The results give considerable credibility to the preceding conclusions of this 
study derived from collinear models. 

1.4 Overview of the Contents and Results 

In chapter 2, the basic concepts and assumptions commonly applied in the 
analysis of vibrational energy transfer are reviewed. A brief historical 
review of vibrationally inelastic collision models is first presented that 
provides a guide to a number of more complete review papers. The concepts 
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and physical factors that control the rate of energy transfer are then dis- 
cussed, and the several theoretical approaches from which the rates may be 
estimated are evaluated as they apply to the objectives of this study. The 
semiclassical approximation is shown to be the most suitable approach to the 
objectives previously stated. 

Having established the primary theoretical approach to be one in which 
the incident particle path is computed classically, chapter 3 examines the 
features of the interaction potentials that determine both the internal mole- 
cular motion and the classical trajectory. Arguments are presented showing 
that a Morse-oscillator/rigid-rotor description of the molecule is adequate 
for the purposes of this study and that the classical trajectory may be deter- 
mined from just the short-range repulsive forces between colliding nuclei. 

In chapter 4, a collinear, semiclassical, collision model for predicting 
V-T transition probabilities from arbitrary initial states is developed and 
evaluated. Comparisons are made with equivalent, exact, fully quantum mechani- 
cal solutions obtained from the literature for a broad range of collision 
parameters, molecular types, and initial states. While similar comparisons 
have been made before, they have been less complete and limited to harmonic 
oscillator models of the molecule initially in the ground state. This work 
includes a more extensive variation of collision parameters and tests the 
application of the semiclassical approximation to an enharmonic Morse oscil- 
lator in several elevated initial states. In the past, there has been a 
variation of opinions on the best method of compensating for the lack of 
energy conservation in the semiclassical approximation. The comparisons of 
this study show that the correction is nearly the same from all methods sug- 
gested and no clear choice of the best method is possible — nor is a choice 
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necessary for they are all adequate in the range of collision energies of 
interest. Finally, the comparisons in chapter 4 bring to light the result 
that, while the semiclassical approximation works well for a broad range of 
collision parameters when the molecule is treated as a harmonic oscillator, the 
more realistic enharmonic oscillator model imposes some definite limitations. 
When the anharmonic oscillator is homonuclear and struck by a collision partner 
whose mass is less than either molecular nucleus, the semiclassical approxima- 
tion is very successful. However, its application to heteronuclear molecules 
or to homonuclear molecules struck by a heavy collision partner produces 
anomalous resonances that do not appear in an equivalent harmonic oscillator 
model. These anomalies are partially eliminated when additional coupling 
between the oscillator and the incident particle is introduced. 

With the limitations of the semiclassical approximation established for 
anharmonic oscillators, chapter 5 describes an investigation of the factors 
that influence the prediction of V-T rates for initially excited molecules. 

The capability of several analytic theories for reproducing rate coefficients 
predicted by a more exact numerical model is also evaluated. Unfortunately, 
the most widely used and simplest analytic formula also produces the poorest 
estimate of quantum number dependence. But two slightly less convenient 
analytic models are found to reproduce the more exact predictions for well- 
defined and easily identifiable ranges of conditions. Both favorable analytic 
models are based on a collinear semiclassical description of the collision. 

The validity of the collinear collision models used in the previous 
chapters is evaluated in chapter 6 using a three-dimensional semiclassical 
model developed for that purpose. A complete model is first constructed that 
allows an arbitrary number of coupled vibration-rotation states to be included. 
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It provides a basis for evaluating the accuracy of some approximate 
formulations in which the combined effects of the degenerate projection states 
associated with each rotational quantum state are decoupled and treated col- 
lectively. The ability to decouple the degenerate states makes the objective 
of this part of the study possible. An "effective Hamiltonian" approximation 
is found to be the most useful, and it is applied to a study of the influence 
of rotational energy transfer on the rate of vibrational excitation. The 
results show that the effects of rotations can be segregated into three 
classes. For molecules like hydrogen or the hydrogen halides that have a 
rotational frequency only a magnitude smaller than the fundamental oscillator 
frequency, the rotational coupling is large and energy transfer can proceed 
via rotation- vibration transitions with very little conversion of translational 
energy. The behavior of molecules with these properties is further separable, 
depending on the initial angular momentum, but the use of a collinear colli- 
sion geometry is physically unrealistic in any case and the corresponding 
analytic rate formulas are consequently of little value. On the other hand, 
the third and much larger class of molecules, in which a multitude of rota- 
tional levels is contained in each vibrational state, is not influenced by the 
accompanying rotational motion induced in a three-dimensional collision. Rate 
coefficients obtained from a collinear collision model then reproduce all of 
the physical features contained in the three-dimensional results and predict 
an essentially identical dependence on vibrational quantum number when com- 
pared with the net vibrational transition rate summed over all final rotational 
states. Correspondingly, the predictions for such molecules are also shown to 
be insensitive to the initial rotational state of the molecule. These results 
lend considerable credibility to the results from collinear collision models 
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and to the analytic solutions that depend on them. Finally, chapter 7 
summarizes the new aspects of the results of this study and presents some con- 
siderations for additional theoretical and experimental study. 
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CHAPTER 2 


CONCEPTS IN VIBRATIONAL ENERGY TRANSFER 

In the preceding chapter, methods for describing the macroscopic behavior 
of a vibrational relaxation process were discussed. The microscopic aspects 
of the process were contained in a thermally averaged rate coefficient, , 

of undefined nature. The remainder of this study concentrates on the physical 
factors that affect theoretical models used to evaluate it. 

This chapter provides a general discussion of the concepts leading to a 
theoretical model. A brief historical review is given first that provides a 
commentary on some pertinent publications describing the various concepts in 
greater detail. The latter part of this chapter defines the controlling 
d 3 mamic and molecular parameters affecting and reviews the general con- 

siderations leading to a choice for the fundamental theoretical approach to be 
applied in the remainder of this study. 

2.1 Historical Summary and Review Literature Guide 

In the 1930 's and before, the anomalous absorption and dispersion of 
ultrasonic waves propagating in gases were the principal phenomena motivating 
the study of vibrational energy transfer. All fluids absorb ultrasonic waves 
through shear viscosity losses and, in most cases, through heat conduction. 

By those mechanisms, they all display a corresponding dispersion- However, 
molecular fluids (i.e., those with a capacity for internal-energy storage) 
have an additional absorption and dispersion originally accounted for by 
heuristically introducing a "bulk viscosity" into the Navier-Stokes equations 
that describe the process. Investigators soon recognized, however, that the 
artificial viscosity was a manifestation of internal-energy absorption in the 
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molecule. As early as 1928, Herzfeld and Rice^ explained the origin of the 
additional absorption and dispersion conceptually in terms of collisional 
energy transfer between the translational and internal degrees of freedom at a 
finite rate. A few years later, Oldenberg^^ discussed molecular collisions 
qualitatively to show the inelasticity of rotational and vibrational motion. 

In 1931, Zener was the first to give a detailed mathematical treatment for 
vibrationally inelastic collisions. His theory was based on a quantum mechan- 
ical perturbation method, referred to as a "distorted wave approximation," 
applicable to low collision energies where the transition probabilities are 
small. Then, in 1936, Landau and Teller^ published their historic paper in 
which the properties of the rate coefficient were explored, again from a more 
conceptual point of view. They used partially intuitive arguments (with no 
reference to the earlier work) to show that equation (1.2) described the rudi- 
mentary dependence of ^v,v-l ^ they obtained the equally important 

dependence of on temperature, given by log k^ T“^/^. 

Later, interest shifted away from theoretical work addressed to ultra- 
sonic relaxation phenomena, but increased in the study of a very similar 
inelastic molecular collision problem. The accommodation coefficient, related 
to the energy transfer between gas molecules and a solid surface, was studied 
extensively, first by Jackson and his coworkers^^ and later by Lennard-Jones 
and his coworkers. The paper by Jackson and Mott^^^ became particularly 
noteworthy because it provided a simplified mathematical derivation of the 
"distorted wave approximation" still referred to in modern texts. 

Only a few contributions to the field of vibrational relaxation followed 
until the early 1950’ s. For example, one-dimensional treatments by Zener and 
other early investigators were extended to three-dimensional collisions. 
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After 1950, Takayanagi^^ introduced the "modified wave number" approximation 
designed to reduce the numerical labor in three-dimensional problems and 
applied it to rotational transitions in hydrogen. Meanwhile, Schwartz, 

Slawsky, and Herzfeld^^ published their well-known paper in which vibrational 
transitions were treated with the distorted wave approximation. Their closed- 
form analytic formulas for resonant and nonresonant transition probabilities 
in collinear collisison have become the most widely used means of estimation 
until recent times. Their formulation, often referred to as the SSH theory, 
was later extended to three-dimensional collisions but for a nonrotating 
molecule. 

Experimental methods for measuring vibrational relaxation times were also 
developing. Improved rate data were obtained from measurements in jets^^»^^ 
and behind shock waves^^»^^ at other than room temperatures- The existing 
perturbation theories were not always applicable to the analysis of these new 
experimental techniques, however, because the theories were limited to low- 
energy collisions pertaining mostly to near- room temperatures. In 1958, 
Kerner^^ obtained a nonperturbative exact solution to the Schrodinger equation 
for a harmonic oscillator in the presence of a time-dependent forcing function. 
The only constraint on the forcing function was that it be linear in the 
oscillator coordinate. Kerner’s solution was subsequently applied by Treanor^^ 
in a semiclassical treatment of high-energy collinear collisions, thereby 
achieving an analytically exact formula in closed form for the transition 
probabilities of a harmonic oscillator at all collision energies. 

The 1960 ’s brought on a deluge of publications concerned with vibrational 
relaxation phenomena that has persisted to this day. So much experimental 
information became available that Millikan and White^® were able to correlate 
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empirically vibrational relaxation times for a large number of diatomic mole- 
cules in terms of the fundamental oscillator frequency and the reduced mass of 
the collision pair. Their correlation was only modestly guided by theory, how- 
ever. During this period, theoretical studies were stimulated by a rapidly 
growing computer technology. Consideration of exact numerical solutions to the 
collision problem became a reasonable occupation; but vibrationally inelastic 
solutions were awkward until Secrest and Johnson^^ developed a numerical method 
of "amplitude density functions” that allowed one-dimensional scattered wave 
functions. to be obtained efficiently. Their methods have since been extended 
to treat three-dimensional collisions with vibrational and rotational 
inelasticity. 

The early 1960 *s also marked the appearance of some review articles of 
modern interest that describe the various theoretical approaches in detail. 

One of the first was the chapter by Herzfeld^^ on "Relaxation Phenomena in 
Gases.” His discussion is based primarily on the application to ultrasonic 
absorption and dispersion, but he gives a clear account of the early theoreti- 
cal approaches in which the fundamental collision parameters are described. 

He later provided a more complete account in textbook form. ^ Cottrell and 
McCoubrey^^ took a slightly more modern approach in their book by dealing 
with the quantum mechanical aspects of the collision process in greater detail. 
However, their discussion remains physically descriptive and valuable as an 
introduction to the theoretical approximations leading to analytic solutions. 
Takayanagi^^ provided the first comprehensive review of the theoretical aspects 
of vibrationally and rotationally inelastic collisions covering the period up 
to 1963. With the rapid developments in the field following 1963, Takayanagi^^ 
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published a second » equally comprehensive, review covering the developments to 
1965. 

The reviews published before 1965 preceded the time when exact quantal 
solutions for vibrationally inelastic collisions were obtainable with suffi- 
cient ease and confidence to serve as a basis for testing the approximate 
methods. Prior emphasis was directed toward the comparison of approximate 
theories with experiment as a test of their validity. A more recent review by 
Rapp and Kassal**^ was therefore addressed, in part, to an evaluation of the 
earlier theories by comparing them with exact numerical solutions. With the 
greatly increased detail of Information about the collision dynamics in view, 
both from fully quantum mechanical and fully classical numerical solutions, 
Rapp and Kassal evaluate many of the assumptions contained in the approximate 
collision models and provide a useful guide to their range of applicability. 
Their article also deals with some aspects of modem interest such as the 
transfer of vibrational energy between oscillators and the effects of oscil- 
lator anharmonicity. However, at the time of Rapp and Kassal*s writing, the 
new class of vibrational relaxation processes had not quite impacted the 
theoretical community. Their emphasis therefore centered on the dynamics of 
oscillators in or near the ground state. 

The beginning of this decade brought in widespread efforts to deal with 
the new requirements in the analysis of vibrational energy transfer. Rich and 
Treanor^^ presented a comprehensive review of the aspects of vibrational 
relaxation in gasdynamic flows. Their discussion is devoted mainly to non- 
equilibrium flow processes and hence to an application of vibrational rate 
theories, but they also provide a detailed description of many aspects of 
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vibrational energy transfer that motivated this study and their review serves 
as an introduction to the new class of relaxation processes. 

With numerical investigations now a practical and popular approach to 
molecular physics, this decade begins the era in which molecular collisions 
can be studied in much greater detail. Many of the new methods in "numerical 
physics" are described by Secrest^^ in a recent review of their application to 
rotational and vibrational energy transfer. Most methods are motivated by the 
need to reduce the numerical labor and expense. For example, activity in the 
use of fully classical calculations for reactive and vibrationally inelastic 
collisions has flourished with the development of "Monte Carlo" or random- 
selection methods for averaging the results of many trajectories and orienta- 
tions.^® Another new concept is the semiquantal approximation, developed 
independently by both Miller^® and Marcus.®® Also termed the "classical 
S-matrix" theory, the principal distinction of the semiquantal approach from 
earlier semiclassical methods is that, in the former, all degrees of freedom 
are treated classically but with the quantum-mechanical principle of super- 
position subsequently applied. While the method allows pure quantum effects 
such as tunneling, selection rules, and interference to be studied and has 
given reasonably accurate predictions for vibrationally inelastic collisions,®^ 
its limitations are not yet fully charted. Two reviews of the theory have. ; 
recently been published by Miller.®^’®® 

A final milestone that has contributed significantly to the results of 
this study is the success of several efforts to average the combined contribu- 
tions of degenerate rotational states and thereby make the study of vibration- 
rotation interactions feasible. Forthcoming discussions in this study demon- 
strate the futility of a complete treatment of the problem. The first 
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successful solution to the problem was presented by Rabitz, who formulated 
an "effective Hamiltonian" that nullifies the contribution of individual pro- 
jection quantum states before the equations of motion are solved. Following 

his work, McGuire and Kouri^^ proposed a "j -conserving" approach of some simi- 

z 

larity to the effective Hamiltonian approximation. A different method of 
reducing the rotational aspects of the problem has been studied by Pack and 
his coworkers, who treat the rotational motion in a "sudden approximation," 
well known in its basic form from numerous modern textbooks on quantum mechan- 
ics, The relationship of all of these methods for decoupling the internal 
angular momentum of the molecule has recently been examined by Secrest.^^ 

These methods and their application represent a large part of the current 
activity in studies of vibrationally and rotationally inelastic collisions. 

2 . 2 General Cons ider at ions 

Considered in the following paragraphs are some of the general concepts 
and controlling parameters that form the basis of most theoretical models for 
collisions involving a diatomic molecule. The general features of several 
theoretical approaches are then reviewed to guide the choice of a method best 
suited to the intentions of this study. 

2.2,1 Modes of Energy Transfer in Diatomic Molecules 

Binary collisions involving a diatomic molecule are not yet treated in 
general terms entirely from first principles. Ab initio approaches to the 
many-body problem presented by three or more nuclei and their attendant elec- 
trons are still intractable on present-day computers. Fortunately, when elec- 
tronic transitions are not of concern, an adequate collision model does not 
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require an explicit description of the coupled nuclear and electron motion but 
Instead relies on the nearly instantaneous adjustment made by the electrons to 
the nuclear motion. The problem then reduces to one of describing only the 
dynamics of the nuclear motion during a collision but requires Interaction 
potentials independently obtained by some less rigorous means. 

In this study, our interest is further restricted to collisions only of a 
diatomic molecule with a heavy structureless particle such as an atom in its 
ground electronic state. The inelasticity of the collision is then confined 
to the internal rotational and vibrational energy modes of just one molecule 
and we avoid the complexity of dealing with the exchange of internal energy 
between molecules. Energy transfer still occurs to any of several internal 
inodes, however, as figure 2.1 illustrates. A characteristic of most diatomic 
molecules is the widely separated vibrational eigenenergies (heavy line levels 
in fig. 2.1), each with a manifold of closely spaced rotational states (light 
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Figure 2.1.- Modes of energy transfer in a diatomic molecule. 
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line levels in fig. 2.1). Thus, a collision with sufficient energy to excite 
vibrational motion in the molecule will simultaneously excite many rotational 
states as well. The arrow labeled V-R-T (vibration-rotation-translation) in 
figure 2.1 typifies those kinds of transitions. The internal energy change 
appearing in the molecule will be reflected as a change in the translational 
energy of the colliding pair. While a complete description of V-R-T energy 
transfer is complex, certain limited paths for energy transfer are often the 
dominant mechanism and they can then be treated separately. For example, 
near-resonant transitions between vibration- rotation states (V-R in fig. 2.1) 
may be dominant in some molecules with a suitable initial condition. In this 
case, any energy traded with translation appears only as an elastic deflection 
after the encounter. In another situation, the small amount of energy required 
to induce a rotational transition within the same vibrational state makes the 
exchange of rotational and translational energies (R-T in fig. 2.1) probable 
at collision energies where the vibrational state of the molecule may be 
ignored. The molecule is treated as a rigid rotor in such circumstances. The 
analyses of these limited cases involving rotation are usually simpler than 
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(a) Three-dimensional encounter. 
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(b) One-dimensional collinear encounter. 

Figure 2.2.- Collisional motion. 
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cannot be applied in a collinear collision model because b or are not 



the general case, but a three-dimensional collision geometry is still required 
as shovm in fig. 2.2(a). In contrast, an even simpler but more restricted 
treatment of the collision is one in which no rotational transition occurs. 
Translational energy is exchanged only with vibration (V-T in fig. 2.1). 

Events of this nature can occur in a three-dimensional encounter of arbitrary 
orientation because transitions always take place with a probability less than 
unity (making no rotational transition also probable) , but a frequently used 
approach that drastically reduces the complexity of the problem is to assume 
that the most efficient producers of vibrational transitions are collisions 
with trajectories along the intranuclear axis of a nonrotating molecule. The 
corresponding one-dimensional collinear geometry is illustrated in fig- 
ure 2.2(b). As it turns out, we shall find in this study that the collinear 


km , = ( 8 kT /.„)'/2 j ( 2 . 2 ) 

•'o 

where unsubscripted k is Boltzmann’s constant and y is the reduced mass of 
the collision. If m^^ denotes the mass of nucleus i, then, using the nota- 
tion in figure 2, 

“c> 

y = 7 T T 7 (2.3) 

(m^ + mb + m^,) 

Methods of performing the integrals in equations (2.1) and (2.2) are 
discussed in greater detail in chapters 5 and 6 where specific applications 
are made. At this point, one only needs to recognize that obtaining the tran- 
sition probability, P(E,b), from a semiclassical treatment or the differential 

v^v’ 

cross section, da/dS^, from a full quantum mechanical treatment is the funda- 
mental problem. Once either of these results is obtained, a calculation of 
the corresponding rate coefficients is relatively straightforward. 

2.2.3 Controlling Variables in Vibrational Energy Transfer 

For a collision to affect the vibrational motion of an oscillator, the 
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the general case, but a three-dimensional collision geometry is still required 
as shown in fig. 2.2(a). In contrast, an even simpler but more restricted 
treatment of the collision is one in which no rotational transition occurs. 
Translational energy is exchanged only with vibration (V-T in fig. 2.1). 

Events of this nature can occur in a three-dimensional encoimter of arbitrary 
orientation because transitions always take place with a probability less than 
unity (making no rotational transition also probable), but a frequently used 
approach that drastically reduces the complexity of the problem is to assume 
that the most efficient producers of vibrational transitioiis are collisions 
with trajectories along the intranuclear axis of a nonrotating molecule. The 
corresponding one-dimensional collinear geometry is illustrated in fig- 
ure 2.2(b). As it turns out, we shall find in this study that the collinear 
collision model is surprisingly useful for predicting the rate of vibrational 
energy transfer when many rotational states are associated with each vibra- 
tional mode. Unfortunately, the incomplete nature of the collinear model 
makes it awkward to obtain a rate coefficient from the detailed collision 
dynamics calculated, whereas rate coefficients evolve naturally from the 
results of a three-dimensional model, as the following discussion demonstrates. 

2.2.2 Rate Coefficients from the Collision Dynamics 

The relationships between a thermally averaged rate coefficient and the 
results of a microscopic collision model depend somewhat on the theoretical 
description used to model the collision. For example, a semiclassical formu- 
lation in which the incident particle trajectory is obtained from the classical 
equations of motion deals with collisions characterized by the parameters E 
and b, where E is the initial, relative, kinetic energy in a center-of-mass 
reference frame and b is the impact parameter measuring the minimum distance 
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between mass centers of each collision partner that would occur if the rela- 
tive path were not deflected. The outcome is a transition probability, P(E,b) 
between each pair of internal quantum states, v and v’ , included in the molecu- 
lar model. A necessary step in obtaining the corresponding rate coefficient is 
to first produce the total cross section, o(E) . In the semiclassical frame- 
work, a total cross section may be generated by repeating the calculation for 
a sufficient range of impact parameters to evaluate the integral 

r 

a(E) = P(E,b) 27Tb db (2.1) 

v>v’ I v^v’ 

•'o 

Likewise, a fully quantum mechanical description of both the molecule and 
incident particle motion follows a similar procedure. However, the incident 
particle path is not localized so the collision must be characterized instead 
by E and a final scattering direction, 0,. The outcome is a differential 
cross section da/d^, which may be computed for all elemental solid angles and 
integrated over the sphere by an expression similar to equation (2.1). In 
contrast to any three-dimensional method, equation (2.1) or its equivalent 
cannot be applied in a collinear collision model because h or ^ are not 
included variables. Thus, obtaining a(E) from a collinear description 
requires some kind of "steric factor" on an "effective hard-sphere cross sec- 
tion" to be introduced. These additional artifacts are discussed in greater 
detail in chapters 5 and 6 where collinear collision models are evaluated. 

Once the cross section is in hand, the desired rate coefficient for 
V v' transitions is obtained by averaging the energy-dependent total cross 
sections over a thermal energy distribution characterized by a kinetic temper- 
ature T. Remembering that E is the relative kinetic energy in a center-of- 
mass frame, the result is® 
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v^v' 

where unsubscripted k 
the collision. If 
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is Boltzmann's constant and y 
denotes the mass of nucleus i, 



is the reduced mass of 
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Methods of performing the integrals in equations (2.1) and (2.2) are 
discussed in greater detail in chapters 5 and 6 where specific applications 
are made. At this point, one only needs to recognize that obtaining the tran- 
sition probability, P(E,b), from a semiclassical treatment or the differential 

v^v’ 

cross section, da/d^2, from a full quantum mechanical treatment is the funda- 
mental problem. Once either of these results is obtained, a calculation of 
the corresponding rate coefficients is relatively straightforward. 


2.2.3 Controlling Variables in Vibrational Energy Transfer 

For a collision to affect the vibrational motion of an oscillator, the 
disturbing force created by the impacting particle must vary in a time period 
that is comparable with or less than the normal oscillator period. Otherwise, 
a slowly applied disturbance simply allows the oscillator to adjust adiabati- 
cally, leaving its final condition unaffected by the encounter. Conversely, 
an impulsively applied force will severely disturb the phase of the oscillator 
motion and efficiently upset the proportions of energy in vibration, rotation, 
and translation. These conditions can be expressed more explicitly by denot- 
ing a representative time interval in which the colliding pair interact as 

and letting represent the fundamental oscillator frequency. We then 
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argue that energy transfer will occur with increased efficiency as becomes 

less than the oscillator period 1 /Vq, that is, 0. One way to evaluate 

the magnitude of is by noting that the interaction distance for collisions 

influenced mainly by repulsive forces can be loosely measured in terms of a 
range parameter L, Then ^ 2L/u, where u is the average relative colli- 
sion speed. If the relative kinetic energy is E and y is the reduced 
mass, then the relative collision speed is /(2E/y) and the efficiency of 
energy transfer will increase as 



We should therefore expect the quantities related to the efficiency of energy 

transfer, such as P(E,b) and a(E) , to increase with E and to decrease as 

v^’ v->v^ 

the oscillator frequency, interaction range, or nuclear masses are made 

larger. Note that the impact parameter, b, could also serve as a measure of 

the interaction range when it exceeds L, Collisions at increasing b will 

become more adiabatic with an accompanying decrease in P(E,b). 

v->v’ 

The preceding relationships are sometimes described in terms of an 
**adiabaticity parameter** 

The larger g becomes, the more adiabatic is the encounter. Values of ? 
below unity lie in the **sudden’* region. Similar ideas are applied to rota- 
tional motion as well simply by replacing with the fundamental rotational 

frequency, v^. Note that in many molecules, so that collisions in 

some energy range may be adiabatic with regard to vibrations while sudden in 
respect to rotations. 
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The effects of varying ? can be made slightly more quantitative by con- 
sidering its sequel, the "resonance function," R(C)» a measure of the effi- 
ciency of energy transfer. If AE(?) is the average energy transferred in a 
collision characterized by 5, we can define 

R(0 = AE(?)/AE(5 = 0) (2.6) 

where AE(? = 0) is the energy transferred in the sudden limit. Although the 
computation of R(?) requires a solution of the detailed equations of motion, 
the outcome for atom-molecule collisions appears approximately as 

R(C) e"^ (2.7) 

Hence the efficiency of energy transfer can be expected to decrease exponen- 
tially as the encounter becomes more adiabatic. This feature will manifest 
itself in the following chapters by the use of semilogarithmic plots for all 
quantities related to R(C) when plotted as functions of the variables con- 
tained in 5* 

2.3 Theoretical Methods for Modeling Collision Dynamics 

The preceding discussion made use of simple conceptual arguments to char- 
acterize the collisional transfer of energy, but an estimate of the amount of 
energy transferred can be obtained only from a detailed solution of the equa- 
tions of motion. The motion is customarily described in the literature using 
one of three levels of quantization: (1) a fully classical treatment in which 

quantization is imposed artificially on the internal energy of the molecule 
before (and sometimes after) the collision, (2) a semiclassical approximation 
in which the path of the incident particle is obtained classically but the 
molecule response is handled quantum mechanically, and (3) a fully quantum- 
mechanical formulation in which all members of the system are represented in a 
quantum-mechanical wave equation. (A fourth intermediate level of quantization 
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might also be given as the semiquantal method developed by Miller^^ and 
Marcus^^ in which the phase distortions of the motion are obtained classically 
but subsequently treated quantum mechanically. The validity of the semi- 
quantal approximation in treating vibrational energy transfer is still a topic 
for study, but it appears to give satisfactory results for the few examples 
examined^^*^®”*^^ and has several advantages worthy of consideration. At the 
time this study was begun, however, the implications of the semiquantal 
approximation were not clearly established, precluding its further consldera- 
tlon here.) To choose the most suitable theoretical method for meeting the 
objectives of this study, we now briefly consider the general features of each 
of the three principal methods of approach. 

2.3.1 Classical Collision Theories 

A large number of fully classical calculations for collisions involving 
the vibrational and rotational motion of a diatomic molecule have been carried 
out in recent times . Modern results have shown that, when the 
oscillator zero-point energy is included (a quantum limit to the minimum vibra- 
tional energy), classical calculations for the total transfer of vibrational 
energy reproduce the equivalent quantal predictions quite well. This is not 
too surprising if one notes that the transfer of vibrational energy does not 
depend on any pure quantum effects such as tunneling or wave interference. 
However, a kind of quantum effect is ignored in a classical treatment when a 
continuum of energy is transferred to the molecule without restriction to 
discrete quantum increments. Collision energies in a moderate thermal range 
excite only a few vibrational quanta so that the partitioning of energy into 
widely separated quantized levels, excluded in a classical treatment, may have 
some influence on the molecular motion. The validity of a classical 
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description that ignores these effects is not clear. For rotational motion, 
quantum selection rules impose clearly identified limitations on the path of 
energy transfer throughout the Internal states of the molecule. The effects 
of these limitations on the total energy transferred are not very severe in 
molecules where the rotational energy spacing is small and continuum-like but, 
again, the criteria for treating the rotational motion classically are not 
obvious. Finally, a classical description of the collision d 3 mamlcs reveals 
only the total energy transferred to the molecule, but it provides no rigorous 
description of the manner in which the energy is partitioned among quantized 
internal states. Since the objectives of this study pertain specifically to 
the rate of energy transfer to individual internal states, a classical colli- 
sion model would require considerable interpretation to produce such results. 
Hence the use of a fully classical description of the collision dynamics does 
not appear to be suitable for this study and no further consideration has been 
given to it. 

2.3.2 Quantum- Mechanical Theories 

Fully quantum-mechanical calculations of vibrational and rotational 
energy transfer have also been abundant in recent literature. 

As the opposite extreme to classical treatments, they contain a complete 
description of the energy deposition and provide an exact basis for comparison 
with more approximate methods. The difficulties associated with a full quan- 
tal formulation lie mainly in the mathematical and numerical requirements to 
obtain a solution. The radial motion of the system is governed by a linear 
second-order differential equation with at least two independent variables. 

A numerical solution involves matrix manipulation and quadrature integration. 
While general methods for dealing with these numerical aspects have 
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become highly developed, their incorporation into a complete algorithm for the 
collision process is not an inviting task* Even putting that inconvenience 
aside, the experience of others^^”^^ has shown that the computing time neces- 
sary to reach a complete solution varies as the ovbe of the number of coupled 
molecular internal states included in the calculation* For calculations 
involving both vibrational and rotational states, this cubic dependence has 
been the primary factor restricting comprehensive studies of the energy trans- 
fer for all but a few special molecules like H2- To impose similar limita- 
tions on this study would yield results little different from previous work. 

Several approximations to the quantal formulation have been devised which 
allow solutions to be obtained with less effort. The recently developed seml- 
quantal method^^>^® is one example, but a much older and more easily applied 
approximation is the “distorted wave** approach first proposed by Zener^»^^ and 
later reformulated by Jackson and Mott.^^^ We briefly mention the distorted 
wave approximation here to show that it too is not the best choice for the 
purposes of this study even though it retains much of the **exactness** of a 
full quantal solution. The approximation has been applied in the past both to 
collinear collision models^^ and to three-dimensional rotational models. 
Generally, it is a perturbation method that may be carried to arbitrary order, 
but its greatest advantage is realized by retaining just the first-order term. 
As with all first-order perturbation solutions, the results are accurate only 
when transition probabilities and interactions between nonadjacent states are 
small. This study of excited-state transitions, in which single- quantum 
transitions are expected to be large and multiple-quantum transitions to be 
important, would then require at least a second-order theory for accuracy. 

The distorted wave approximation is therefore not an attractive choice here* 
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2.3.3 Semiclassical Collision Model 


A semiclassical or "Impact-parameter” method for treating the collision 
dynamics has several attractive features that promote its selection as the 
primary theoretical approach in this study. For one, it retains all the 
quantum-mechanical aspects of energy transfer within the internal states of 
the molecule that are absent in a classical treatment while avoiding most of 
the mathematical difficulties associated with a full quantal solution. In the 
semiclassical approximation, the Schrodinger equation describing the collision 
may be reduced to three, coupled, time-dependent , first-order, linear, differ- 
ential equations — one describing the molecular wave- function dynamics and two 
classical trajectory equations for the relative motion of the incident par- 
ticle. When the trajectory and wave equations are uncoupled, the trajectory 
equations may be exactly integrated analytically for collinear encounters and 
approximately for spherically symmetric interactions. Clearly, these reduc- 
tions relax the numerical requirements considerably. Furthermore, the reduc- 
tion to obtain first-order differential equations affords a second, and per- 
haps more significant, advantage to the semiclassical formulation because the 
computation time then varies with the square of the number of coupled channels, 
making feasible the calculation of much larger sets than might be considered 
with a full quantal method. 

From another viewpoint, the ability to separate the motion of each colli- 
sion partner in a semiclassical formulation allows analytic solutions to be 
obtained for the final state of the molecule if some further approximations 
are made. The usefulness of these more approximate analytic solutions was 
discussed in chapter 1, but their accuracy requires validation. The valida- 
tion is done most convincingly by comparing the analytic predictions with 
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numerical solutions in which the trajectories are also obtained classically. 

By choosing a general semiclassical approach, such comparisons may be included 
as part of the results of this study. However, before a classical description 
of the incident particle motion can be justified, several restrictive criteria 
must be met. The following section delineates the pertinent criteria that 
apply to vibratiotially inelastic collisions. 

2.4 Basic Criteria for the Semiclassical Approximation 

A detailed derivation of the semiclassical equations of motion is pre- 
sented in appendix B. In this section, we shall discuss only the general cri- 
teria necessary to justify a semiclassical formulation. 

The conditions for the validity of semiclassical theory have been exam- 
ined for numerous applications in the past by many authors; but recently, 

Delos et reported a careful and detailed reexamination that revealed 

some of the implications and restrictions of a semiclassical collision model 
in much greater depth. They show that the classical trajectory equations may 
be obtained in two fundamentally separate ways: one based on a classical 

wave-packet description involving the correspondence principle, and one that 
makes no reference to a conventional classical picture but is based on an 
extension of the usual WKB approximation. Their work was motivated by the 
observation that semiclassical models work well even at collision energies too 
low to justify a localized wave-packet description of the incident particle. 
Indeed, they found, by comparison with the second method for obtaining the 
same classical equations, that criteria based on the localized wave-packet 
concept were overly restrictive. Specifically, in the classical wave-packet 
picture, the correspondence principle leads to the classical trajectory equa- 
tions only if 
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(X/L)1/2 « 1 (2.8) 

where A is the de Broglie wavelength and L is the range of the collision 
interaction. On the other hand, the WKB approach yields the classical trajec- 
tory equations by requiring the weaker condition 

X/L « 1 (2.9) 

Thus, Delos et aZ. conclude that, even when a classical picture involving the 
correspondence principle is not justified, the classical trajectory equations 
may still produce a reasonable description of the final molecular state. They 
suggest that the classical trajectory equations are implicitly more fundamen- 
tal to the collision dynamics than just in the correspondence limit. The com- 
parisons of semiclassical and quantum-mechanical collision models for harmonic 
oscillators by Rapp and Kassal^^ and for enharmonic oscillators in chapter 4 
support these conclusions. For example, in chapter 4, typical threshold col- 
lision energies for a single vibrational quantum transition correspond to 
X/L ^ 1, which violates even equation (2.9), and yet the semiclassical and 
full quantal predictions appear identical. Delos et at. warn that by the same 
rule, however, a classical interpretation of the computed oscillator dynamics 
during intermediate times in the collision should not be given unwarranted 
value for conditions outside equation (2.8). Since all comparisons to date 
have examined only the predictions of final oscillator states, no conclusions 
can be inferred about intermediate times. These warnings suggest some caution 
when coupling the classical trajectory to the molecular dynamics, for example. 

Delos et at. continue by showing that equation (2.9) is not the only cri- 
terion necessary to validate a semiclassical theory. They state two further 
requirements that, in combination, demand that the elastic collision trajec- 
tory for all channels (internal molecular states) be approximately the same. 
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These criteria are most restrictive when several potential energy surfaces are 
considered and when the system is allowed to cross from one surface to another; 
but the applications here, where reactive collisions and electronic transi- 
tions are not considered, require only a single Interaction surface. These 
additional criteria then reduce to the stipulation that the difference in 
diagonal matrix elements defined by the instantaneous state of the molecule 
and by the interaction potential be small. Thus, if is the Instantaneous 

diagonal matrix element for state k, the criteria requiring similar elastic 
trajectories may be expressed by 


V' - V* 
nn kk 

v’ + V' 
nn kk 


« 1 


( 2 . 10 ) 


With simple exponential interactions of the type most commonly used (see ch. 3 
for examples) , varies with k primarily as a result of vibrational 

anharmonicity . When the molecule is treated as a harmonic oscillator and the 
Interaction potential is also linear in the oscillator coordinate, is 

identical for all k and equation (2.10) is satisfied precisely. Conversely, 
we can expect the results of a semiclassical treatment of anharmonic oscil- 
lators with nonlinear interactions to compare differently with their full 
quantal counterparts than found in similar comparisons using harmonic oscil- 
lators. On the other hand, varies only weakly with k for most molecules 

because the vibrational anharmonicity is generally small. Hence a semiclassi- 
cal model of anharmonic oscillators can still be expected to retain a large 
measure of accuracy provided the differences remain small. 


With the preceding criteria in mind, we shall adopt a semiclassical for- 


mulation throughout this study. However, since the criteria given by equa- 
tions (2.9) and (2.10) are not quantitatively explicit in establishing the 
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range of collision parameters and molecular properties that may be suitably 
applied in such a method, much more explicit statements regarding the validity 
of a semiclassical collision model will be obtained from the results in chap- 
ter 4. Before proceeding with numerical solutions, however, the intramolecu- 
lar potential determining the molecular dynamics and the collision Interaction 
potential determining the perturbing forces on the molecules must first be 
modeled. Chapter 3 describes these potential models and the considerations 
leading to them. 
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CHAPTER 3 


INTRAMOLECULAR AND COLL IS ION- INTERACTION POTENTIALS 

A rudimentary aspect of any collision model is the description of forces 
acting between elements of the system. In this study, three separate nuclei 
constitute the system and, in principle, the forces acting on any one of them 
depend on the relative positions of all three. In practice, however, a self- 
consistent potential surface is rarely known, except for the simplest systems. 
Instead, the potential surface is usually constructed in a semlempirical 
manner using simpler concepts. To that end, we follow the conventional tech- 
nique of considering only independent pairwise interactions. The potential 
surface is then separable into three additive components, each dependent on 
only the distance between two nuclei. Two of the components include the inci- 
dent particle as one nucleus and are classified here as "collision- interaction 
potentials." The third component is between the two bound molecular nuclei 
and termed here the "intramolecular potential." Both types are modeled indi- 
vidually below. 

3.1 Intramolecular Potential 

Since only pairwise interactions are considered, the intramolecular 
potential is independent of the disturbance to the molecule brought by a col- 
lision. The potential model is therefore based on the spectral properties of 
an undisturbed molecule. 

3.1.1 Vibrational Anharmonicity 

The importance of vibrational anharmonicity in the intramolecular motion 
has been emphasized several times in previous chapters and its inclusion has 
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been stated as a basic feature of this study, Anharmonicity will manifest 
itself in the molecular model by the appearance of second- and higher-order 
terms in the expression for vibrational eigenenergies . A minimum example is 
then 

Ev/h = + I) - (3.1) 

The first term in equation (3.1) is the familiar result obtained for a harmonic 
oscillator potential of the form 

Vo(r) = j Uo03e(r - (3.2) 

where is the reduced mass of the oscillator, r is the internuclear 

separation, and r^ is the equilibrium separation. Numerous potential func- 
tions will produce an anharmonic term identifiable with the second term in 
equation (3.1). In fact, any potential function with a higher-order depen- 
dence on r - r^ than equation (3.2) will do. However, two highly desirable 
additional features of the potential function are (1) that it realistically 
represent a real molecule for all separations to dissociation and (2) that it 
be of an analytically convenient form. Again, several potential functions 
have been proposed that fit these requirements (see ref. 71, ch. 5, for 
examples). One example that has received perhaps the greatest attention since 
its conception is that proposed by Morse^^ in the form 

V„(r) = (3.3) 

In the absence of molecular rotation, Morse obtained eigenenergies approxi- 
mated by equation (3.1) plus higher-order terms that are clearly negligible 
for all diatomic molecules. The potential constants are related to the oscil- 
lator frequency and anharmonicity by 

Dq = ha)^^/4a)^Xg (3.4) 
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and 


a = ( 2 yQ 0 )gX 3 /h) (3.5) 
Morse went further to show that if the steady-state wave function is separated 
according to 

Rv(r) 

<Kr,9,4>) T" 

where is a normalized spherical harmonic function in the polar 
angles 0 and <J), then an analytic solution to the radial wave equation is 
obtained for integer values of v given by 


where 


R^(r) = (z) 

z = k’ 

= [ab'r(v + l)/r(k' - v)]^/2 


and L^(z) is the Laguerre polynomial 


(3.7) 


T^' / = V (-1)”^ r(k’ - v) 

/ j m! (v - m)! T(k' - 2v + m) 


r(k* - v) ja 


m=0 


(3.8) 


The parameters k' and b’ are defined in terms of the oscillator frequencies 
by 

k’ = 2(2 pqDq) = (Og/ojgXg (3.9) 


and 


b' = k' - 2v - 1 (3.10) 

The Morse potential is not as accurate as other more recent models (e.g. , 
the function proposed by Hulburt and Hirschfelder^®) , but its analytic conven- 
ience has made its appearance in the literature seemingly second in popularity 
only to the harmonic oscillator. Figure 3.1 illustrates a comparison of the 
Morse function for the hydrogen molecule with a more accurate numerical 
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Figure 3.1- Intramolecular potential functions for H 2 (taken from ref. 74, 

based on ref. 75). 

representation obtained by the Rydberg-Klein-Rees method^^*^^ of analyzing 
spectroscopic information. Hydrogen is an example in which the Morse potential 
appears at its worst because of the high degree of anharmonicity associated 
with Figure 3.1 suggests that the Morse potential model should reproduce 

the molecular properties accurately for vibrational energies at least up to 
^ However, a comprehensive comparison for H 2 is difficult because 

high- lying vibrational energies in are not available. They are available 
for CO, however, where they have recently been measured up to v = 37 using 
laser-spectroscopy techniques . Figure 3.2 compares vibrational energies 
of CO, obtained from the laser measurements and defined by terms up to sixth 
order, with the second-order Morse expression. As expected, values of 

< Dq/ 2 are reproduced by the Morse function with reasonable accuracy. We 
should note at this point, however, that eigenenergies are not as sensitive to 
variations in the intramolecular potential as are the wave functions or their 
overlap integrals used in collision theory. Thus the measure of accuracy 
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Figure 3, 2.- Vibration eigenenergies of CO (Dq = 11 ,108 eV). 

suggested by figures 3.1 and 3.2 is not entirely representative of the corre- 
sponding accuracy given to a collislonal model. Nevertheless, in view of the 
other inaccuracies inherent in the collision model, the Morse function has 
been adopted here as a sufficiently realistic Intramolecular potential model. 

3.1.2 Vibration-Rotation Coupling 

A remaining question regarding the molecular model pertains to its rota- 
tional characteristics. In a classical picture, the vibrational and rota- 
tional motion of the molecule are clearly coupled. A quantum-mechanical 
description of the molecule will produce corresponding coupling terms in the 
wave functions and eigenenergies associated with vibration and rotation. The 
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task is to determine the degree of coupling that must be retained for the 
purposes of this study. 

Pekeris^^ has solved the steady-state radial wave equation for a Morse 
potential with the rotational terms included. His results can be reformulated 
to better show the influence of vibration- rotation coupling in the following 
way: to second order, the radial wave function obtained by Pekeris is similar 

in form to equation (3.7). Part of the difference lies in the parameters, b’ 
and k’ , which depend now on the rotational quantum number £. Denoting the 
il-dependent parameters as and k^, they are related to the values b’ and 
k* given previously for a nonrotating molecule by 

kji = k' (1 + (3.11) 


A - < 

1 + c„ 


- 2v - 1 


( 3 . 12 ) 



and 


= £(£ + l)2h2/(y^are3Do) 


(3.14) 


In this extended notation, the radial wave function becomes 

R^j^(r) = (3.15) 

in which z = k e as before and = [ab£F(v + l)/r(k£ - v)]^/2, 

equation (3.15) shows, the wave function is distorted by in the coordi- 

nate 2 . 

To evaluate the effect of vibration-rotation coupling on the molecular 
model, equation (3.15) must be used in place of equation (3.7) when computing 
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the appropriate overlap integrals appearing in the collision theory formula- 
tion. Obviously, that kind of evaluation can be made only in retrospect after 
the collision theory has been formulated in detail. However, at this point, 
we are at least able to indicate the extent of the coupling by computing the 
magnitudes of the corrections to b* and k’ and the distortion of the wave 
function introduced by A^. Choosing the properties of CO as an example and 
defining £ < 40 as the range of interest, we find that, for £ = 40, 
kj^/k’ = 0.997, b^ = 0.96k - 2v - 1, = 1.05 

The overlap integrals are sensitive to small changes in the wave function so 
the effect of these small corrections is not clear; but one can see that the 
influence of rotational coupling is not a major aspect of the vibrational 
motion. 

The initial choice of a method for describing the molecular properties in 
this study was guided by a less sensitive but readily evaluated measure of the 
vibration-rotation coupling. The eigenenergy terms associated with vibration- 
rotation coupling were compared with those resulting from vibrational anhar- 
monicity. The objective was to decide if vibration-rotation coupling should 
be included to maintain a consistent degree of accuracy in the eigenenergy 
expression. The energy expression corresponding to the second-order radial 
solution given by equation (3.15) was shown by Pekeris^^ to be 


( 

l\ / 

i\2 r / 

1 \ 

= ‘^e\^v + 

2 ) - 

h^j + £(£ + l)Bg - £(£ + 1) 

' + |j + Dg£(£+1) 


(3.16) 


We can then ask for what values of Si is 

,2 
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Again, using spectroscopic constants for CO, the rotational coupling terms 
equal or exceed the enharmonic terms when £(£ + 1) > 623 [v + (1/2)]. The fol- 
lowing table lists some sample values. 

TABLE 3.1.- ROTATIONAL QUANTUM STATES WITH SECOND-ORDER CORRECTIONS 
COMPARABLE TO THE ANHARMONIC CORRECTION FOR CO 


V 

z 

(v + i) 

0 

17 

0.006 

1 

30 

.009 

2 

39 

.012 

5 

58 

.031 

10 

80 

.062 

20 

113 

.124 


Thus, for £ > 40, the second-order rotational terms are comparable only with 
enharmonic terms for the first few vibrational states. The third entry in 
table 3.1 is a measure of the anharmonicity correction that may be compared 
with unity. For the first few vibrational states, it is much less than unity 
and we may therefore conclude that, while a completely consistent molecular 
model should contain vibration-rotation coupling terms, their influence on the 
molecular dynamics is not expected to be large in any case. At the time the 
molecular model was formulated for this study, the small amount of added com- 
plexity introduced by the vibration-rotation coupling terms appeared to be 
greater and a Morse-oscillator/rigid-rotor description was chosen to ensure 
analytic progress. The eigenenergies are then expressed by 

= CD^^v “ ^e^e (3.17) 

and the steady-state radial wave function is given by equation (3.7). In 
retrospect, the coupled expressions of Pekeris given by equations (3.15) 
and (3.16) would have required additional computational effort, but they do 
not increase the complexity of the formulation substantially. 
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3. 2 Collision- Interaction Potential 


The outcome of any scattering event depends strongly on at least some 
features of the interaction potential; yet the shape and magnitude of interac- 
tion potentials are poorly known for all but a relatively few simple cases* 

The potentials between elastically scattered atoms are generally well estab- 
lished from both theory and atomic-beam experiments, but the interactions 
influencing inelastic collisions involving diatomic molecules are still an 
active subject for research and computation. The topic has been discussed 
extensively in relation to the transport properties of gases®^ and, more 
recently, in relation to scattering events.®^ To circumvent the complexity of 
the subject, we develop in this section an empirical model of the interaction 
forces based on the general nature of the interaction. 

The nature of the interaction forces depends greatly on the modes of 
energy transfer and on the internal energy states that participate in the 
collision dynamics. Here we are interested only in nonreactive interactions 
between collision partners in their ground electronic states. Even then, ah 
'in'it'Lo calculations for three-body systems of the type considered are diffi- 
cult and still incomplete for even the simplest system, H 2 -H (e.g. , ref. 82b). 
The Born-Oppenheimer separation of electronic and nuclear motion is generally 
used, but the complexity associated with electronic coupling between charge 
clouds of three nuclei has limited present accomplishments. A few cases have 
been obtained using ’’self-consistent field theory” where the number of elec- 
trons is a minimum (H^-H (ref. 82), H 2 "He (refs. 83, 84), H 2 “Li“^ (ref. 85)), 
but similar calculations for heavier nuclei are either less rigorous (HF-HF 
(ref. 86)) or not available. Experimental determinations appear to be 
equally difficult. The interaction potential cannot be measured directly but 
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must be implied from some colllsionally dependent observable property in which 
the potential function is only implicitly contained. The uniqueness of the 
potential thus determined is often in serious doubt and the accuracy of its 
details are consequently diminished. Generally, the collision- interaction 
potential, however obtained, appears to be the greatest source of uncertainty 
in the calculation of vibrational inelastic collision dynamics for most mole- 
cules larger than H 2 - 

Despite our inability to accurately define the interaction potential, we 
can at least describe its qualitative features with some confidence. To sim- 
plify the description, the potential may be separated into two major compo- 
nents: an average spherical component that depends only on the separation of 

the molecular mass center and the incident particle and an anisotropic compo- 
nent that accounts for variations with the molecular orientation relative to 
the direction of the incident particle location. Figure 3.3 is a qualitative 
sketch of both components. 



Figure 3.3- Characteristics of nonreactive intermolecular potentials. 


46 



3. 2,1 Spherical Features 

The spherical component may be further divided, though somewhat arbitrar- 
ily, into several types of forces. They are commonly referred to as (a) short- 
range repulsive forces that occur principally when two nuclei are close enough 
for their electronic charge clouds to overlap, (b) long-range attractive 
forces that result from large separations in which the charge distributions of 
each nucleus are independently distorted, and (c) intermediate forces that 
simply refer to the transition region between the preceding extremes. 

The short-range forces do not lend themselves easily to theoretical sim- 
plification but, fortunately, the collision dynamics of vibrationally inelastic 
events are not particularly sensitive to their precise shape. Usually, the 
short-range force gradient at closest approach is the single most important 
feature to the collision dynamics. Generally, theoretical attempts to model 
the short-range forces have led to sums of terms like 

VgR *=» A(x) (3.18) 

where x is the separation distance, A(x) is slowly varying in x, and L 
is a constant range parameter. Commonly, A is also taken as a constant and 
effective values of A and L are deduced from experiment by simple comparison 
with predicted results.®^ The results of ab-initio calculations have also 
been reproduced surprisingly well by equation (3.18) in interactions between 
two atoms, but combination rules for cases where three interacting nuclei are 
present are still a topic of discussion. 

Long-range attractive forces, on the other hand, are more easily approxi- 
mated from well-founded physical arguments. At separations large compared to 
the intranuclear separation of the molecule, the molecule appears as an inde- 
pendent multipole with some polarizability. The resulting multipole 
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interaction forces are the same both in quantum-mechanical and classical 
treatments with the general result appearing as a series of terms in the form 


''LR 



(3.19) 


n 

Often the dominant terms are induction forces (dipole- induced dipole) and 
London dispersion forces (similar to dipole- induced dipole) , which both vary 
as 


V »» -C/x® 


(3.20) 


The single term represented by equation (3.20) has commonly been used to 
approximate the influence of long-range forces. 

Numerous empirical representations for the entire potential have emerged 
from considerations similar to those just described. An immediate choice 
would appear to be the Buckingham potential 

Vg = A - C/x® (3.21) 

but this form has the unrealistic property of reaching a maximum for small x 
and becoming infinitely negative as x ^ 0. A more realistic formula that 
emphasizes the same long-range force dependence is the well-known Lennard- 
Jones 12-6 potential: 



The zero-potential separation, explicitly appearing in equation (3.22), 

provides a useful measure of the equivalent "hard-sphere” radius that will be 
required in the collinear collision model application described in chapter 5. 
Similarly, the well depth, D, is also an explicit parameter. The Lennard-Jones 
potential gained early popularity in the analysis of transport properties in 
gases and its effective constants, and D, have been evaluated for many gas 
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mixtures using viscosity data and virial coefficients. Unfortunately, the 
mathematical form of equation (3.22) is inconvenient for the calculation of 
interaction overlap integrals in a collision model where the oscillator is 
treated quantum mechanically. For that reason, we shall seek a more conven- 
ient representation, but we can use the Lennard- Jones potential as a basis for 
comparison. 

A potential representation consisting of exponential terms is particularly 
convenient in the mathematics of enharmonic oscillators. One such representa- 
tion is the Morse potential 

v„ . D - 2D (3 23 , 

where the separation at the potential minimum, x^, is related to the Lennard- 
Jones zero-potential separation by x^ = 2^/^p^. Equations (3.22) and (3.23) 
can then be made to yield similar potentials in the region of the potential 
well by a suitable choice of L. However, the shape of the potential well is 
usually not important to the collision dynamics at energies sufficient to 
induce vibrational transitions. The threshold energies for vibrational tran- 
sitions are near E ~ - 0.2 to 0.5 eV, while apparent well depths for 

many-electron molecules like N 2 , CO, O 2 , etc., are typically D - 0.01 eV. 

Thus, E/D >> 1 in most cases of interest and the influence of the potential 
well shape is negligible. Instead, the choice of range parameter, L, will be 
dictated by the greater necessity for properly matching the potential gradient 
where V - E. This latter requirement for gradient matching suggests that, 
unless the potential shape is correct in the region of short-range forces, 
the effective range L will depend on E. Later comparisons of theoretical 
and experimental rate coefficients will show that such behavior is obtained 
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and the thermal range of applicability for a given set of potential parameters 
is correspondingly limited. 

Arguments similar to those regarding the potential well region allow us 
to further simplify the interaction potential representation (for the purposes 
of this study) by also neglecting the long-range forces. The role of the 
long-range attractive forces is basically twofold during a collision event — 
they accelerate the incoming particle by an amount related to D/E and they 
induce early transitions between levels whose energy spacings are the order of 
D. Only rotational transitions within a given vibrational state are affected 
in the latter role. In the first role, the effect of long-range forces will 
be apparent only at energies very near threshold and hence only at very low 
temperatures. However, this study is not directed at low- temperature applica- 
tions in particular. In their second role, long-range forces would be impor- 
tant if we were interested in pure rotation-translation energy transfer or in 
the detailed final rotational state of the molecule after a collision involv- 
ing V-R-T energy transfer. We shall find, however, that the net rate of vibra- 
tional energy transfer is not particularly sensitive to the rotational dynam- 
ics during the early or late stages of the encounter. The rotational state of 
the molecule having the greatest influence on its final vibrational condition 
will be the rotational state occupation occurring at the time of closest 
approach when the short-range forces are dominant. Thus, we can justifiably 
neglect the long-range forces entirely for these purposes and adopt the simple 
and analytically convenient repulsive potential 

V = A (3.24) 

The effects of the potential well on vibrational energy transfer have 
been investigated in greater detail by others. Their findings support 
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the conclusion that the well may be neglected when E/D or hoj^/D » 1. 

There are examples, of course, where the well depth is larger (e,g., for 
H^-Li"^, ho)g/D-0,2). For those cases, the entire interaction potential must 
be accurately represented. 

3.2.2 Anisotropic Features 

The anisotropic potential component is responsible for rotational distur- 
bances to the molecule. Figure 3.3(b) illustrates the typical magnitude of 
anisotropy effective during the collision. The magnitude is gauged by noting 
that the equilibrium separation of most diatomic molecules is near r^ == 0.1 nm, 
while the distance of closest approach during a collision will be only 
slightly less than the zero-potential radius, p^. For most common molecule- 
atom interactions, 0. 3 nm. Figure 3.3(b) is drawn for the ratio 

Po/^e = 3 with a radius p^ centered on each molecular nucleus. The equi- 
potential appears mostly spherical with relatively small aspherical components. 
The small anisotropic terms of most calculated interaction potentials confirm 
these observations. 

When rotational motion is considered, it is most conveniently described 
by a coordinate system containing an angle, 6, that defines the rotation of 
the molecular axis relative to the. position of the incident particle. Fig- 
ure 3.4 illustrates such a coordinate system. The subsequent mathematics then 
appear in a convenient format if the potential is expressed as a series 
expansion of Legendre polynomials, Pj, in the form 

V(x,r,6) = ^ Uj(x,r) Pj(cos 6) (3.25) 

In most cases, db-invt'lo potentials are represented by equation (3.25) using 
only two or three terms. 
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Figure 3.4~ Three-dimensional encounter nomenclature viewed in the time- 

dependent plane defined by the location of the incident particle A and 
intramolecular axis BC. 

The typically small anisotropy of most interactions makes two simplifying 
assumptions reasonable that greatly reduce the complexity of a three- 
dimensional collision model. One is operative if the incident particle motion 
is treated classically. An approximate trajectory may then be obtained using 
only the spherical component of the potential. Such an assumption neglects 
out-of-plane deflections and reduces the translational motion to two dimen- 
sions. Averaging the initial configuration is then greatly simplified. A 
second simplification, in keeping with our original notions, is that the aniso- 
tropic potential can be approximated by two additive spherical potentials, 
each centered on a molecular nucleus. In concept, the use of pairwise poten- 
tials ignores the second-order mutual interaction between nuclei and it also 
omits the shielding or shadowing of one nucleus by the other as viewed from 
the position of the incident particle. However, if the anisotropy is small, 
these second-order corrections will be smaller. In addition, the interaction 
potential decreases rapidly with internuclear separation so that for most 
angles, 6, the incident particle will be interacting mainly with only one 
molecular nucleus at a time. Thus, shielding effects should be relatively 
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unimportant. For the purposes of this study, the three-dimensional interac- 
tion potential is therefore simulated by two noninteracting short-range poten- 
tials, each given by equation (3.24) and each centered on a molecular nucleus. 
The final step in representing the potential is to express it in the form given 
by equation (3.25). 


3.2.3 Three-Dimensional Interaction Potential Model 

Part of the nomenclature to describe a three-dimensional encounter is 
defined by figure 3.4. The figure lies in a time-dependent plane containing 
both molecular nuclei and the instantaneous position of the incident particle. 
An additive repulsive force between each nucleus will then produce the inter- 
action potential 


V(x^,,x^) = A(e + e (3.26) 

where A and L are considered identical for both interacting pairs. The 
internuclear separations may be expressed in terms of mass-centered variables 
as 


- f 

■4 


cos 6 + 


1/2 


1 + 2(1 - Y) cos 6 + (1 - Y)2('i 


1/2 


(3.27) 


where Y = [m^/ (m|j + molecular mass ratio and m^ > m^. For the 

potential to appear nearly spherical, r/x must remain small over the entire 
trajectory. We have argued that it does remain less than unity since its 
largest value at closest approach is only r/x - 1/3. Equation (3.27) may 
then be expanded in a uniformly convergent power series of r/x, giving to 
first-order, 
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(3.28) 


(3.29) 


X|, = x^l - Y ^ cos 6 + . . .^ 

Xj, = x|l + (1 - Y) cos 6 + . . .J 

The corresponding first-order representation of equation (3.26) is then 

V(S,r,«) = A «] 

Figure 3.5 demonstrates that the first-order potential, equation (3.29), is a 
reasonable approximation of the additive spherical potential, equation (3.26), 
for the typical range of parameters used here. 

The potential model given by equation (3.29) is now in a form that can be 
conveniently handled in the framework of a three-dimensional collision model. 
It may be transformed to a form equivalent to equation (3.25) by use of the 
expansion 


±z cos 6 


= J (2J + iKil)*^ Ij+l/2^^^ Pj(cos 6) (3.30) 



Figure 3.5.- Variation of a pairwise repulsive potential and its first-order 
approximation with orientation angle. Calculations were done for 
r/L = 5. Solid curves are the infinite-order variation from equa- 
tion (3.26); dashed curves are the first-order approximation, 
equation (3.29). 
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(see ref. 73, p. 445) where Ij+ 1/2 is the modified spherical Bessel function 
of the first kind and Pj is a Legendre polynomial as before. By reducing 
the notation with the definition 

ij(z) = /^72? <3- 31) 

equation (3.29) Is transformed to the series representation: 

V(x,r,6) = A e"^/^ J <23 + l){^j[Yr/L] + - Y)r/L]| Pj(cos d) 

(3.32) 

As expected, '^j(z) decreases rapidly with J for typical values of z so 
that only a few terms contribute to the summation. The correspondence between 
equations (3.25) and (3.32) is obvious. 

3.2.4 Collinear Interaction Potential Model 

The discussions of previous sections have noted that a one- dimensional 
collinear collision model will be useful in obtaining analytic descriptions of 
the rate of vibrational energy transfer. A brief development of the required 
interaction potential expression is therefore included here. 

The most effective collinear collision to induce vibrational motion will 
be an encounter with the lightest nucleus in the molecule. Collinear encoun- 
ters are therefore limited to those where 6=0. Equation (3.28) then 
reduces to 

V(i,r) - A (3.33) 

Typically, r/L - 5 while Y is always less than unity so that the second 
term in equation (3.33) is always much smaller than the first. The interac- 
tion potential for one- dimensional collinear collisions is therefore taken to 
be 
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V(i,r) - A (3.34) 

This form of the potential or a linearized version of It have been used for all 
colllnear vibrational energy transfer theorles.*^^***®***® It Is applied In the 
following chapter to study the applicability of a semlclasslcal collision 
theory to enharmonic oscillators. 
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CHAPTER 4 


A COMPARATIVE EVALUATION OF THE SEMICLASSICAL APPROXIMATION 

We have shown in chapter 2 that, while a semiclassical collision model 
offers several advantages in a study of vibrational energy transfer, the cri- 
teria for its application are not explicit- Hence they offer marginal help in 
evaluating the validity of the semiclassical approximation for a specific 
set of conditions- We have also observed that the criteria, expressed by 
equations (2.9) and (2.10), suggest that the validity of a semiclassical model 
will be influenced by the oscillator anharmonicity . However, previous com- 
parative evaluations of semiclassical theories^ ^ have been only for harmonic 
oscillators. Furthermore, the results, although favorable, have not been 
entirely conclusive because, in addition to the absence of anharmonicity, the 
semiclassical formulation is typically only one of several approximations con- 
tained in the comparisons, while all the corrections known to improve the 
semiclassical predictions are not always included. The primary purpose of 
this chapter is therefore to compare the vibrational transition probability 
predictions from a semiclassical model for enharmonic oscillators with a 
comprehensive set of solutions from an equivalent, fully quantum mechanical, 
collision model- The onVy difference in the two models is the treatment of 
the incident particle motion. 

The physically incomplete nature of a semiclassical treatment requires 
some interpretation and correction, however, to facilitate its correlation 
with more exact collision models. Hence, in this chapter, we shall also deal 
explicitly with the corrective aspects of a semiclassical treatment. For 
example, a well-known weakness of semiclassical theories is the inherent lack 
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of energy conservation. Several methods of compensation have been suggested 
that aim at interpreting either the classical trajectory energy^® or veloc- 
ity^^ in terms of corresponding values averaged over the collision. Compari- 
sons described here between these semiclassical predictions and the exact 
quantum-mechanical calculations show that, while such an Interpretation is 
necessary to correct the semiclassical predictions, the results are insensi- 
tive to the choice of method in the energy range of practical interest. 

Regardless of the corrections for energy conservation, the conventional 
semiclassical treatment will be shown to fail badly in some cases. The 
failures appear in the form of anomalous resonances that occur only in anhar- 
monic oscillator models and are caused by an incomplete account of the oscil- 
lator compression and recoil during impact. Within the usual semiclassical 
framework, the classical trajectory is computed, assuming that the oscillator 
remains in a pure eigenstate having a fixed average separation of its nuclei. 
In reality, the oscillator is compressed by the impact and enters a mixed- 
state condition in which the average internuclear separation oscillates with 
frequency components from each of the excited states. To include this 
behavior in the semiclassical theory is not equivalent to conserving energy, 
but it has the effect of introducing an oscillator "feedback" on the classical 
trajectory. The effect can change the entire nature of the results in some 
cases. Extremely heteronuclear or enharmonic molecules, such as the hydrogen 
halides, will be shown as members of the class strongly affected. 

We also simplify the collision geometry used here by confining it to 
colllnear encounters. Direct comparisons with the fully quantum mechanical 
results of reference 68 are thus made possible. Calculations have been 
presented in the literature of more realistic three-dimensional encounters, 
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but mainly for harmonic oscillators initially in the ground state. They have 
also been made in chapter 6 for anharmonic oscillators in excited vibrational 
states. However, a three-dimensional collision geometry introduces many addi- 
tional complexities, as we shall demonstrate in chapter 6, and little would be 
gained by including it in this chapter. 

In the paragraphs to follow, a multistate semiclassical formulation 
requiring numerical solution is assembled first that Includes modifications of 
the standard treatment to account for the effect of oscillator response on the 
classical trajectory. The model is entirely equivalent to the fully quantvnn 
mechanical model in reference 68, except for the classical treatment of the 
incident particle motion. The accuracy of a first-order perturbation theory 
used by Mies^^ for anharmonic oscillators is also evaluated. As expected, 
the first-order theory is suitable only where the transition probabilities 
are small; but it must also be limited to cases where the oscillator feedback 
effects are negligible. Such cases will be shown to pertain mainly to heavy 
homonuclear molecules impacted by lighter collision partners. 

4.1 Semiclassical Model for Collinear Collisions 

A full description of the semiclassical formulation for a general colli- 
sion geometry is given in appendix B. In this section, only the results 
pertinent to collinear collisions are recalled in detail. 

The collinear collision geometry is shown in figure 4.1 for a structure- 
less particle of mass, m , impacting a diatomic heteronuclear molecular with 

a 

nuclear masses, m^ and m^. The impacted oscillator nucleus, m^, extends from 
the molecular mass center by a distance yv y where y “ ^ 

three-body center-of-mass reference frame is taken in which the relative 
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Figure 4.1.- Collinear collision geometry. 


collision speed is u. (Barred symbols identify the incident particle vari- 
ables to be evaluated classically and later interpreted as average values.) 
Except for the notation, this configuration is identical to those used by 
Mies^^j^^ and in references 45 and 68. 

The interaction potential to be used here is of the same form given in 
reference 68 , namely , 

V(*) = A 

where L and A are constants. The potential in terms of mass center and 
oscillator coordinates defined in figure 4.1 then becomes 

V(S,r) - A (4.1) 

The Hamiltonian for the three-body system is given by 

- i;r + 2^ 

where the symbols with subscript o refer to oscillator quantities and the 
other symbols denote incident particle variables. The oscillator reduced 
mass is = m^m^/ (m^ + m^) and the collision reduced mass is 
y = + m^ + m^). 
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4, 1.1 Incident Particle Motion 


The application of Ehrenfest’s theorem to the incident particle dynamics 
guides the formulation of its equation of motion with the quantum mechanically 
averaged quantities properly included. In terms of quantum mechanically 
averaged variables, the equations of motion are 

d<p> _ /9^\ _ /9y(x»r)\ 

dt \8x/ “ \ 3x / 

= /M\ . i <p> 

dt \3p/ y ^ 

By separating the total wave function according to 

H^(x,r,t) = (f)(x,t)T|^(r,t) 


and treating the incident particle classically so that 


<p> -> yu 


and 


<x> X 





the equations of motion for the incident particle combine with equation (4.1) 
to give 


y 


d^x A -x/L / yr/L\ 
^ “ L ® <f / 


where the bracket notation implies 

4'(r,t)dr (4.2a) 


For the purely repulsive exponential potential used here, the potential con- 
stant, A, influences only the distance of closest approach, a quantity of no 
direct consequence to the transition probabilities. It may be removed by a 
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transformation suggested by equating the potential energy at closest approach 
for a stationary oscillator with the initial kinetic energy, that is. 


A - i 
'^kk 


where E = (1/2) yu^ is the semiclassical relative collision energy before 
interaction, is the time- independent diagonal matrix element defined by 


and is initial oscillator stationary-state eigenfunction. Note 

that Xq is not the distance of closest approach for a nonstationary oscil- 
lator. However, if a new interaction coordinate is defined as z = x - x , 

o’ 

the incident particle motion is unaffected and may then be described by 

y § = f e~^/^R(t) (4.3a) 

where 

R(t) . (4.3b) 

The variable R(t) represents the quantum-mechanical average effect of the 
oscillator motion on the classical path. It is a measure of the oscillator 
distortion during impact and subsequent "ringing" afterward. 

To compute the trajectory classically, we must first describe the oscil- 
lator motion to obtain R(t). The usual practice at this point has been to 
consider the oscillator fixed in its initial pure eigenstate so that 
i|;(r,t) ^ Then R(t) ^ 1 for all time and the classical equation of 

motion is reduced to the equation for a constant energy trajectory: 


E -z/L 
IE2 “ L ® 


(4.4) 
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Equation (4.4) can be Integrated analytically**® so that the interaction poten- 
tial 

V(t,r) = i ^yr/L (4.5a) 

can be written explicitly in terms of time by use of the result: 

(4.5b) 

In such an approximation, the transformation parameter becomes the 

distance of closest approach. 


4.1.2 Oscillator Motion 

The unsteady motion of the oscillator is treated in the usual way by 
expanding its time- dependent wave function i/^(r,t) in terms of stationary- 
state Morse oscillator eigenfunctions according to 

ip(r,t) c^(t)e"^“*"^ij)^(r) (4.6) 

n 

where (o = E /h and E is the nth state eigenenergy. Continuum states 
n n n 

are neglected. We showed in chapter 3 that, for the Morse oscillator, 

“n - “e(" + I) - + W ■ 

The dynamical wave function iKr,t) describing the oscillator response 
during a collision is the solution of the time-dependent Schrodinger equation: 


it, 

3t 

where, from chapter 3, 


[- 


h2 32 

2y 3r2 
o 


+ V^(r) + V(t,r)h)(r,t) 


(4.7) 


V^(r) = 

The solutions are invariant with the equilibrium separation r^ and it may be 
set equal to zero. The remaining potential parameters are equated to the 


63 


familiar spectroscopic constants and according to equations (3.4) and 
(3.5). 

The solution of equation (4.7) is reduced in a standard way (see appen- 
dix B) to a set of linear, coupled differential equations for the expansion 
coefficients defined in equation (4.6). Denoting 


''”3 


(4.8) 


and incorporating the form of the interaction potential in equation (4.5a), the 
the coefficients in equation (4.6) vary in time according to 




dt 


E -z(t)/L 

r~ ® 

kk 


> c.(t)e “ 3 V . 

^ 3 nj 


(4.9) 


The probability that an oscillator, initially in state k at t = will 


reside in state n at t = +<», is then P(E) = |c («>) I ^ with the initial 

k->n 

conditions |c.(-«>)p = 6, . , where 6, . is a Kronecker delta. 

' 3 ' kj’ kj 

The matrix elements given in Integral form by equations (4.2b) and (4.8) 
may be evaluated analytically. If, for convenience of notation, we 
define a = y/aL and 3 = x then 


V =3“ I-(-g--- 

nj a 


j 


i - j) V ^ N A+j -n r(l + g + j - £)r(3 - g- l- j- n+jp 

n! ^ Z!(j - Z)!r(l + a + j - n - Jl)r(3 - 2j + £) 


£=o 


where r(g) is the gamma function^^ with argument C, and the normalization 
constants are 

K = fa(3 - 1 - 2 m)- )l /2 
m L r(3 - m) J 

To stay within the maximum exponent constraints imposed by most computers, the 
evaluation of matrix elements with large indices requires the ratios of gamma 
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functions to be reduced to products of algebraic terms and a residual gamma 
function with an argument less than unity. 


4.1.3 Coupling of the Oscillator Motion and the Classical Trajectory 
— z/L 

The term e in equation (4.9) may be evaluated either from equa- 

tion (4.5b) or by the coupled integration of equation (4.3), depending on 
the treatment of R(t). The former case ignores any coupling between the 
oscillator and the classical path, thus assuming R(t) = 1 for all t. The 
latter case requires evaluation of R(t) in terms of the expansion coeffi- 
cients — a task easily done by combining equations (4.2) and (4.6) to give 


Ji m 

where k again denotes the initial state. 

Equation (4.10) characterizes the classical nature of the quantified 
oscillator motion. The motion will become oscillatory as soon as a mixed- 
state condition is produced during the collision and will remain so afterward. 
Near closest approach, large transient excursions of R(t) occur, reflecting 
the oscillator compression and recoil. 


4.1.4 First-Order Perturbation Solutions 

From a practical viewpoint, the convenience of an analytical solution 
warrants even the coarsest assumptions, provided the limits of applicability 
are understood. This study attempts to confirm those limits for a first-order 
perturbation analysis applied to enharmonic oscillators in initially excited 
states. We shall see that the perturbation solutions are quite successful 
within their intended limits and will serve as useful approximation in many 
cases. 
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An analytical solution of equation (4.9) may be obtained If the motion of 
the classical particle Is described by equation (4.4). For an Initial state 
k, the perturbation method further requires that |cj^(t)p ^ 1 and 
|c„(t)|2 « 1 throughout the duration of the collision. Then only the 
Initial and final states are coupled, allowing equation (4.9) to be written In 
the Integral form: 


|Cn(» 




- V 
E nk 


h V. 


kk 


f sech2(|0exp ^ J r^j^(t' )dt’ dt 

— “ L o J 


(4.11) 


with 


V, 


kk 


.ech^(f) 


Equation (4.11) may be Integrated to glve^^ 


P(E) = 
k^ 




(4.12) 


where 


L(o) - 0), ) V - V, , 

n k ^ _ yLu nn kk 

and M(1 +ig,2,i2X) is the confluent hypergeometric series with complex argu- 
ments. The necessity of complex algebra may be avoided when computing the 
modulus |m( 1 +ig,2,^2X)| by noting its relation to the Coulomb wave func- 
tion with zero index. The result is 


|m( 1 + ig,2,i2X)| = 0^(-g,X) 


(4.13) 


where 


‘o <-*•»> -£ ^ 

£=1 
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with Ai = 1 and A 2 = -g. The remaining coefficients are obtained from 

£(£ - l)Aj^ = -2gA^^j - A ^__2 

4.1.5 Numerical Solution Methods 

Solutions to the coupled set of equations (4.3) and (4.9) were obtained 
by first separating equation (4.9) into a separate set for each complex com- 
ponent and adopting the equivalent of a multiple-state, close-coupling 
approach. Numerical integration was accomplished with a polynomial extrap- 
olation algorithm originally developed by Bulirsch and Stoer^^ and given in 
FORTRAN by Gear.^^ Fifth-order polynomials and a required accuracy of one 
part in 10® seemed to optimize the calculation of a selected test case and 
allowed a complete encounter to be computed in 50 to 1000 steps, depending on 
the collision energy and the number of coupled states. Solutions were started 
with the molecule in a pure eigenstate and with the incident particle at a 
distance such that the interaction potential had a value 10"^ times the esti- 
mated value at closest approach. The calculation was terminated at an equal 
distance after the encounter. All values of |c^(t)p were sufficiently 

constant at termination. The closure relation ~ ^ used 

n 

throughout the encounter to monitor accuracy. If the initial state quantiam 
number were k, states from n = 0 to 2k were included at the maximum ener- 
gies considered to ensure that the solution was unaffected by neglected states. 
Of course, in the energy range where the perturbation theory was successful, 
as few as two states were adequate. The computing times required to obtain 
all the matrix elements and to integrate the dynamics of a twelve-state model 
was approximately 0.1 sec/step on a single precision (14-digit) CDC-7600 
computer. 
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4.2 A Comparison With Fully Quantum Mechanical Solutions 

The availability of tabulated results for exact quantum mechanical calcu- 
lations^® over a broad range of collision parameters provides an excellent 
opportunity to evaluate the semiclassical approximation in this application. 
The extent of the examples covered is characterized by the range of the mass 
parameter m = m^m^/m^(m^ + m^ + cases chosen, m varies from 

0.006 for Br 2 “H collisions to 3.7 for HBr-He collisions. (Reference 68 labels 
one data set as Br 2 ~H 2 , but uses a mass parameter corresponding to Br 2 “H.) A 
full range of oscillator frequency and anharmonicity is also represented. 
Figures 4.2(a) to (f) compare the predictions of the semiclassical theory and 
its first-order approximation to a sampling of the results in reference 68 for 
the homonuclear oscillator cases. The semiclassical transition probabilities 
are plotted as functions of the normalized initial kinetic energy of the 
incident particle, E/ho)^. The probabilities from reference 68, hereafter 
referred to as "exact," are shown at energies displaced according to a trajec- 
tory symmetrization scheme (to be discussed). In the paragraphs to follow, 
the comparisons in figures 4.2(a) to (f) are used to evaluate the validity of 
several methods of compensating for the lack of energy conservation in the 
semiclassical approximation and to demonstrate the influence of coupling 
between the recoiling quantum-oscillator and the classical incident -particle 
motion. 


4.2.1 Energy Conservation and the Classical Parameters 

The absence of energy conservation in the semiclassical approximation 
requires an interpretation of the initial relative kinetic energy E assigned 
to the classical trajectory. It may be considered an effective value. 
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(a) H 2 (k = 0) - He collisions; o: k n = 0 1; (b) H 2 (k = 2) - He collisions; o: k n = 2 3; 

O:k^n = 0->-3. ^:k->“n=2->-5. 

Figure 4.2.- Comparisons of semiclassical and quantum-mechanical^® transition probabilities of homo- 
nuclear molecules. All calculations are done for L = 0.02 nm. Open symbols denote points tabu- 
lated in reference 68 and plotted according to equation (4.15). The curves are semiclas- 

sical multistate solutions using classical trajectories coupled to the oscillator motion via 

equation (4.3). The curves are the same without coupling via equation (4.4). The curves 

curves are first-order perturbation solutions given by equation (4.12). 






(c) H 2 (k = 0) “ H collisions; coupled and uncou- 
pled solutions are superimposed, 
o:k“>-n = 0->l; 0:k->n = 0^3, 


(d) Br 2 (k = 0) - H collisions; coupled and uncoupled 
solutions are superimposed, o: k->-n = 0->-l; 

□ :k->^n = 0-^2; Ork-^n^O-^^S. 


Figure 4.2.- Continued. 



(e) N 2 (k = 0) - (N 2 ) collisions; o: k ->• n = 0 -► 1; (f) N 2 (k = 2 or 5) - (N 2 ) collisions 
°:k-»-n=0-»-2. (N 2 ) denotes N 2 treated o:k->n = 2-»-3; ci:k->n=2 

as a structureless particle. A; k n = 5 -► 6. 


Figure 4.2.- Concluded. 



averaged over the trajectory from the true initial value Ej^ to the final 
value E^, when, the molecule undergoes a transition from state k to n. If 
total energy is conserved, E^, E^, and E^ are related by 

\ " ^n 

No formal guidelines are available, however, for simply relating E to the 

exact energies and E^. Perhaps the closest one can come is with the 

method described in reference 42, where the formulation of a linearized 

quantum-mechanical approximation is compared to its semiclassical counterpart. 

Expressions for the transition probabilities given by both approximations 

become similar if 1 is defined by the average velocity u = (u^ + u^)/2. 

Another approach is taken by Mies,^^ who uses the intuitively appealing 

arithmetic energy average E = (E + E, )/2. Combined with equation (4.14), 

n ic 

the total energy can then be related to the average energy according to 

= E + h((o^ + 0)^)/2 (4.15) 

where h (to + (o,)/2 is the oscillator energy averaged over the transition, 
n K. 

Occasionally, even the geometric average E = (E,E has been suggested. 

X n 

Equation (4.15) was chosen here for the comparisons in figures 4.2, where it 
is shown to be generally successful. It correlates the predictions of both 
theories for all initial states, transitions, and mass ratios tested and 
appears applicable for all energies E from threshold up to at least the 
first probability maximum. Figure 4.2(c) shows correlation beyond the first 
maximum. Note however, that when the effect of the oscillator motion coupled 
to the classical trajectory is distinguishable, the coupling must be included 
to preserve the accuracy of equation (4.15) (e.g., see fig. 4.2(a)). The other 
averaging methods are no less accurate, however. Table 4.1 reveals that all 
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the averaging methods described give essentially the same results and that, 
within the range of these comparisons, the best choice cannot be selected. 
Differences in the three methods (or apparently any other method) will only 
become distinguishable at values of where, from thermal considera- 
tions, is beyond the range of practical interest. From a pragmatic view- 

point, equation (4.15) is attractive because, unlike the other averages, it 
provides an energy transformation, E^ - E, Independent of E and allows the 

energies E or E to be written explicitly In terms of E with simple 
iC n 

algebraic form. These features are convenient for additional manipulation 
such as thermal averaging. 


TABLE 4.1.- A COMPARISON OF SYMMETRIZATION METHODS APPLIED TO ANHARMONIC H 2 


Transitions 
k n 

(E - E)/ho) 

1 e 

Observed 
Figs. 2(a)-(b) 
Note (a) 

® ' i + V 

“ ■ I '“n + “k’ 

Note (b) 

E = (E E, )l/2 

n k^ 

Note (b) 

0-1 

1.0 

0.97 

0.99 

1.01 

0-2 

1.3 

1.41 

1.46 

1.50 

0-3 

1.7 

1.83 

1.94 

2.05 

2-3 

2.7 

2.75 

2.76 

2.78 

2-4 

3.1 

3.14 

3.20 

3.26 

2-5 

3.5 

3.51 

3.65 

3.80 

5-6 

5.0 

L ■ ■ 

5.03 

5.06 

5.08 


Note (a) The observed energy difference between the semiclassical and 
exact results for a given probability near threshold. The semiclassical 
results include oscillator feedback. 

Note (b) Computed for E^o)^ =6. 
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4.2.2 Influence of Oscillator Response on the Classical Motion 

The discussion to this point has been confined to homonuclear molecules. 
Figures 4.2 indicate that coupling oscillator motion and the classical trajec- 
tory has a noticeable effect only for the most anharmonic molecule, H 2 , and 
then only when struck by a relatively heavy particle. He. However, semiclas- 
sical calculations for heteronuclear cases are much more sensitive to the 
oscillator response. In the customary semiclassical formulation, the incident 
particle dynamics are related only to its distance from the mass center of the 
molecule, and no regard is given for the location of the impacted nucleus 
(e.g., see eq. (4.4)). In an extreme heteronuclear case where the impacted 
nucleus is extended to a distance similar to the distance of closest 
approach, the incident particle can spatially overlap the impacted nucleus 
i without constraint. Of course, even the approach to this extreme situation 
signals the failure of the assumptions leading to equation (4.4). 

The hydrogen-halides represent examples of diatomic molecules whose 
heteronuclear properties strongly influence the incident particle motion, with 
the effects further augmented by the accompanying large anharmonicity . As an 
example, figure 4.3 illustrates the behavior of HBr-He collisions, where H is 
the impacted nucleus. Similar results were obtained for HCl-He and are assumed 
to be characteristic for all hydrogen-halide-like molecules. Both the semi-’^ 
classical numerical solutions and the analytical theory predict an anomalous 
resonance at low energy when the classical trajectory is obtained from equa- 
tion (4.4). The resonance is a combined result of an improper trajectory and 
the oscillator anharmonicity since similar calculations treating the molecule 
as a harmonic oscillator behaved normally and in accordance with corresponding 
quantum-mechanical solutions.^® Considerable care was exercised in verifying 
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Figure 4.3.- Comparison of semiclassical and quantum-mechanical transition 
probabilities for HBr(k = 0) - He collisions. The impact is between H 
and He. The notation is the same as in figure 4.2; o denotes 
k n = 0 1. 

the resonance as a real solution of the theoretical model rather than a numer- 
ical artifact. The similar behavior of the analytical solution supports the 
conclusion that the effect is a real consequence of the model used. When 
oscillator motion is included via equations (4.3) and (4.10), the resonance 
disappears and the solution is more in accordance with the quantum-mechanical 
results for single-quantum transitions. However, multiple-quantum probabil- 
ities such as Pq 2 still display a low-energy anomalous resonance. The inter- 
polation of E for single-quantum transitions is also shown to be less 
accurate, but equation (4.15) still performs well near threshold. The results 
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suggest that the Interference between the oscillator and the incident particle 
is not fully accounted for, but if it were, equation (4.15) would apply. 

The effects of the oscillator motion are not generated simply by large 
excursions of R(t) during the collision. Figure 4.4 compares the time- 
dependent variation of R(t) with and without the effects coupled to the 
collision dynamics for two extreme cases: (a) H 2 “H collisions, where the 
excursions of R(t) are largest, but the effect is negligible; and (b) HBr-He 
collisions, where the excursions are smaller, but a phase shift is introduced 
in the oscillator motion that severely alters the remaining oscillator 
response. 



Qt/2L 

Figure 4.4.- Transient oscillator effects on the interaction potential; R is 

defined by equation (4.3b). Curves denote the potential term with 

the oscillator motion coupled to the classical trajectory; curves 

denote the potential term without coupling. 
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4,2.3 Applicability of First-Order Perturbation Theory 

Figures 4.2 and 4.3 amply demonstrate the fact that, when the oscillator 
is anharmonic, not only must be small for the first-order perturbation 

theory to apply but the transient motion of the oscillator must also have no 
significant effect on the collision dynamics. When the effects of the oscil- 
lator motion were negligible, the perturbation approximation fails when 

Further increases in E cause to exceed unity before teaching 

a maximum — a familiar feature of equivalent harmonic oscillator models. In 
the application of a perturbation approximation to anharmonic oscillators, a 
probability maximum may appear for P ^ 1, but the maximum is always an 
artifact of the semiclassical perturbation approximation and signals the 
failure of the theory due to the neglect of R(t) variations. These conclu- 
sions are not surprising, but they further constrain the anharmonic oscillator 
perturbation theory to heavy homonuclear molecules such as N 2 , O 2 , and the 
halogens. First-order perturbation calculations for slightly heteronuclear 
molecules such as CO also require careful attention. Just as for H 2 , colli- 
sions of CO with lighter particles (e.g., CO-He collisions) were unaffected by 
the oscillator motion. However, figure 4.5 shows that perturbation calcula- 
tions of collisions with heavier particles (e.g., CO-Ar collisions) display 
large errors due to the coupled oscillator motion and also due to an increased 
coupling of the nonadjacent oscillator states that are not included in the 
perturbation approximation. As figure 4.5 shows, the two effects compensate 
each other for the initial state chosen. In the CO-Ar case, the perturbation- 
theory errors are not accompanied by anomalous probability maximums in the 
energy range of practical interest. 
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Figure 4.5.- Semiclassical transition probabilities for CO(k = 1) - Ar colli- 
sions; L = 0.02 nm. 

4 . 3 Suinmary 

The semiclassical approximation has been applied to vibrational transi- 
tions induced in anharmonic oscillators by collinear collision with inert 
atoms. Multistate numerical solutions have been compared with exact quantum- 
mechanical calculations of an equivalent collision model for a wide range of 
initial molecular states and collision partners. The comparisons allow a 
comprehensive assessment of the semiclassical approximation for the anharmonic 
oscillator model. The semiclassical predictions accurately reproduce the 
quantum-mechanical transition probabilities for all initial collision energies 
from the threshold to at least the first probability maximum if either the 
semiclassical collision velocity or energy is interpreted as a simple average 
of the exact initial and final values. The accuracy of the correlation 
between theories is not sensitive to the choice of averaging method. 
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The semiclassical approximation, in its usual form where the classical 
trajectory is computed independently, was found to be applicable to heavy 
homonuclear molecules such as N2, O2, and the halogens on impact with lighter 
partners. Lighter homonuclear molecules such as H2 showed poorer agreement 
when impacted by a heavier collision partner. Heteronuclear anharmonic mole- 
cules such as the hydrogen-halides displayed anomalous resonances at low 
energy that do not appear in their harmonic counterparts. The accuracy of the 
semiclassical approximation for light or heteronuclear anharaonic molecules 
was significantly improved by coupling the effects of the time-dependent 
average motion of the recoiling oscillator to the classical trajectory. 

A convenient, analytical, first-order, perturbation analysis for anhar- 
monic oscillators was found to be accurate for small-transition probabilities, 
but only if the effects of the oscillator motion on the classical trajectory 
were unimportant. The analytical approximation is therefore not applicable 
to significantly anharmonic and heteronuclear molecules and must be applied 
with care for slightly heteronuclear molecules such as CO. 
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CHAPTER 5 


VIBRATIONAL QUANTUM NUMBER DEPENDENCE OF ENERGY-TRANSFER RATES 

In chapter 1, we discussed the sparseness of experimental and theoretical 
efforts to describe the rates of vibrational energy transfer from Initially 
excited vibrational states. The objectives of this chapter are therefore 
to examine the factors Influencing the prediction of such rates and to 
evaluate the validity of several analytic formulas In popular use for 
estimating their quantum number dependence. However, the approach to be 
taken must be limited by pragmatic considerations. For example, accurate 
V-T rate-coefficient calculations by any theoretical model are clouded by 
uncertainties In the shape and magnitude of the Interaction forces between 
colliding pairs for all but a few simple cases. Thus, we can examine only 
the qualitative features that are not masked by Interaction potential 
uncertainties. Furthermore, even an extended collision model must retain 
some approximations, particularly regarding the collision geometry. If It 
Is to remain computationally practical In the prediction of thermally 
averaged rate coefficients. Hence, attention Is confined here to the col- 
llnear semlclasslcal treatment described In chapter 4 that accurately repro- 
duces all the main characteristics of vibrational energy transfer to Ini- 
tially excited oscillators, but offers the advantage of being further reduc- 
ible to yield closed-form analytic solutions. The analytic solutions are of 
particular Interest because of their practical Importance In the numerically 
cumbersome analysis of macroscopic nonequilibrium processes for which rate 
Information for several modes of energy transfer must be economically 
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provided. The complete semlclasslcal model, requiring numerical solu- 
tion, is applied both to an examination of the qualitative nature of upper- 
state transitions and as a basis for evaluating the accuracy of the analyt- 
ical solutions. 

In the sections to follow, the features of the collision model that 
appear most important to the dynamics of a vlbratlonally excited oscillator 
are first discussed, followed by a description of several approximations, 
each of which retains one or more of the features considered. Approximate 
values of the interaction potential parameters and their range of uncertainty 
are then estimated by comparing the predicted ground-state rate coefficient 
with a comprehensive set of experimental values. Collisions of CO with He 
are chosen as the example because of the abundant data available. The 
implied potential parameters are then used to compare the numerical model 
with some experimental excited-state rate coefficients and with the analytic 
predictions. Finally, the effects of multiple-quantum transitions from 
excited states on a vibrational relaxation process are considered both for 
molecules like CO, where the effect is secondary, and for molecules like 
the halogens, where the effect can be dominant. 

5.1 Collision Model 

5.1.1 Features Influencing the Excited-State Collision Dynamics 

As the quantum number of the Initial oscillator state is Increased, 
several aspects influencing the oscillator dynamics and Its interactions 
with the Incident particle become Increasingly Important. For example, 
the wave functions that describe vlbratlonally excited eigenstates become 
more extended in the oscillator coordinate. Consequently, when the 
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oscillator is distorted by a collision, the wave-function overlap is greater 
not only with adjacent eigenstates but with more remote states as well. This 
feature is reflected by the increased magnitude of the matrix elements 
dynamically coupling the eigenstates which, in turn, accounts for the greater 
probability of V-T energy transfer through both single- and multiple-quantum 
transitions. Furthermore, the increased coupling of nonadjaaent states 
during the collision can affect the final occupation of states adjacent to 
the initial state and thereby influence the rate of single-quantum transi- 
tions. Thus, a calculation of the oscillator dynamics from an excited 
initial state must include multiple-state interactions at collision energies 
where they are normally unimportant for oscillators in the ground state. 

The degree with which multiple-quantum transitions influence the 
oscillator dynamics during a collision depend, in part, on the form of the 
interaction potential. A common practice, often used to simplify the 
analysis of ground-state oscillators, is to consider the oscillator motion 
small compared to the range of interaction and linearize the interaction 
potential in the oscillator coordinate. In a harmonic oscillator, this 
treatment has the effect of equalizing all the diagonal matrix elements and 
forbidding multiple-quantum transitions. The occupation of nonadjacent 
oscillator states is then possible only through a sequence of single-quantum 
steps during the collision. Nonlinear interaction terms remove these 
restrictions and modify final-state occupations in two related ways. First, 
all the nonadjacent states are directly coupled, thereby increasing their 
accessibility. Second, the diagonal matrix elements are no longer equal, 
leading to additional phase distortions in the quantum-mechanical oscillator 
motion which modify the probability of transition. The additional phase 
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shifts depend on the product of the difference between diagonal matrix 
elements and the strength of the interaction. They appear explicitly in a 
semiclassical impact parameter treatment described by Bates and applied to 
enharmonic oscillators by Mies.^^ The formulation has been reviewed in 
section B.3 of Appendix B. 

All the foregoing effects are amplified when oscillator anharmonicity 
is included. Nonadjacent states become coupled even for linearized inter- 
actions and the larger difference between the diagonal matrix elements 
creates phase distortions that can become a significant fraction of the 
unperturbed oscillator period. Mies^^>^^ has shown the influence on transi- 
tion probability predictions to be large even for oscillators initially in 
the ground state. A second, and in some cases greater, effect of anharmon- 
icity is its influence on the variation of eigenenergies with quantum number. 
Since transition probabilities and the related rate coefficients are known 
to depend on the amount of energy transferred, a lowest-order effect of 
oscillator anharmonicity may be demonstrated by simply inserting anharmonic 
oscillator eigenenergies into a harmonic oscillator theory such as that given 
by Schwartz et at. The results deviate substantially from the simple 
Landau-Teller relation for the rate coefficients given by equation (1.2) and 
repeated here as 

k (T) = m k (T) (5.1) 

m,m-i i>u 

where k (T) denotes the rate coefficient for transitions from state 
m,m-l 

m to m-1 and is a function of the kinetic temperature T. However, the 
simple ad hoc insertion of anharmonic eigenenergies into a harmonic oscilla- 
tor model is not always a sufficient means of accounting for anharmonicity. 
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The Influence of anharmonicity on the interaction matrix elements, which in 
turn effects both the magnitude and phase of the oscillator motion, is often 
so great that an anharmonic oscillator model must be used from the start. 
Fortunately, oscillator anharmonicity and nonlinear interaction potentials 
present only a slight increase in computational difficulty, particularly if 
a Morse oscillator and an exponential form of the interaction are adopted. 

The necessary matrix elements are then conveniently expressed in closed 
algebraic form just as for harmonic oscillators. 

Finally, an oscillator potential creating anharmonicity also admits to 
the existence of continuum states. We shall neglect their contribution to 
the energy transfer process, however, since they are energetically inacces- 
sible by a large margin for the combinations of collision energies and 
initial states considered here. Although their effects have not been evalu- 
ated, the occupation of continuum states is presumed to be as small as the 
nearby bound states, and no bound states near the continuum were found to 
influence the d 3 mamics of any states at the quantum levels of interest. 

5.1.2 Aspects of the Semiclassical Numerical Model 

To obtain V-T rate coefficients, we calculate the associated transition 
probabilities using the semiclassical collision model described in chapter 4, 
wherein the trajectory is constrained to collinear encounters. One of the 
penalties of using a semiclassical approximation was shown to be that total 
energy is not conserved; but as shown in chapter 4, the effects of that 
omission are easily and accurately compensated for by interpreting the 
relative collision energy as an average of the known initial and final 
values. A far more severe limitation of the semiclassical theory was 
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found to be its incomplete treatment of the interaction when the oscillator 
is very heteronuclear (e.g., hydrogen-halides). Such cases are avoided 
here and have presented numerical difficulty in exact treatments 

The implications introduced by a restriction to collinear encounters are 
not as well understood, but the restriction is necessary if the quantum^ 
number dependence of thermally averaged rate coefficients is ever to be 
obtained in a reasonable computing time. Clearly, a more realistic approach 
would Include a three-dimensional collision geometry in which simultaneous 
rotational transitions are coupled with the vibrational motion, but the large 
number of rotational states that become accessible at collision energies 
sufficient to cause vibrational transitions would make our objective of 
studying thermally averaged rate coefficients Impractical for all but a few 
special molecules, like On the other hand, so long as the rotational 

eigenenergies of the undisturbed molecule are well described by a rigid-rotor 
model (suggesting that the rotational and vibrational motions are separ- 
able) , and the molecular properties are such that vibrational energy is 
traded mainly with translation, the disparity between a collinear and a 
three-dimensional theory is not expected to be very sensitive to the initial 
vibrational quantum number. Chapter 6 provides greater insight into the 
necessary conditions. For our purposes, the collinear predictions are 
normalized according to the ratio k^ ^_^/mk^ thus avoiding the predic- 
tion of absolute rate coefficients and, hopefully, much of the absolute 
error associated with the collinear restriction. Such a ratio also absorbs 
the lowest-order quantum-number dependence suggested by equation (5.1). 
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5,1.3 Thermally Averaged Rate Coefficients from a Collinear Semiclassical 
Model 

With the possible exception of molecular beam analyses, the applications 
of an inelastic collision model usually require results in the form of a 
thermally averaged rate coefficient. A general formulation of the averaging 
integral is well known, but here the restriction to collinear trajectories 
and the use of a semiclassical approximation require some special considera- 
tion. Generally, the rate coefficient for a kinetic temperature T may be 

written in terms of the energy parameter e = E /kT and an energy-dependent 

m m 

cross section a(Em), where n denotes the final quantum state and E 
nrti 

is the relative kinetic energy before a collision with an undisturbed 
oscillator in a pure eigenstate m. The rate coefficient is then similar to 
equation (2.2), that is, 

oo 

k (T) = C I o(E^)e e ® de (5.2) 

m,n J mm ^ 

0 

where the average thermal speed is C = (SkT/iry)^/^ and y is the reduced 
collision mass defined by equation (2.3). A further requirement for the 
collision model is that it conform to the detailed balance relations. 
Originating with the reciprocity theorem, the requirements of detailed 
balance propagate through three levels of microscopic detail, giving the 
general physical relations for spinless nondegenerate collision partners as 


P(Em) = P(En) (5.3a) 

nt>n n“>m 


E a 
m 


(Em) = E a(E„) 
umi n->ni 


k 

m,n 


-k(d /kT 
(T) e “ 


= k 


n,m 


-ftw /kT 
(T) e “ 


(5.3b) 


(5.3c) 
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where P 


Is the transition probability from state m to n and hu Is 

m 

the oscillator energy of state m. 

The colllnear collision geometry produces semiclassical transition 
probabilities that behave according to equation (5.3a), but the restriction 
to a zero impact parameter leaves the cross section required by equation (5.2) 
undefined. One common solution is to adopt an effective hard-sphere cross 
section and compute the inelastic cross section according to 


c^(Em) = ^ 


m->n 


m^ 
umi 


and 




n-»in 


n-^ 


Equation (5.3b) requires that 


(5.4a) 

(5.4b) 


I 

a 

o 


[1 + 


h (oj - 
n 


oj )/E ]a 
m no 


thus suggesting that the ”hard-sphere" size must depend on the collision 
energy and transition in question! This contradiction results from the 
colllnear approximation, but the error is negligible when |h (co^ - o)^ ^) |/E^ « 1. 
When the ratio approaches unity, the transition probability Is typically so 
small that the Integral In equation (5.2) Is unaffected. 

Equation (5.2) must be modified further to compensate for the lack of 
energy conservation Inherent In the semiclassical approximation. This 
discrepancy is easily and accurately corrected by interpreting the semi- 
classical relative collision energy E or speed u as an average of the 
Initial and final values. The results in chapter 4 demonstrate that, while 
the correction can be large, the method of averaging has no apparent effect 
on the outcome for vibrationally Inelastic collisions at all*' energies from 
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threshold up to the limits of practical interest. For convenience, we use 
an arithmetic energy average. Denoting the total energy as E^, the 
semiclassical approximation is brought into close agreement with an equiva- 
lent quantum-mechanical calculation by the interpretation 


E = - h(o) + 0) )/2 (5.5) 

T m n 

Combining equations (5.2), (5.4), and (5.5) then gives the thermal averaging 
prescription for a colllnear semiclassical collision model: 


‘' 41 ^ 


= a C e 


ho) /2kT 
mn 


/■’<!> f l^l) 


-e 


de 


(5.6) 


where e = E/kT and co = o) - co . To make the satisfaction of equa- 

mn m n 

tion (5.3c) by (5.6) more obvious, the lower integration limit in equa- 
tion (5.6) has been set to zero even though the independent variable trans- 
formation from e to e via equation (5.5) produces a limit of 
m 

±|ha)^/2kx|, depending on the sign of o)^^. The negative limit may clearly 
be reset to zero, but even when the limit is positive, the probability 
threshold is nearly twice the limit, so that again setting it to zero has no 
effect on the integral. 


5.2 Analytic Approximations 

Of the many analytic approaches appearing in the literature (see ref. 45 
for a partial summary) , three that stand out in their application for 
estimating the V-T rate coefficient variations with quantum number are 

(a) the semiempirical formulas for Morse oscillators of Keck and Carrier, 

(b) the perturbation treatment of Morse oscillators developed by Mies,^^ and 

(c) the exact sol 4 tion to a linearly forced harmonic oscillator obtained by 
Kemer.^^ Each approach retains one or more of the aspects of special 
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Interest to this application. They share the common feature that all 
Incorporate colllnear collision geometry and all are based on an exponentially 
repulsive Interaction potential (later referred to as potential I) of the form 

Vj(x) = A (5.7) 

where the coordinates are defined In figure 4.1. 

5.2.1 Keck-Carrler Formula for Anharmonic Oscillators 

The formula obtained by Keck and Carrier®** is an adaptation of the 
distorted-wave harmonic oscillator theory of Schwartz et at. for a Morse 
oscillator. It Includes an empirical fit to the numerical solution of an 
Integral equation for the "adiabaticity factor" and provides a particularly 
simple formula for estimating single-quantum transition rates from an 
arbitrary initial state. Keck and Carrier made no claim for the suitability 
of their formula in applications beyond a demonstration of the role of 
vibrational nonequilibrium in a dissociating gas; but the formula was subse- 
quently applied by Bray®® in a pioneering and detailed calculation of a 
vibrational relaxation process for anharmonic oscillators, apparently because 
of its simplicity and for lack of a better estimate. For similar reasons, 
the Keck-Carrler formula has since gained widespread use in the detailed 
analysis of upper-state kinetics in lasers. Its consideration here is 
motivated primarily by the number of kinetic models that incorporate it. The 
Keck-Carrier formula can be written in a form similar to equation (5.1) as®® 

/l - X \ F 

k fT) =m(.j ~)^K „(T) (5.8) 

m,m-l yi - mx^yF^ 1,0 
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where F Is obtained from the empirical formula 
in 


. i (3 - 


47rr|/3 


(5.9) 


in which 


n = -w „ ,L(y/2kT)l/2 


(5.10) 


The transition frequency w _ = oj - 03 , is computed for a Morse 

mjm—i m m — 1 

oscillator as done previously by the expression 


u) = 0 ) /m + - 0 ) X (m + •“) 

m e\ 2/ e e\ 2/ 


(5.11) 


5.2.2 Mies Perturbation Solution for Anharmonic Oscillators 

The closest approximation to the numerical model used here is a 
semiclassical first-order perturbation treatment developed by Mies.^^ It 
properly includes the effects of anharmonicity but, by the nature of first- 
order methods, it neglects the influence of states other than the designated 
initial and final states. Furthermore, to obtain an analytical solution, the 
classical path must be computed independently from the motion of the oscilla- 
tor. The theory is therefore applicable only to single-quantum transitions 
in which the transition probabilities are small compared to unity. As 
demonstrated in chapter 4, the independent classical path further restricts 
its application to nearly homonuclear oscillators such as CO (and, of course, 
all homonuclear molecules) colliding with atomic particles of lighter mass 
than either of the molecular nuclei. The appearance of a probability 
maximum signals the failure of the theory. Despite these shortcomings, we 
shall see that the Mies solution still provides a more useful approximation 
of the numerical predictions than the other analytic formulas investigated. 
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A convenient form of the Mies result was given by equation (4, 12) and is 
rewritten here for the transition probability from state m to n, where 
n - m ± 1 as 

2 


P(E) 

nr>n 


V - “I 

mn 27TgyLu ^ 

“ V~h sl^CTTg) 

_ mm 


(5.12) 


As with the numerical model, equation (5.12) produces energy-dependent 

transition probabilities, while a temperature-dependent rate coefficient is 

desired. No analytic solution of the integral equation (5.6) with P(E) 

nH-n 

given by equation (5.12) is apparent, but a reasonably accurate technique 

(labeled the "method of steepest descent") for obtaining an analytic 

approximation^^ is based on the well-defined maximum contained in the 

integrand of equation (5.6). The value of e at which the maximum occurs is 

determined primarily by the exponential arguments. The remaining function 

is slowly varying over the range of the integrand and may be evaluated at 

the single value locating the peak. The exponential argument is then 

expanded to second order about the peak and the term integrated analytically. 

In this application, the notation is simplified with the substitutions 

e = ho) /2kT and n = L(y/2kT)^/^. The exponential nature of 

mn mn mn 

equation (5.12) is also simplified by noting that, in the energy range where 

the perturbation analysis is applicable, the transition period t^ = 

is typically less than the effective collision period t = 2L/n. Thus, 

c 

TTg = t /t >1 and sinh(TTg) (l/2)e’'^§. Equations (5.6) and (5.12) are then 
c p 


combined to give 


k(T) = a C 
m,n o 


V 2 
mn , n 

— — 4tt 

V e 
mm mn 




— — — 1 / 9 

2 -e +2Trrie ' ^ - 


de (5.13) 
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The integrand peak is located at 


e = [p(irto L)2/2kT]^/3 
p mn 


(5.14) 


Using the procedure described, the approximate solution to equation (5.13) 
becomes 2 

\ mm mn/ ^ ^ tr tr 


-3e +e 
p mn 


|l + erf [(3Ep/4)l/^] | 


(5.15) 


where g = e /tt and X = e ^(V - V )/(Tre V )- The error function in 

®p p p p mm nn ' mn nm 

equation (5,15) is close to unity for most cases. Equation (5.15) has 
log k a as expected and satisfies equation (5.3c). The temperature 

at which a given collision speed is coincident with the peak of the integrand 
in equation (5.13) defines the most effective speed at that temperature; this 
temperature will also be useful and can be identified from equation (5.14) 
as 

T = Mu^/(27rk(a3 (l) (5.16) 

p mn 

Comparisons of the approximate integration in equation (5.13) with exact 
numerical integrations show that the approximate method is most accurate at 
low temperatures. The first-order perturbation formula, equation (5.12), is 
most accurate at low energies, thus further contributing to the accuracy of 
equation (5.15) at low temperatures. 


5.2.3 Kerner Solution for Linearly-Forced Harmonic Oscillators 


The final analytic formula to be considered is an exact solution 
obtained by Kerner^^ for a harmonic oscillator that undergoes a forcing 
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function linear In the oscillator coordinate. That condition may be satisfied 
In situations where r/L « 1 In equation (5.7). The potential may then be 
linearized according to 

V(x,r) = A e'^^^^Cl - yr/L) (5.17) 


Kerner's solution was applied by Treanor^^ in a semiclassical collinear 
approximation using equation (5.17). Within the framework of the collision 
model, the resulting formula exactly calculates the probability of transitions 
between arbitrary states with the interaction of all states included. Thus, 
it can be applied at high collision energies where the interactions of more 
than two states influence the oscillator dynamics. Despite the approximate 
nature of the harmonic oscillator model, wherein direct multiple-quantum 
transitions and the unbalanced coupling of higher and lower states caused by 
anharmonicity are excluded, the Kerner solution remains useful because it 
offers the only analytic means for estimating transition probabilities at 
high energies. Examples will be shown where multiple-quantum transitions 
and oscillator anharmonicity are not dominant, allowing accurate prediction 
by the Kerner solution. 

Kerner and Treanor write the probability for transitions between two 
arbitrary states m and n as 



m!n! 



)^(n - j)!j!(n - j)! 



(5.18) 


where J is the lesser of the quantum numbers m and n. The parameter 
is the energy absorbed by a classical harmonic oscillator divided by one 
quantum of vibrational energy. For a collinear collision and the interaction 
of equation (5.17), Rapp^® obtains 
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= 2(2ira)Lyii)^ e (5.19) 

where u is the reduced mass of the oscillator and o) is the oscillator 
o 

frequency. The accuracy of the model, when applied to highly excited 

oscillators, is substantially improved if the effective oscillator frequency 

is corrected for anharmonicity for each initial state m according to 

u) = 0 ) (1 - 2x m) , Without the correction, the excited-state rate coefficients 
e e 

would simply behave according to the Landau-Teller relation, equation (5.1), 
at low temperatures where the effective values of are all less than 

unity and give k /mk, _ < 1 for large E . An inconvenience of the 
Kerner formula is its incompatibility with the approximate integration method 
of equation (5.6) for obtaining a rate coefficient. A simplified version of 
equation (5.18), assuming E^ << 1, permits an approximate analytical 
solution. However, the calculations are then restricted to a thermal range 
where multiple-quantum effects are insignificant and the theory loses its 
advantages over perturbation solutions. In the comparisons to follow, we 
have therefore resorted to a numerical integration of equation (5.6) when 
the Kemer solution is applied. 

5.3 Comparisons with CO-He Experiments 

In this section, the ability of the theoretical model to reproduce 
experimental rate coefficients is tested. Unlike past comparisons of 
vibrational rate coefficients with theory, we now have access to at least 
one set of experimental values for excited initial states.^ To test the 
consistency of the theory and experiment for all vibrational states, however, 
the effective interaction range L and the hard-sphere cross section 
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are determined from the abundant collection of measurements dominated by 
transitions between the ground state and first vibrational state. The 
Interaction parameters required to match the ground-state experiments are 
then applied In comparisons with the excited-state rate measurements. 

5.3.1 Effective Interaction Potential Parameters 

The computational convenience gained from the simplified Interaction 
potential 1, equation (3.7), justifies Its use, but as a consequence of Its 
simplified form, the predicted rate coefficients cannot be expected to 
reproduce the experiments at all kinetic temperatures. Transitions Induced 
In an oscillator depend to a large extent on the potential gradient near 
the distance of closest approach; while in a collinear collision, the distance 
of closest approach is determined solely by the coordinate where the poten- 
tial magnitude equals the Initial kinetic energy of the collision. Thus 
both potential features are Important. However, the magnitude of a purely 
repulsive potential, such as equation (5.7), and that of a more realistic 
potential with an attractive well may be the same at the closest approach 
distance, but have a significantly different gradient. Consequently, where 
collisions are averaged over a range of energies, the predicted variation 
of rate coefficients with kinetic temperature will be different for the 
two potentials. By matching theory and experiment in several thermal 
ranges, and by using more than one potential form, an indication of the 
degree of uncertainty in rate coefficients attributable to potential errors 
can be obtained. For that purpose, we consider a second potential given by 

(x -x) /L (x -x) / 2L 

Vj.j(x) = De -2De (5.20) 
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Potential II is a Morse-type interaction with an attractive well of depth 
-D at coordinate As with equation (5.7), the exponential form allows 

matrix elements to be calculated analytically. 

Predictions by the numerical anharmonic oscillator model with the 
oscillator initially in the first eigenstate m * 1 are compared with 
experiment in figure 5.1. When potential I, equation (5.7), is used, the 
rate coefficients are independent of the magnitude A, so that only the range 
L requires specification. Similarly, the predictions using potential II 
are independent of x^, but require both L and D to be specified. The value 
D/k = 100 K is representative of well depths inferred from viscosity 
measurements.®^ The two potential gradients are different by about 20 percent 
at closest approach for the typical conditions considered. Figure 5.1 
demonstrates the expected results. No unique set of potential parameters 
reproduces the experiments over the complete thermal range, but the more 
realistic potential II comes the closes. The required values of L fall 
between 0.02 and 0.03 nm, depending on the thermal range considered. 

As an interesting aside, note that the low-temperature departure of 

“ 1 / 3 

the experimental rate from a variation proportional to T ' is also 
followed by the theory using simple repulsive potentials. As Shin®^ points 
out, these low- temperature departures do not necessarily depend on weak 
attractive forces normally omitted from the interaction potential; they even 
occur with a repulsive potential when the thermal averaging integration is 
done accurately for low collision energies. We know, however, that real 
interaction potentials usually contain an attractive component and it will 
augment this low- temperature behavior. 
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Figure 5.1.- A comparison of experimental rate coefficients for C0(M = 1)-He 
transitions to the ground vibrational state with predictions from the 
numerical model of chapter 4. The solid and long-short dashed lines 
were computed using the repulsive interaction potential I, equation 
(5.7). The short-dashed line was computed using the Morse interaction 
potential II, equation (5.20). Hard-sphere collision cross sections 
were chosen for each potential to match the experiment at T = 1000 K. 
Experimental values are from: o reference 16, • reference 17a, 

^ reference 17b, ♦ reference 17c. 


5.3.2 Comparisons with Excited-State Rate Measurements 

Normalized rate coefficients, predicted for initially excited CO at 

T B 300 K, are compared in figure 5.2 with the room temperature measurements 

of Hancock and Smith. The parameter k /mk is much less sensitive 

m,m-i 1,0 

to interaction uncertainties than the absolute rate coefficients and varies 
in a simple, nearly linear manner with initial-state quantum number ra. 

The nearly linear quantum-number dependence, increasing with m at room 
temperature, is predicted for all the interaction potentials examined and is 
believed to be an accurate description of the real behavior. As figure 5.2 
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INITIAL-STATE VIBRATIONAL QUANTUM No„ m 

Figure 5.2.- A comparison of experimental rate coefficients at T = 300 K for 
C0(m)-He transitions from vibrational states m to m - 1 with predictions 
from the numerical model using repulsive potential I, equation (5.7). The 
excited-state data are from reference 18 and have been normalized using 
the experimental value of Milliken^^ ’ ^ ^ (fig- 5.1). 

shows, the experimental excited-state values compare favorably in magnitude 
with the predictions, but their trend is inconsistent with a linear extrapo- 
lation to m = 1. A highly nonlinear extrapolation is contrary to any 
prediction of the collision model at any temperature. Although the 
collision model contains many simplifications awaiting refinement, the 
behavior implied by the experimental rates appears also to require further 
verification and extension. In the interim, the theoretical predictions 
of excited-state rates seem to be qualitatively reasonable and self consistent 
despite their quantitative uncertainty. Unfortunately, their verification 
by experiment remains inconclusive. 


5.4 An Evaluation of the Analytic Approximations 

The computational expense of the numerical model makes it impractical as 
a general means of estimating excited-state rate coefficients. Instead, it is 
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used here as a basis for evaluating the more convenient but less complete 
analytic formulas. The predicted rate coefficient variations with quantum 
number for several models are illustrated in figure 5.3 for two extreme 
temperatures. The differences in the various models depend strongly on the 
kinetic temperature, but they all predict a simple mono tonic change with 
quantum number. The analytic approximations are therefore more clearly 
evaluated by choosing the highest initial quantum number of practical interest 
and then comparing the predictions for a range of temperatures. For CO, 



T=300K 


T=3000K 


INITIAL-STATE VIBRATIONAL QUANTUM NO., m 


Figure 5.3.- The CO(m)-He rate coefficient dependence on quantum number 
predicted by several collision models. The solid lines represent the 
enharmonic numerical model, the long-short dashed lines represent the 
Kerner harmonic oscillator solution, equation (5.18), and the dashed 
lines are from the formula of Keck and Carrier, equation (5.8). The 
potential range L = 0.02 nm was used in all cases. 


Rich et have shown that energy transfer from vibrational levels as 

high as the twentieth can influence the net energy balance in an electrically 
excited CO laser system. Choosing m = 20 as an example, the single-quantum 
rate coefficients predicted by all the collision models are compared 
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was chosen so that 


in figure 5.4. The independent parameter (hco^/kT)^^^ 
predictions by the Keck-Carrier formula, equation (5.8), appear as a nearly 
straight line. A comparison of the rates from the numerical model using 

T.K 



Figure 5.4.- A comparison of excited-state rate coefficients for 

CO(m = 20)-He predicted by several collision models. The potential 
range was L = 0.02 nm in all cases. 

potentials I and II shows the moderate sensitivity of k ,/mk, to the 
^ ^ m,m-l 1,0 

form of the potential for one potential range L at all temperatures. Not 
shown is the great sensitivity of the magnitude of k^ ^_^/mk^ ^ to other 
potential ranges at any temperature. Note, however, that the qualitative 
nature of the predictions are undisturbed by the form of the potential and 
are therefore considered realistic. As expected, the Mies solution, 
equation (5.15), accurately reproduces the numerical results at low tempera- 
tures, but fails at higher temperatures where multiple-state interactions 
begin to affect the single-quantum transitions. The departure is signaled 
when transition probabilities approaching unity influence the thermal 
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averaging Integral, equation (5.6). Since CO is not very anharmonlc , the 
Kemer harmonic oscillator model, equation (5.18), frequency-corrected for 
anharmonlclty at m = 20, works well over the entire thermal range. Note 
that the anharmonlc correction must be included, however, as all predictions 
are significantly above the result stated by equation (5.1) for a single- 
frequency harmonic oscillator. Finally, figure 5.4 shows that the Keck- 
Carrier formula, equation (5.8), is too crude an approximation for large 
initial quantum numbers. 

The degree of oscillator distortion caused by the collision of a light 
helium atom with a CO molecule has an insignificant effect on the classical 
trajectory. This fact is made evident by the small difference at low 
temperatures between the numerical model (where the effect is Included) and 
the Mies solution (where it is neglected). An example in which the coupling 
is larger is illustrated in figure 5.5 for C0(m = 20) -Ar collisions. In this 
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Figure 5.5.- A comparison of excited-state rate coefficients for 
C0(m = 20)-Ar. Potential I, equation (5.7), was used with 
L = 0.02 nm in all cases. 
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situation, none of the analytic models do well at low temperatures because 
the effects of oscillator distortion on the classical path modifies the 
transition probabilities even near threshold. The small corrections are 
then greatly amplified by the thermal averaging Integral at low temperatures. 

The small anharmonlclty of CO(x^ » 0.0062) has Influenced the preceding 
examples mainly by altering the energy spacing between excited eigenstates. 
Anharmonlclty also modifies the absolute magnitude of the rate coefficients, 
but that effect Is not apparent In k ./nk, . An example In which the 

anharmonlclty Is large Is Illustrated In figure 5.6 for H^Cm = 10) -He 
(x^ = 0.0268). In this case, the frequency-corrected harmonic oscillator 
model Is Inaccurate at all temperatures. The large spacing between elgen- 
energles In suppresses the onset of multlstate Interactions 
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Figure 5.6.- A comparison of excited-state rate coefficients for 
HgCm “ 10)-He. Potential I, equation (5.7), was used with 
L B 0.02 nm In all cases. 
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at high temperatures, making the Mies solution an accurate reproduction of 
the numerical results over the entire thermal range. The difference In mass 
between the He and H nuclei produces only moderate coupling between the 
compressed oscillator and the classical path (e.g. , see fig. 4.2). 

As the preceding comparisons Indicate, one cannot generally choose a 
single analytic model for estimating excited-state rate coefficients that 
Is applicable to all collision pairs. The situations where a model should 
not be used are easier to Identify. Clearly, the Keck-Carrler formula, 
equation (5.8), Is too approximate In all the examples. The Kemer harmonic 
oscillator solution, equation (5.18), with anharmonlclty-corrected frequen- 
cies Is reasonably accurate unless the anharmonlclty Is large. The Mies 
anharmonlc oscillator solution, equation (5.15), Is a poor approximation 
when multiple-state Interactions become Important. Finally, no analytic 
model based on the semlclasslcal approximation will be realistic when the 
oscillator dynamics have a significant Influence on the classical path of 
the Incident particle. This restriction limits all the models considered 
to collision pairs In which the mass of the Incident particle Is not 
significantly greater than the mass of the Impacted nucleus and to oscilla- 
tors that are not extremely heteronuclear. 

5.5 Multiple-Quantum Transitions 

In the preceding section, only transitions to an adjacent state were 
examined. Here, we Investigate the relative Importance of multiple-quantum 
transitions, particularly for oscillators in highly excited states. The 
probabilities of multiple-quantum transitions are compared in figure 5.7 
both for CO(m)-He collisions in which the oscillator is initially in an 
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EFFECTIVE TEMPERATURE. Tp. K FROM EQ. (5.16) 
3(X) 1000 30CX) 5000 



Figure 5, 7.- Multiple-quantum transition probabilities for CO(m)-He 

collisions using the enharmonic numerical model with potential I and 
L = 0.02 nm. The effective temperature Tp locates the most effective 
collision speed contributing to the thermally averaged rate coefficient 
at the temperature designated. 


excited state and in states near the ground state. The collision speeds 
contributing most to the thermally averaged rate coefficient at a selected 
temperature are indicated by the effective temperature T^. In the thermal 
range considered, multiple-quantum transitions to the ground state are always 
improbable compared to single-quantum transitions from the first vibrational 
level, but the situation is clearly different when the oscillator is 
initially in the twentieth quantum state. However, thermally averaging the 
transition probabilities in figure 5.7 reduces the apparent importance of 
multiple-quantum transitions in a relaxation process. Figure 5.8 illustrates 
the resulting rate coefficients for two potential ranges, using potential I 
and values of obtained from the experimental match in figure 5.1 at 

T « 1000 K. The amplified uncertainty caused by the interaction potential 
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Figure 5.8.- Multiple-quantum rate coefficients for C0(m)-He. Potential I 
was used in the anharmonic numerical model. The hard-sphere cross-section 
values Oq for each potential range are those required to match the 
experimental rates in figure 5.1 at T = 1000 K. 


and its influence on the Implied value of is most obvious, but the 

qualitative features are again consistent for both potential ranges. For 
oscillators like CO, multiple-quantum transitions provide a significant 
path for energy transfer only at very high temperatures according to these 
predictions. 

A temperature marking the onset of competitive multiple-quantum 
transitions is the characteristic vibrational temperature of the oscillator, 
here defined as 0 = hu /k (for CO, 9 = 3122 K) . An oscillator in which 

V e V 

multiple-quantum transitions will dominate the relaxation process can then be 

identified if 0^ is small compared to the thermal range of interest. One 

extreme example is Brg for which 0^ = 465 K. Since the anharmonicity is 

also small in Br„(x = 0,0033), the Kerner harmonic oscillator model has 
2 e 
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been used to obtain the Br^-He rate coefficients displayed in figure 5,9. 
Two- and three-quantum transitions from the tenth vibrational level are 
shown to be significant even at room temperature, and the temperature 
dependence of the single-quantum rate (m n = 10 9) is inverted by 

multiple-state interactions when compared with the dependence shown in 
figure 5.8. The high probability of multiple-quantum transitions in this 
case contributes to the extremely fast and thermally insensitive relaxation 
rates measured in the halogens and it destroys the concept of a single 
"relaxation time" that is independent of the nonequilibrium state of the 
process for molecules of this type. 

T. K 
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Figure 5.9.- Multiple-quantum rate coefficients for Br 2 (m)-He predicted 
using the Kerner harmonic oscillator solution, equation (5.18), with 
L = 0.02 nm. 

5 . 6 Concluding Remarks 

We have relied on a collinear semiclassical model for vibrationally 
inelastic collisions entirely for pragmatic reasons. The collinear geometry 
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affords an economically reasonable means of estimating V-T rate coefficients 
for excited molecules and the semlclasslcal approximation Is easily reduced 
to practical analytic solutions. While these simplifications clearly 
obviate the quantitative accuracy of the calculations, no serious omission 
other than rotational coupling Is apparent that would modify their qualitative 
nature, even In the presence of uncertain interaction potentials. Unfortu- 
nately, an attempt to confirm the predicted features through experimental com- 
parison was Inconclusive. However, the experimental conditions that would 
test the model most severely can at least be Identified. For example, the 
choice of collision partner has a large Influence on the rate coefficient sen- 
sitivity to initial quantum number. Note the large deviations of k ,/mk, „ 
from unity in figure 5.5 for CO-Ar compared to those in figure 5.4 for Co-He. 
Furthermore, the increased oscillator distortion caused by heavy atom impact 
requires a more complete description of the interaction than needed for light 
atoms. From another viewpoint, the lesser sensitivity of some features of the 
prediction to uncertainties can guide the choice of experimental variables to 
be emphasized. In particular, an apparently universal feature of the V-T 
excited-state rate predictions is their monotonic low-order variation with 
quantum number. Once this feature is confirmed, the experimental emphasis 
can be shifted to the less predictable variations with kinetic temperature. 
Finally, a comparison of the estimates using various potential parameters 
suggests that a self-consistent set of experimental rates for both high and 
low initial quantum numbers contains much more information that defines the 
interaction potential than ground-state rates alone. 

Comparisons of the analytic and numerical rate coefficients graphically 
delineate the suitable range of application for each analytic model. However, 
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the utility of an analytic approximation can also depend on the physical 

properties of the application. For example, the Kerner harmonic oscillator 

model, with anharmonically corrected frequencies, predicts k /mk 

m,m-l 1,0 

with surprising accuracy for many molecules; but before the model can be 
economically applied, an analytical solution to the thermal averaging 
integral, including Kerner 's transition probability formula, awaits develop- 
ment. Even with that solution in hand, one must be concerned with the effect 
of anharmonicity for each molecule treated by the model. On the other hand, 
Mies’ solution for enharmonic oscillators, equation (5.15), fails at high 
temperature. At those conditions, however, many nonequilibrium processes 
are insensitive to the V-T rates of excited states, either because the vibra- 
tional state population distribution is nearly Boltzmann or because the 
process is controlled by some separate energy-transfer mechanism. At lower 
temperatures, the model accurately deals with a broader range of oscillators 
because anharmonicity is rigorously included. Collision partners for which 
the theory fails are poorly treated by all the analytic solutions based on a 
semiclassical approximation. Similarly, the frequently used formula 
developed by Keck and Carrier is useful because of its simplicity, but the 
additional computation required by the Mies solution is not prohibitive. 

The series function $(-g,A) converges rapidly and the matrix elements 
may be computed in advance. 

Calculations of multiple-quantum transition rates from excited states 
validate the assumption most often made in kinetic models of nonequilibrium 
processes: they can usually be neglected. As before, at very high tempera- 

tures where multiple-quantum transitions become competitive, a nonequilibrium 
process is usually not controlled by excited-state V-T rates, while 
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ground-state transitions are still dominated by single-quantum steps. 

Molecules with closely spaced vibrational energy levels, such as the halogens, 
are notable exceptions requiring a more careful analysis. 
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CHAPTER 6 


EFFECTS OF ROTATIONAL TRANSITIONS ON VIBRATIONAL ENERGY TRANSFER 

The preceding chapters have utilized a collinear collision model in which 
vibrational motion was the only form of energy transferred during a collision, 
and, indeed, in many cases it is. However, the use of a one-dimensional col- 
lision geometry renders the model incomplete in the sense that cross sections 
and rate coefficients cannot be obtained directly, but require artificial 
three-dimensional parameters such as steric factors or hard-sphere cross 
sections that must be estimated by some other means. 

From a more physical point of view, collisions with sufficient energy to 
induce vibrational transitions will simultaneously cause numerous transitions 
among the more closely spaced rotational states, thus Invoking an additional 
energy sink not represented in the collinear models. This chapter investi- 
gates the influence of rotational motion on the net rate of vibrational energy 
transfer. It is motivated by a need to assess the validity of the collinear 
models since they remain the most practical means for predicting vibrational 
rate coefficients in a kinetic analysis. With that motivation, emphasis is 
directed here toward the net rate of vibrational energy transfer summed o^er 
all final rotational states, rather than individual vibration-rotation transi- 
tional rates, since only the former can be compared with the collinear 
predictions. 

In the sections to follow, a semiclassical three-dimensional collision 
model is first developed in detail. It is followed by a description of some 
approximations that significantly reduce the number of coupled states neces- 
sary to obtain a complete solution, thus making numerical results practical. 
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The collision dynamics are then studied for molecular types that represent all 
extremes in the role of rotational motion on vibrational transition rates. 

The role of rotational coupling is summarized and categorized in conceptual 
terms in the final section and the common characteristics of all related V-R-T 
mechanisms are discussed. 

6.1 Vibration- Rotation Collision Model 

A three-dimensional, vibration-rotation collision model is formulated 
here within the semiclassical framework described in appendix B and based on 
the following underlying concepts: the incident particle is considered struc- 

tureless while the target is a diatomic heteronuclear molecule. A natural 
parameter for measuring the heteronuclear nature of the molecule is its mass 
ratio 

Y = vxj (mj^ + m^) 

where and m^ are the nuclear masses and m^ > mj^. The inertial proper- 
ties of the molecule are modeled by a Morse-oscillator/rigid-rotor description 
(discussed in ch. 3). The collision geometry is described classically 
with the associated coordinates shown in figure 3.4 as viewed in a rotating 
interaction plane containing all three nuclei. The classical trajectory is 
assumed to be dominated by the spherically symmetric component of the interac- 
tion potential and is determined from just those terms. Consequently, the 
relative path remains in a single plane and the subsequent formulation is 
greatly simplified. These concepts are reflected in an illustration of the 
collision geometry shown in figure 6.1, where it is viewed from a space-fixed 
position in a center-of-mass reference frame aligned with the plane of the 
trajectory. 
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from a fixed position in the center-of-mass reference frame. The x%y’ 
plane contains the molecular mass center and the incident particle path. 
Only the portion of the molecule above the plane is shown. 


The detailed formulation of a semiclassical model divides logically into 
three main parts — one defining the nature of the interaction potential 
between the colliding pair, one detailing the quantized motion of the molecule 
in response to a time-dependent disturbance, and one describing the classical 
motion of the incident particle that produces the disturbance. The descrip- 
tion to follow proceeds in the same order. 


6.1.1 Interaction Potential 

A comparison of the predictions from the collinear model described in 
chapter 4 with those of a three-dimensional model will be most meaningful if 
the interaction potentials are similar. We therefore make use of the expo- 
nentially repulsive and pairwise additive interaction discussed in chapter 3 
and write, with reference to figure 3.4, 

V’(x^,x^) = (6.1) 
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To conveniently separate the molecular dynamics from the classical motion, 
however, equation (6.1) must be expressed explicitly in terms of coordinates 
describing each of the motions separately. To this end, we have shown in 
chapter 3 that equation (6.1) can be represented to first-order in r/x by 

V'(5,r,S) - A .-*/!■ <r/W cos « ^ ^-(l-YHr/L)cos ^^ 2) 

where the relative angle, 6, is related to the space-fixed coordinate angles in 
figure 6.1 by 


cos 6 = sin 0 cos((|) - 0) (6.3) 

Equation (6.2) expresses the potential in the desired explicit form. However, 
to aid the following separation of dynamical equations into coordinates that 
describe the motion of each collision partner, we shall temporarily denote the 
potential by an equivalent notation V' (x,r,6) e V' (q,t) where q = q(r,0,(()) 
locates the molecule in configuration space and t emphasizes the temporal 
dependence of the trajectory coordinates, x(t) and fi(t), appearing in equa- 
tions (6.2) and (6.3). 


6.1.2 Quantized Molecular Dynamics 


In appendix B, we show that the quantum-mechanical equation of motion may 
be reduced to a set of differential equations — one for each bound eigenstate 
of thfe tthdlsturbed molecule — in terms of the probability amplitudes associated 


with each bound state. The result is 


ih 


dt 


= Yj Cn(t)e’'^“j “n)t<j|v»(^,t)|n> 


(6.4) 


n 

where the bracket notation refers to 


<j |V (q,t) |n> 



(q)V' (q,t)i|<jj(q)dq 


(6.5) 
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integrated over all q space and is an eigenfunction of the undisturbed 

molecule determined by the intramolecular potential, V^Cr). As in chapter 4, 
the intramolecular potential is modeled by the Morse function: ‘ 

Vi(r) - (6.6) 

Combined with a rigid- rotor description of the molecular rotational motion, 
the eigenfunctions are then determined by throe quantum numbers |vJlm> 

and have their eigenfrequency given by 

(6.7) 

The remaining task, leading to solutions for equation (6.4), is then to evalu- 
ate the matrix elements defined by equation (6.5) in terms of the fundamental 
collisional and molecular parameters. 

The numerical labor of solving equation (6.4) will be reduced greatly if 
the matrix elements can be factored into a time-dependent term obtained classi- 
cally and state-dependent terms containing all the quantum-mechanical spatial 
integrals. The spatial integration may then be completed independently and in 
advance for transitions between all eigenstates in the basis set. A step in 
that direction is taken by writing 

<j|v'(q,t)|n> = U(b,t)Vj^(t) (6.8) 

where, from the potential given by , equation (6.2), the right-hand terms are 

U(b,t) = A e“x(b,t)/L {.6.9) 

and 

a. Y cos 6 

V’ (t) = <j|e^^ ln> (6.10) 

^ i=l,2 

with Oj = Y and - -(1 - Y). Clearly, the function U(b,t) can be 
obtained classically, but the expression for Vjjj(t) must now be developed to 
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further isolate its time dependence. A useful aid will be the expansion of 
z cos 5 

e in terms of Legendre polynomials Pj(cos 6) given by equation (3.30), 

In the present notation, the expansion is 

J 

/ 

( 6 . 11 ) 


a. 7 “ cos 6 \ / 

‘ ^ ■ Z (t^) l) 

J=0 ^ f 


Pj(cos 6) 


where ij(a^r/L) is a spherical Bessel function of the first kind.^^ 

Recalling that |n> refers to a vibration-rotation state with quantum 
numbers |v£m> and that the molecule is represented as a Morse-oscillator/ 
rigid rotor, the eigenfunctions represented in equation (6.10) may be factored 
in the traditional manner according to 


The factorization, in conjunction with equation (6.11), then allows the matrix 
elements to be represented as a sum of products given by 




j=o 


(2J + 


v’ V 


(t) 


(6.13) 


where 



(6.14) 


is a radial matrix independent of time and . 

= <A’m' |Pj(cos 6)|jlm> (6.15) 

is a spherical matrix element containing the trajectory coordinate ?2(t) via 
equation (6.3). The remaining discussion is divided into separate develop- 
ments of each element. 
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6. 1.2.1 Radial Matrix Elements 


The radial term defined by equation (6.14) appears in expanded form as 

First note that when the molecule is homonuclear (i.e., 7 = 1/2), the radial 
term will be nonzero only for even values of J and that it will emphasize 
even J values for most heteronuclear molecules where Y is nearly 1/2. 

This property directly affects the probabilities of rotational transitions. 

To evaluate the spherical Bessel functions in equation (6.16) in terms of the 
molecular coordinate r, we note that its arguments lie typically in the range 
0 < a^r/L < 5. Larger values are suppressed by the vibrational eigenfunctions 
that approach zero at large r. In this range of arguments, the spherical 


Bessel functions may be calculated with rapid convergence by the ascending 


- — 


, . z^/2 (z^/2)^ 

(2J + 3) 2!(2J + 3)(2J +5) 


n (2j + 1) 
j=0 


(z^/2) 


■ ] 


The radial matrix element may then be written as 


^(J) ^ (2J + 1)Y^ ^ (Y^/2) ^ ^ 

v*v J n! n 


/i *\/ \ ^ n"4“*J 


^(2n+J) (6 18 ) 

V V 


n (2J + 1) n =0 
j=0 


n (2J + 2k + 1) 
k=0 


where represents the remaining simplified integral 




Equation (6.19) is similar to the integrals that define dipole matrix elements 
associated with vibrational band intensity calculations and, as such, it has 
been solved exactly for m = 1 and 2 using Morse oscillator 
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eigenfunctions In principle, the integral can be solved exactly for any 
Integer value of m, but the complexity of the solution increases and becomes 
impractical for powers larger than m = 2* An iterative numerical method has 
also been developed,^®® but it becomes laborious for large quantum numbers. 

In appendix C, we derive an approximate closed-form solution to equa- 
tion (6.19) with the result: 

-^ (^) Vv «•“) 

^ s=0 


where X = k' - 2 - (v' + v) , and F(X + s + 1) is a gamma func- 
tion. Coefficients Dg and contain elements of the Morse eigenfunc- 

tions (defined in appendix C). The derivation of equation (6.20) depends only 
on the provision that X » 1, which is easily obtained because k* >> 1 for 
all diatomic molecules. Typical values of k’ range from 37 for to 161 for 
CO. Corresponding errors in from equation (6.20) are 1.5 percent for 

and 0.3 percent for CO. Hence equation (6.20) is sufficiently accurate in 
this application, although some numerical difficulties arise that can limit its 
use. For example, the terms of the summation in equation (6.20) alternate in 
sign so that, for large k' and large vibrational quantum numbers, small dif- 
ferences between extremely large terms cause the loss of all significant digits. 
The CDC-7600 computer with 28 digits in double precision allowed equation (6.20) 
to be evaluated with at least 3 accurate digits for vibrational quantum num- 
bers less than 12 when k’ = 161 and for all quantum numbers to the continuum 
when k' = 37. In the few cases when equation (6.20) could not be used, the 
integral was evaluated with a standard Gauss-Laguerre quadrature algo- 

very high accuracy. (The author is indebted to D. G. Galant, 

NASA Ames Research Center, for developing the algorithm.) 
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6.1. 2.2 Spherical Matrix Elements 

The spherical term defined by equation (6.15) incorporates the rigid- 
rotor eigenfunction introduced in equation (6.12). The eigenfunctions are 
known to be spherical harmonics that satisfy^^^» 

L2Yj^^( 6,<|)) = h2£(Jl + 1)Y^^(0,(|)) 

and 


where L is the orbital angular momentum operator and is a component. 

Note that the Legendre polynomial in equation (6.15) can also be represented 
as a spherical harmonic according to 


allowing the spherical matrix element to be written as 

'•27T ^TT 
'0 •'0 




The solution to equation (6.21) follows common procedures in the mechanics of 
quantized angular momentum^®^ > and it is described in detail in appendix D. 
The result takes the form 




( 6 . 22 ) 


where the time- independent term is 




(j) 


* m^ £m 


. , ,'(J-a)/2 /(j + a)!(j - m )T ; (’■’ ■' ‘Y*’ •’ 


(6.23) 


Nonzero values are obtained only if 


m = m* 
J ± m 


even J 


(6.24) 
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The bracket symbols in equation (6.23) are Wlgner 3~j symbols that account 
for the vector coupling between the angular momentum states, |il’m’> and |)lm>, 
and the trajectory orbital momentum designated by J. They impose further 
constraints associated with nonzero values of ixa. listed col- 

lectively as 


- il| < J < £' + £ 


]m| < il , |m' I < ii' , |m| < J 

V + + J even 


(6.25) 


6.1.2. 3 Complete Factored Matrix Elements 

Factorization of the time-dependent terms in the complete matrix elements 
defined by equation (6.8) can now be accomplished by defining a time- 
independent matrix element (with an unprimed symbol) as 


V 

V 


I 


£'m' vilm 


^ , 


(2J 


V V 


*^£'m' £m 


(6.26) 


and writing the complete matrix element as 

<v'£'m|v’(q,t)|vJlin> = U(b, t)e"^ Vjl’in’vllni 

The time-dependent terms, U(b,t) and ^2(b,t), in equation (6.27; are then 
the functions required from the classical trajectory. Generally, they repre- 
sent both an induced force and an induced phase shift in the molecular motion, 
where the magnitude of the latter depends on the rotational transition con- 
sidered. However, the primary quantal properties of the collision dynamics are 
determined by the time- independent matrix elements given by equation (6.26). 

For example, the summation limits in equation (6.26) reflect the constraints 
imposed by equations (6.25) and lead to the only significant selection rule 
associated with vibration- rotation transitions. Recall from equation (6.16) 
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that* for homonuclear molecules, nonzero only when J is even. Thus, 

nonzero matrix elements and transition probabilities will '•'occur for homonuclear 
molecules only when - jlj is even, according to equation (6.26), Simi- 
larly, for heteronuclear molecules, transitions of even will domi- 

nate, although odd- increment quantum changes are allowed in those molecules. 

Other selection rules regarding intermultiplet transitions, m^ to m, are also 

« 

Inplled by equation (6.25) and their accompanying zeros may be observed in 
table 6.1, where a typical transition matrix is listed. However, the role of 
these degenerate states in determining the observable behavior of a kinetic 
process is usually lost due to subsequent averaging and hence become important 
only in the computational aspects of the collision model. 


TABLE 6.1- TIME- INDEPENDENT TRANSITION MATRIX FOR CO 

V = 0, a - 2 



2 

1 

0 



-1 

-2 



mm 






3 






2 






1 






0 




0 

-.065 

-1 

0 

-.097 

0 

.156 

0 

-2 

.022 

0 . 

-.135 

0 

.142 

-3 

0 

.034 

0 

-.189 

0 

-4 

-.006 

0 

.044 

0 

-.322 


Note: Shaded elements depict the symmetry of the matrix. 


Figures 6.2 and 6.3 illustrate some general properties of the time- 
independent matrix elements and allow one to reach some early conclusions con- 
cerning the nature of vibration- rotation transitions during a collision. For 
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e 6.2- Variation of 
tion (6.26) with an; 
v' =1, and L = 0.0, 


10 " 



> 







I 



A1 

(b) CO 


Figure 6.3.- Concluded. 

example, figure 6.2 depicts the variation of matrix elements with projection 
state quantum number. These features will be of Interest later when methods 
for reducing the number of coupled states in the calculation are sought. At 
this point, we simply note that transitions from the corners where m = ±il, 
m* ■ ±Z*, and m* , m have the same sign are the dominant route by which energy 
is transferred. As I and V increase, the dominance also increases so that 
In the limit of large £, A’, all other m ±Z and m' ^ ±£' states may be 
ignored. 

Another basic aspect of the matrix elements is illustrated in fig- 
ures 6.3(a) and (b) where for both H 2 and CO — two extremely opposite molecular 
types — only states with small differences in angular momentum, Ail (where 
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All E |£' - Jl|) are shown to share significant coupling. This feature persists 
for all values of Z, as shown by the sets of symbols in figures 6.3, and it 
also holds for all vibrational transitions. Note that the emphasis of small 
AZ transitions is independent of any degree of resonance that may occur between 
the initial and final states of the transition since equation (6.26), from 
which these matrix elements are obtained, contains no reference to the eigen- 
energies of the transition. 

Nevertheless, resonance enhances the probability for transition. Its 
effect, in combination with the preceding small AZ constraint may be 
anticipated by examining the relative vibration- rotation eigenenergies 
depicted in figure 6.4 for both H 2 and CO. Note that while the transition 
vZ Z^ = 0,8 1,0 in para-H^ is nearly resonant, AZ = 8 is too large and 

the coupling between these states (fig. 6.3(a)) will be very small. Much 



Figure 6.4.- Vibration-rotation eigenenergies for para-H^ and CO. 
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(b) CO (NOTE: Not all rotational levels are shown). 

Figure 6.4- Concluded. 

larger will be required in CO to approach resonant transitions while the 

restriction to small A£, will make the rotational aspect of a vibrational 
transition insignificant. From these observations, we may conclude that, while 
transition probabilities are enhanced by resonance, vibrat'ion-^'Potat'ion tvansi- 
tions with smalt A£ witt dominate the intramoteoular energy transfer process^ 
regardless of resonance! 

Finally, we see from figures 6.3 that, although odd A£ transitions do 
not occur in homonuclear molecules where Y = 1/2, a slightly heteronuclear 
molecule like CO (Y = 0.43) will allow odd A£ transitions with only moderate 
suppression, as indicated by figure 6.3(b). Hence, the selection rule regard- 
ing even A£ transitions applies strictly to homonuclear molecules. 

The foregoing derivation of the molecular equations of motion has explic- 
itly identified the required trajectory functions. We can now proceed to 
evaluate those functions classically in terms of the fundamental collision 
parameters. 
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6.1.3 Classical Trajectory 


The classical, two-body, central- force equations of motion are obtained 


In appendix B as 



dS2 _ ^ 

dt ~ -2 


(6.28) 

(6.29) 


where the coordinates refer to figure 6.1 and V(x) is a spherically symmetric 
version of the interaction potential. The total energy, E, and the correspond- 
ing initial speed, u = /2E/y, are "effective" values averaged over the trajec- 
tory (as discussed in ch. 4). 

In this chapter, we shall not include the dynamic influence of the mole- 
cule on the incident particle motion as in chapter 4. To do so in three 
dimensions would require a great deal of artificial approximation and would 
demand extensive comparisons with exact quantum-mechanical solutions for vali- 
dation, as in chapter 4. Such an investigation is outside the intent of this 
chapter. Instead, we use the results in chapter 4 as a guide to those colli- 
sion parameters where the dynamic coupling is unimportant and confine the 
examples used here to only those cases. 

The form of the interaction potential necessary to determine a trajectory 
is obtained here by spherically averaging the model potential given in equa- 
tion (6.2) over all coordinates of the molecule in its initial state. Thus, 

V(x) = <ilv'(q,t)|i> (6.30) 

where |i) denotes the initial state. By use of equation (6.27), the matrix 
element above may be rewritten in the more workable form: 

V(x) = A (6.31) 
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This method of averaging the potential produces a single trajectory for each 
initial state of the molecule and it is independent of the numerous possible 
final states. Transition probabilities computed using other methods of aver- 
aging (e.g. , involving a final state designation) were found to be only negli- 
gibly different from those obtained with equation (6.31) when the energy trans- 
ferred inelastically is small compared to the total energy. Eqiaatlon (6.30) is 
an approximation consistent with the adopted classical equations of motion 
since their accuracy is likewise contingent on the requirement that the energy 
traded inelastically remain small. 

At this point, we have the sufficient formulation to obtain a full 
numerical solution of equation (6.4). However, the numerical Integration of 
a coupled set of differential equations often proceeds with much less labor if 
all equations relax at intrinsically similar rates. The selection of step 
intervals in the independent variable (t in this case) is then controlled by 
the unanimous behavior of all dependent variables rather than the conflicting 
behavior of several dependent-variable subsets. (The amplitudes, c^(t), and 
the trajectory coordinates, x(t) and f2(t), are conflicting subsets with dissinh- 
ilar relaxation rates in this case.) While a set of equations cannot always be 
idealized in such a manner, the numerical effort is reduced here significantly 
by using approximate analytic solutions of the classical trajectory equations 
and solving only the molecular equations of motion numerically. 

An analytic solution of the trajectory equations is possible only because 
we haye adopted an exponential interaction potential with convenient analytical 
properties. In most cases, however, nonexponential potential functions, such 
as the Lennard-Jones potential, eqxiation (3.22), can be represented by an 
exponential function over the essential regions of interaction with acceptable 
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accuracy. A development of analytic solutions to the trajectory equations is 
described in detail in appendix E. The solutions are based on first-order 
approximations developed by Hansen and Pearson^ and later extended by 
Stallcop.^^^ Here we outline the procedure only briefly as it has been 
.adapted in this application to obtain the functions U(b,t) and fi(b,t). 

The function U(b,t) defined by equation (6.9) may first be converted to 
a more convenient form by noting that at closest approach t = 0, x = x^, and 
dx/dt = 0. Equation (6.28) may then be solved for V(x^) to obtain 

A v^. = E[1 - (b/xo)2] (6.32) 

Equation (6.9) then becomes 

U(b,t) [1 - (b/x„)2] (6.33) 

ii 

When b = 0, we also have the solution to equation (6.28) 

(6.M) 

Hansen and Pearson^ assumed that a solution to equation (6.28) for nonzero 
impact parameters will have a similar form and defined a function, a|j(b,t), so 
that, for all values df b, 

e ^o)/^ = sech^ [a^ (b, t)ut/2L] (6.35) 

We show in appendix E that an expansion of both sides of equation (6.35) 
about t = 0 then gives, to first-order, 

a^(b,0) = [1 - (b/Xo)2(l - (6.36) 

Furthermore, exact numerical solutions of equation (6.28) show that a|^(b,t) 
varies so slowly with t that equation (6.36) may be used for all t. Fig- 
ure 6.5 illustrates the accuracy of equation (6.36) for large impact param- 
eters where the error is greatest and at a collision energy near threshold. 
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Figure 6.5.- Exact and approximate (eq. (6.33)) trajectory functions for non- 
zero impact parameters; E/AV^^ = 10"^. 


A similar figure for energies 10 times greater would appear almost identical. 
Thus, we have the required analytic approximation: 

U(b,t) [1 - (b/Xo)2]sech2[a^(b,0)ut/2L] (6.37) 

where a|^(b,0) is obtained from equation (6.36) and x^ is the positive root 
of equation (6.32). For small impact parameters where the peak contribution 
to vibrational transition cross sections is made, x^ may be approximated by 
the analytic expression obtained from equation (6.32) for b = 0; namely 




Q ^ Xo(b = 0) = L In 




(6.38) 


Figure 6.6 Indicates the range of b/L where equation (6.38) may be applied. 
Finally, for very large b/L where b/x^ 1, the classical path follows a 
straight line, as figure 6.6 indicates, and although equations (6.36) 
and (6.37) are still sufficiently accurate, an alternate and more accurate 
approximation is 
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U(b,t) 



A exp [-/(b/L)2 - (ut/L)2j 


(6.39) 



Figure 6.6.- Distances of closest approach. 


The remaining trajectory function to be determined is Q(b,t). Hansen 
and Pearson^®® obtain an analytic expression by expanding ?2(b,t) in a Taylor 
series about t - 0, again where the interaction is greatest* and keeping only 
terms to first order. Thus, 


Q(b,t) = jj(b,0) + t 



+ . . . 


With the aid of equation (6.29), a first order approximation is then 

n(b,t)«^t (6.40) 

Stallcop^^^ explores the error in neglecting higher-order terms and finds it 
to increase with A£. However, the error in cross section is only about 
10 percent for AH = 4 and appears insensitive to collision energy. We see in 
figure 6.3 that the matrix elements and hence cross sections decrease rapidly 
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with AJi and will be of marginal interest for A£ > 4. Thus, the approximate 
function given by equation (6*40) is used for its computational simplicity. 


6.2 Representations of the Results 

The solutions to equation (6.4) provide complex probability amplitudes, 
c^(+»), for all states included in the vibration- rotation set. The final 
transition probability from an initial state |i> to some other state |n> 
is then 




In more explicit notation, we designate the initial state by unprimed quantum 
nuinbers, vilm, and denote the probability 

(6.41) 


P(E,b) = |c 

0 V X. m 


vJinmr’ Jl’m’ 

However, transitions between individual projection states, m and m' , are 
rarely of interest in the analysis of a kinetic process. A more easily mea- 
sured and useful quantity may be referred to as the vibvation--votation transi- 
tion probability, defined as the final probability averaged over all initial 
projection states and summed over all final projection states according to 


P(E,b) = o n I 1 y] P(E,b) (6.42) 

vil'^v'Jl’ m,m’ v£nr>v* £^m’ 

Finally, to compare with a collinear model in which rotational transitions are 
nonexistent, we define a net vihvational transition probability as the total 
probability of finding the molecule anywhere in the manifold of rotational 
states associated with a designated vibrational state. Such a probability is 
computed simply as the sum 

P(E,b) = 2 P(E,b) (6.43) 

vS.-^v' 
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The net vibrational transition probability so defined still depends on the 
initial rotational state, however, and cannot be compared directly with col- 
linear predictions. To evaluate the collinear approximation, we must first 
identify the properties of that must be obtained to produce similar 

cross sections and rate coefficients from either collision model. 

The relation of cross sections to transition probabilities has been dis- 
cussed in chapter 2 and is known to be 

P(E,b)b db (6.44) 

o v£-»-v' 

Equation (6.44) poses no difficulty in a three-dimensional model but, as 

chapter 5 demonstrates, its application to the collinear model requires some 

additional consideration. Since the Impact parameter, b, is not a collinear 

variable, an effective "hard-sphere" cross section, Oq, must first be chosen 

by some independent means and the collinear cross section is then defined by 

0 ( 1 ) = o^Pd) (6.45) 

v>v' v>v' 

The validity of equation (6.45) is contingent on the idea that o^ is 
invariant with both quantum number and collision energy. One test of the 
utility of a collinear model is then to test this contingency by defining a 
parameter equivalent to Oq but obtained from the three-dimensional collision 
theory. To that end, equation (6.44) may be recast into a form defining an 
"equivalent elastic cross section," a®, by the relation 

a(E) = a®P(E,0) (6.46) 

v)l->-v* ° v£->-v' 

where 



(6.47) 
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and is an arbitrary collision radius. Equation (6.47) has been written 

in a manner that also correlates with other concepts used in the methods for 
selecting a hard- sphere cross section, For example, the constant, x^, may 

be considered a radius corresponding to the elastic collision radius used to 
compute gas-kinetic collision rates.® The associated elastic cross section is 
then the constant term, appearing in equation (6.47). When a Lennard- 

Jones interaction potential is adopted, x^ is often equated to the zero- 
potential radius and evaluated from viscosity or viral-coefficient 
measurements • ® ® 

The remaining integral term in equation (6.47) is 


S = 2 


f rp(E,b)/p(i,o) 

•^o ! v£*>v’ J \ c/ \^c/ 


(6.48) 


which may be identified with a ”steric factor" often used as a correction for 

three-dimensionality with collinear collision models. The steric factor must 

be considered invariant in the collinear model, although its magnitude here 

depends on the values of 1 and the transition quantum numbers. Hence, cases 

0 

in which either or S are invariant with both quantum number and colli- 

sion energy show possibility as examples where a collinear collision model is 
applicable. However, the invariance of or S alone is not sufficient to 

declare collinear rate-coefficient predictions valid, as shown in the follow- 
ing discussion. 

According to the results in chapters 4 and 5, a semiclassically deter- 
mined cross section is related to the thermally averaged rate coefficient, in 
the notation of this chapter, by 

e 


k(T) 

v£-^v’ 


C E e j crp) 

a' v£->v' s, ' / 
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where 


e » E/kT , ^ * (SkT/irv)^/^ 

If we now assume that transitions for which the steric factor defined by 
equation (6.48) is invariant are those for which ~ 

contribution of rotational frequencies to the eigenfrequency can be neglected) ^ 
then equation (6.49) is well approximated by 

e I 

k(T) « e ^ I P(e, 0 )(e + |e ,|)e ^ de (6.50) 

•'o v£->v* ^ 

To compare with the collinear rate coefficient that presumably represents the 
average of all initial rotational states, we introduce a rotational state 
population fraction N^j^(T^) that is dependent on a rotational temperature, 

T^. The net vibrational rate coefficient is then approximately 

k(T,Tp « C a® e f \Y, N (T )P(i,0)| (e + |e^, l)e~® de (6.51) 

For cases where T = T^, it may be compared directly with a collinear rate 

coefficient. The additional conditions to be met before a collinear model may 

be applied are then shown by equation (6.51) to be that P(E,0) must be rela- 

vi-^’ ^ 

tively invariant with qwxntum nvirher^ £3 and it must vary with E in the same 
manner as the aotlinear equivalent probability. 

6.3 Numerical Methods of Solution 

Numerical solutions to the set of differential equations (6.4) were 
accomplished using the same algorithm that was applied to the collinear model 
in chapter 4 , namely, a fifth-order polynomial extrapolation technique 
developed by Bulirsch and Stoer^^ and provided in FORTRAN by Gear.^^ The 
method is shown by Hull et aZ. to be generally advantageous both in 
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computer expense and reliability (error) over several other established 
methods, particularly the frequently used Runge-Kutta methods. A typical col- 
lision event could be computed in this application with as few as 50 steps 
while maintaining six-digit accuracy. Most solutions were obtained with 
100 steps and initialized at a distance where the interaction potential was 
10”*^ of its value at closest approach. 

An additional consideration required by the three-dimensional model was 
’ the' economic use of computer storage. The time- independent matrix elements 
V « , given by equation (6.26), were computed prior to the numerical 

V ^m v36m 

integration and stored in memory. However, their designation by a sixfold , 
index set suggested by the associated quantum numbers is impractical even for 
the CDC-7600 computer with its large-core memory. Thus, advantage was taken 
of the zeros and symmetry of the matrix elements (see tablje 6.1 for examples) 
and an index scheme was developed that reduced the sixfold^ matrix to a one- 
dimensional vector with no zeros or duplicate elements. Details of the index 
method and its implementation in the computational procedure are described in 
appendix F. 

An important conclusion, quickly recognized from early solutions, was 
that the number of coupled vibration-rotation states required to achieve a 
convergent solution (i.e., a solution for which the further addition of states 
made no change) was too large to be computed in a defensible time. For 
example, the execution time per step on the CDC-7600 computer was approxi- 
mately N^/2 msec, where N is the total number of coupled states. A con- 
vergent solution usually requires that at least all energetically accessible 
states must be Included. Thus, if I denotes the uppexrmost rotational state, 
there are 2JI+1 projection states for each orbital state £ and hence (£+1)^ 
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total states from Z = 0 to £. Computer time then varies as (£+1)**. For 
collision energies near vibrational threshold, E > so that, from fig- 

ure 6.4, minimum values are £ - 8 for H 2 and £ - 33 for molecules like CO. 

A corresponding minimum computing time for collisions is then approximately 
5 min — an acceptable value that has allowed numerous studies of vibration- 
rotation transition rates in H 2 * However, for molecules like CO where the 
rotational frequency is much less than the vibrational frequency, we could 
expect to require more than 18 hr per case! Obviously, to study such mole- 
cules, we must seek ways to reduce the required number of coupled states. The 
greatest reduction will be achieved by any method for decoupling the projec- 
tion states within a given orbital state and averaging their effects in 
advance of the calculation. Several such methods are discussed in the follow- 
ing section. Their implementation has been a key factor in reaching the 
objectives presented here. 

6.4 Effective Hamiltonian and Other Approximations 

During the course of this study, three primary methods of approximation 
were examined. They are discussed in this section in an order of increasing 
utility to this work. 

6.4.1 Sudden Rotation/Perturbed Vibration Approximation 

In a recent analysis of vibration-rotation coupling in harmonic oscilla- 
tions, Stallcop^^^ drew renewed attention to the concept that rotational and 
vibrational motion can be treated in separate limiting approximations. For 
example, the rotational period, 1/B^, is typically very long compared with the 
collision period, = 2L/u. Consequently, the molecule appears rotationally 
stationary during the collision period and the induced mixing of rotational 
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states is predicted accurately by the "sudden" or "impact" approximation. ^ ^ ® 
The sudden approximation then provides a closed-form integral description of 
the mixed rotational state of the molecule at any time during the encounter. 

In contrast, the vibrational motion typically undergoes several oscillations 
while the incident particle is at close range. Stallcop treats the vibra- 
tional motion in the extreme adiabatic limit for which many oscillations must 
occur during the interaction period, but that limit is too restrictive for our 
purposes. Instead, we may recall, from chapter 4, the broad range of colli- 
sion energies in which a first-order perturbation treatment is successful and 
then treat the vibrational motion accordingly. Dynamically, the vibrational 
and rotational motions are thus decoupled and a complete closed-form integral 
description of the vibration-rotation transition probability may be obtained. 
Unfortunately, all series solutions of the resulting integral equation were 
found to be ill-conditioned and hence not calculable with the significant 
digits carried by available computers. Numerous attempts to restructure the 
formulation or to evaluate the integral by numerical quadrature were also 
unsuccessful for similar reasons. Thus, while the approach is mentioned here 
because of its potential significance as a means of analyzing vibration- 
rotation energy transfer, we were forced to abandon it for the present study. 

6.4.2 Maximum Coupling Approximation 

The broad range of magnitudes covered by the matrix elements, 
for the range of m and m* values suggests that the number of coupled states 
may be reduced by including only the dominant paths of energy transfer. As 
figure 6.2 indicates, the coupling between vibration-rotation states will 
be dominated for large Z by the projection states in which m and m* are 
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maximum and of the same sign. The "maximum coupling" approximation therefore 
simply excludes all states except those contributing to V .. and 

number of coupled rotational states In each vibrational 
set is then reduced from (£+1)^ to 2£ + !• Calculations incorporating the 
foregoing exclusions are discussed in the following sections where they are 
compared with the "effective Hamiltonian" approximation described next. 

6.4.3 Effective Hamiltonian 

An approximate method for decoupling and subsequently averaging the con- 
tributions of all projection states of each orbital state has recently been 
developed by Rabitz. Unlike the preceding two approximations, however, even 
the relative range of parameters for which the method is expected to be accu- 
rate has not been defined in terms of the molecular properties or collision 
parameters nor is a method of defining them apparent. Hence, for the present, 
the approximation must simply be tested by comparison with limited exact cases. 

In effect, Rabitz determines the form of the Hamiltonian required to 
exactly nullify the contribution of individual m and m' states in the matrix 
elements. He begins by assuming an interaction potential with the general 
form ' 1 : r ' : ' : 

00 

V'(x,r,6) = 2 u^(r,x)P^(cos 6) (6.52) 

J=0 

and finally obtains the effective matrix element, analogous to equation (6.5) 
or (6.27) as 
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(6.53) 


<v’i’ lv4(X.r,6) |vA> » + i) (2Jl + x 

00 

2 <v' hT(i^»5c) |v>(2J + l)-l/2 
J-0 ■ 

where is ah "effective potential" corresponding to equation (6.52). 

Equation (6.53) may be converted to terms corresponding to the particular 
potential, given by equation (6.2), through the use of equation (6.11). An 
equivalence is easily identified as 

\J 

Vj(r,i) - A e-='^*'(2J + 1) ij(“i f) 


U' J £> 
lO 0 0> 


so that 


<v' |uj(r,x) |v> = U(b,t)(2J + l)R^f^ 


The resulting "effective" time-independent matrix element is then 


V® 

v'i’vJl 


(-1) + l)(2i + 1)]!/'+ 


v+i .-V /£’ J a\ 

,Z , ^^2TTT ( j 

J=|il'-£| 0 0/ 

and the complete matrix element, analogous to equation (6.27), is 


(6.54) 


<v*&*|v;lv^> = U(b,t)V®,j^,^^ (6.55) 

The use of equation (6.55) in place of equation (6.27) reduces the total num- 
ber of coupled rotational states in each vibrational manifold from (I+l)^ to £. 

Note by comparison with equation (6.27) that the induced phase shift 
associated with the exponential argument, (m* - m)fi(t), in equation (6.27) is 
lost in the effective Hamiltonian approximation and the remaining formulation 
is independent of fl(t). We shall see in the discussions to follow, however, 
that the rotational energy transfer is dominated by small A£ and hence small 
values of m’ - m. Thus, even in a complete solution, the phase shifts 
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associated with J2(t) are subdued and the accuracy of the effective Hamiltonian 
approximation is not threatened by their neglect. 

6.5 Aspects of Convergence 

Before applying the preceding three-dimensional model, we first require a 
criterion for choosing a sufficient set of rotational and vibrational eigen- 
states to ensure convergence. Experience with the collinear model (ch. 4) 
has shown that near threshold only a few states higher than the energetically 
accessible vibrational states are necessary to obtain convergence. Thus, for 
ground-state molecules with E/ho)^ near or slightly greater than unity, only 
three or four vibrational states are often adequate. Similar concepts can be 
applied to the vibrational states in a three-dimensional model. However, the 
fundamental vibrational frequency of a diatomic molecule is always larger than 
its rotational frequency so that, in a three-dimensional model, collisions 
with sufficient energy to induce vibrational transitions (i.e., E/hWg > 1) 
always couple numerous rotational states in each vibrational manifold. The 
best manner of selecting the minimum, and yet sufficient, number of coupled 
rotational states for each vibrational eigenenergy is therefore difficult to 
determine in advance of a calculation. A seemingly logical first choice would 
be to include only the energetically accessible states and exclude all others. 
For E/h( 0 g, “ 1, this method of selection envelops many rotational states in the 
initial vibrational manifold and only a few at the next higher vibrational 
eigenenergy. We shall find, however, that this criterion for selection is not 
only inadequate to ensure convergence but with a semiclassical formulation it 
almost certainly guarantees a nonconvergent solution! The primary purpose of 
this section is therefore to examine the requirements of convergence in 
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detail. The reader should be aware at the outset, however, that the results 
pertain to our semiclassical model without energy conservation. Their rela- 
tionship to the convergence requirements for an energy-conserving collision 
model will only be inferred. 

We test our ability to obtain convergent solutions for the molecular types 
of interest by choosing the worst numerical case in an example (namely, CO-He 
collisions in which a maximum number of basis states is required) and contrast- 
ing the results with those for H^-He collisions, the opposite extreme. 

As figure 6.4(b) illustrates, a basis set that is convincingly convergent 
for CO-He collisions at E/hw^ ^ 1 will include an impractical number of 
states that extends to large To- deal with such situations, we must there- 

fore first determine the most appropriate method of approximation from sec- 
tion 6.4 that will reduce the necessary number of coupled states and that 
allows solutions for large Z to be obtained in a practical computing time. 

6.5.1 Evaluation of the Projection-State Decoupling Approximations 

In keeping with the predominant objective of this chapter to study the 
effects of rotational transitions on the net vibrational transition rates, we 
use as a basis for comparison both the probability, (defined by 

eq. (6.43)), and its component, defined by 

" 0^ I ^v£it^>v’ Jl*m’ (6.56) 

,m’ 

Note that the former is averaged over all initial m states while the latter 
pertains to a specific initial m state. However, both are summed over all 
final £* and m’ states in the vibrational manifold, v* , and are therefore 
both net vibrational transition probabilities. 
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Since the approximations described in section 6,4 are methods for decou- 
pling the influence of the projection states on the molecular dynamics and sub- 
sequently averaging their effects, the sensitivity of complete solutions to the 
quantum number of the initial projection state is an important aspect in under- 
standing of the relationship of the complete and approximate calculations. 
Figure 6.7 shows that, for CO-He collisions, relatively insensitive 

to the initial value of m, even for large impact parameters where the phase 
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b/L=4.5 
b/L = 9 


(a) A = 6 MeV 
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b/L = 0 
b/L = 6 


(b) A = 1000 eV 


Figure 6.7.- Distribution of net vibrational transition probabilities over the 
range of initial projection states for C0(v£ = 0,3)-He collisions at 
several impact parameters. Only^probabilities for v’ = 1 are shown. 
Fixed collision parameters are E/ha)^ = 1.08, L = 0.02 nm. The basis set 
contained £ = 0-8 rotational states in each vibrational manifold and 
included all corresponding projection states. Symbols denote results 
from: • ■ • a complete solution, o o o the maximum coupling approxima- 

tion. Effective Hamiltonian results appear as a single value. 


shift associated with fi(b,t) is greatest. This result occurs because no 
other dynamical phase interference exists between m states of the same 
eigenfrequency while ?2(b,t) is generally small in the primary region 

of interaction = 0 at t = 0 where the interaction is greatest). In view of 
the foregoing insensitivity to initial m, vibrational transition probabilities 
calculated with an effective Hamiltonian approximate the complete solutions 
more accurately than the maximum coupling approximation. Its similar accuracy 
in reproducing related net vibrational cross sections is indicated by 
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figure 6.8, where the variations of with impact parameter, calculated 

by the two approximations, are compared with complete solutions. Clearly, the 
effective Hamiltonian approximation is superior, particularly for the smaller 
interaction scale factor. A, shown in figure 6.8(b). (The smaller scale fac- 
tor is believed to be more realistic.®^) Similar accuracy is obtained for 
H 2 “He collisions as shown in figure 6.9. One should note, however, that the 
impressive accuracy of the effective Hamiltonian is aided significantly by our 
use of the net probability, » as a basis for comparison. Similar com- 

parisons for detailed vibration-rotation transitions would not appear as 
favorable. Others^^ » ^ ^ ^ ^ have obtained equivalent results for a variety of 
molecular types. 




b/L 

(b) A = 1000 eV 


Figure 6.8.- Variation of net vibrational transition probability with impact 
parameter for C0(vJ2, = 0,3)-He collisions using several approximate solu- 
tions. Only V* = 1 probabilities are shown. The basis set included 
£ = 0-8 in each vibrational manifold for all cases; E/hw^ = 1.08, 

L = 0.02 nm. Symbols denote results from • a complete solution, o the 
maximum coupling approximation, and ^ the effective Hamiltonian 
approximation . 
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EFFECTIVE 




b/L 

Figure 6.9.- Variation of net vibrational_transition probability for 
H 2 (v£ = 0,2)-He collisions; v* = 1, E/hw^ =1.1, A = 303 eV, and 
L = 0.0273 nm. The basis set included £ = 0-10 in both vibrational 
manifolds. Symbols denote the results from: • a complete solution, 

and the effective Hamiltonian approximation. 


6.5.2 Convergence Requirements for Vibration-Rotation Energy Transfer 

The individual vibration-rotation transition probabilities from a com- 
plete solution are shown in figure 6.10 for CO-He collisions at an initial 
kinetic energy just above the vibrational threshold. Rotational states from 
£ = 0 to 10 were included equally in vibrational manifolds, v = 0 and 1. The 
accompanying 484 differential equations and 14,883 dissimilar matrix elements 
exceed a practical upper limit for repetitive computation. And yet, compari- 
son to a similar calculation with £ = 0 to 8 rotational states in each vibra- 
tional manifold shows that convergence in the vibration-rotation probabilities 
^v£^v'£’ from realized. However, this result is not surprising 

because, while all open channels (energetically accessible states) have been 
included in v' = 1, most open channels in v’ = 0 are missing (see 
fig. 6.4(b) for reference). The interesting result is shown in figure 6.11 
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Figure 6.10.- Effect of additional rotational states on the vibration-rotation 
transition probabilities in C0(v£m = 0,0,0)-He collisions. Both cases 
were obtained from a complete solution including all projection states in 
the basis set and vibrational states v = 0, 1. An equal number of rota- 
tional states were included in each vibrational manifold. Collision 
parameters are E/hm^ = 1.08, b = 0. The interaction potential is 
defined by A = 6 MeV, L = 0.02 nm. 
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VIBRATIONAL MANIFOLD 
(a) CO(v£m = 0,0,0)-He 



NUMBER OF ROTATIONAL STATES IN EACH 
VIBRATIONAL MANIFOLD 

(b) CO(v£m » 0,3,3)-He 


Figure 6.11.- Convergence of the net vibrational transition probability for CO-He collisions (l.e., the 
probability has been summed over all I.', m’ according to eq. (6.43), but for a single 
initial m state) . A duplicate set of rotational states was included in each vibrational mani- 
fold. Only probabilities for v' = 1 are shown. Collision parameters are the same as in 
figure 6.8. Symbols denote results from: •, a conq)lete solution, o, the maximum coupling approxi- 

mation, the effective Hamiltonian approximation. 



where we imply from the behavior of value of m that, even 

before all open channels in any vibrational manifold are reached, a convergence 
in the net vibrational transition probability y , is obtained! Similar 

behavior is found for H 2 “He collisions, as figure 6.12 indicates. Note that, 
in this latter case, calculations including all open channels in both vibra- 
tional manifolds were possible and the convergence asymptote is convincingly 
unique. Figures 6.11 and 6.12 also include results from the approximate 
methods of solution and show them to approach similar asymptotes. 
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Figure 6.12.- Convergence of the net vibrational transition probability for 
H 2 (v£m = 0,4,4)-He collisions. Stipulations and symbols are the same 
as in figure 6.9. Collisions parameters are E/ho)^ = 1.5, b = 0, 

A = 303 eV, L = 0.0273 nm. Note that, in this example, cases are shown 
where all open channels are included in both vibrational manifolds. 


The uniqueness of the asymptotes obtained in CO-He solutions is not yet 
confirmed, however. In fact, subsequent calculations using the effective 
Hamiltonian for large but unequal numbers of rotational states in each vibra- 
tional manifold are shown in figure 6.13 to produce drastically different 
results. Note that, even though cases including £ = 0 to 40 in v = 0 but 
only £ = 0 to 10 in v = 2 contain all open channels in both v = 0 and 1, 
the transition probability is two orders of magnitude from a final convergence. 
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Figure 6.13.- Convergence of the net vibrational transition probability in 

effective Hamiltonian solutions for CO(vS. = 0,3)-He collisions v’ = 1, 
E/hoj^ = 1.08, A = 1000 eV, L = 0.02 nm. Cases are represented that 
include large but unequal numbers of rotational states in each vibra- 
tional manifold. Note that cases are also shown where all open channels 
are included in both vibrational manifolds. 


As more rotational states are added to v = 1, a final definite convergence is 
eventually reached where further additions of any kind have no effect. Con- 
versely, figure 6.13 also demonstrates that cases including a duplicate set of 
rotational states in each vibrational manifold obtain a solution near the 
final convergent value with relatively few rotational states. Recall, however, 
that the convergence criterion thus implied pertains only to the net vibra- 
tional transition probability, , and not the individual vibration- 

rotation probabilities 

While there is no attempt made here to construct a mathematically 
definitive argument showing why a duplicate but nonconvergent set of rota- 
tional states in each vibration manifold produces nearly convergent vibra- 
tional transition probabilities, the following conceptual explanation is 
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offered. First, recall that each time- in dependent matrix element, V , . , 

is a measure of the coupling strength between two vibration-rotation states. 

The matrix elements are shown in figure 6.3 to generally emphasize coupling 
only between states with small differences in angular momentum, regardless of 
the angular momentum of either state. This emphasis is primarily a consequence 
of the integrated overlap between rotational eigenfunctions of the undisturbed 
molecule and it is qualitatively unaffected by any vibrational change of 
state. As a result, vibrational transitions also occur predominantly between 
states with small differences in angular momentum. We shall refer to such 
pairs of states as "companion states" to imply that their angular momentum is 
similar (AH is small) but each is a member of a different vibrational mani- 
fold. Thus, if a rotational state has no companion in the basis set of an 
adjacent vibrational manifold, direct vibrational transitions from that rota- 
tional state will be improbable. 

With the foregoing general property of vibration-rotation transitions in 
mind, we can now describe its effect on the energy transfer process within the 
average molecule during a collision. To do so, we shall adopt a point of view 
compatible with our time-dependent semiclassical model and refer to the time- 
variant amplitude modulus instantaneous "occupation" of 

state |v£m>. By "average molecule" we then mean that the progression of state 
occupations in time provides a trace of the average path of energy flux within 
many identical molecules, all experiencing identical collisions. 

Since we consider collision energies mainly above the vibrational thresh- 
old, many rotational states are typically accessible. Early in the collision, 
the energetically accessible rotational states in the vibrational manifold 
containing the initial state become occupied, all with similar probabilities. 
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Near closest approach, vibrational transitions then begin between companion 
states and a corresponding occupation distribution develops in the adjacent 
vibrational manifolds. However, if the manifold receiving energy has been 
given fewer rotational states in its basis set (as it might because a higher 
manifold would have fewer open channels), then not all rotational states in 
the initial-state manifold have companions. For the energy in those rotational 
states to become available for a vibrational transition, it must first return 
to a rotational state in the same manifold where a companion exists. Obviously, 
such occurrences are the artificial consequence of an incomplete basis set and 
they stress the importance of including at least a duplicate rotational set in 
each vibrational manifold (if not a convergent set) to properly reproduce the 
transient dynamics during a collision. 

We may now ask why a duplicate but incomplete rotational basis set in 
each vibrational manifold approximates the convergent solution. Recall that 
it does so only for the total occupation of each vibrational manifold, that is, 
the occupation summed over all rotational states. If all rotational states 
have companions as they do in duplicate sets, the most probable paths are 
available for vibrational energy transfer from each rotational state consid- 
ered and hence no artificial impediment to the energy flow is introduced any- 
where, However, obtaining the correct net rate of vibrational energy transfer 
remains to be questioned. The net rate of energy transfer to all states of a 
vibrational manifold may be written as a sum proportional to the occupation 
distribution among its rotational states. At collision energies exceeding 
the vibrational threshold, the rotational states are strongly coupled 
and the depletion of any rotational state by a vibrational transition 
is rapidly restored, thus maintaining the distribution of rotational state 
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occupations in a time scale short cotiipared to the encounter period. When 
too few rotational states are included in the basis setj the occupation 
distribution is unnaturally constrained and vibrational transitions from ' 
each rotational state occur at an increased rate but from a fewer number of 
states., Since all rotational states are closely coupled, no particular 
preference is given to vibrational transitions from any of them and the net 
vibrational transition rate is only weakly affected, as we have observed. 

Note, however^ that the occupation distribution of rotational state's is 
strongly affected by the completeness of the basis set so that convergence 
in the individual vibration-rotation transition rates cannot be expected from 
an incomplete basis set, as we have also observed. 

Thus, we have rationalized, in the foregoing explanation, the reasons for 
expecting duplicate but nonconvergent sets of rotational states in each vibra- 
tional manifold to closely reproduce the net vibrational dynamics of the mole- 
cule given by a convergent set. We shall incorporate this convergence criter- 
ion along with use of the effective Hamiltonian approximation in all the cal- 
culations of to follow. Conversely, when individual vibration- 

rotation transition probabilities are studied, a duplicate and completely con- 
vergent basis set will always be used. 

6.6 Ambiguities of Unconserved Energy in a Three-Dimensional Semiclassical 
Model 

An example of the final rotational state occupations from a convergent 
basis set is shown in figure 6.14. The distributions in each vibrational 
manifold are plotted on vertically shifted scales to demonstrate their relative 
similarity. However, if we take the total energy of the system to be the 
Initial sum of internal energy and relative kinetic energy, then many of the 
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Figure 6.14.- Vibration-rotation transition probability distributions for 
C0(vil = 0,3)-He collisions from an effective Hamiltonian solution; 

E/ha)g = 1.08, b = 0, L = 0.02 nm, A = 1000 eV. 

rotational states in the v' - 1 manifold that remain occupied aftBr the col- 

lision are energetically inaccessible! (They are closed channels. Since a 
more exact collision model, where total energy is conserved, would leave the 
closed channels in all vibrational manifolds completely empty, the failure of 
our semiclassical model to comply with such considerations is, of course, a 
consequence of the lack of energy conservation when computing the classical 
trajectory, just as it was in the collinear model discussed in chapter 4. We 
showed in chapter 4 that, when the loss of energy conservation had no large 
effects on the molecular dynamics, it could be adequately compensated for by 

use of the average total energy given in equation (4.14) as 

= E + h(o)j^ + 03 j^)/2 (4.14) 
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However, the use of equation (4.14) to compute total energy makes it dependent 
on the transition being considered and, as a result, the concept of open and 
closed channels becomes ambiguous. While this ambiguity poses no practical 
difficulty, a careful interpretation of the semiclassical results must be made 
when comparing them to energy-conserving predictions namely, the probability 
distributions obtained for a single value of E do not correspond to a single 
total energy or single initial kinetic energy and, consequently, they cannot 
be compared directly with a distribution obtained from an energy-conserving 
collision model. But, after transforming the distribution of either collision 
model into the framework of the other, using equation (4.14), the results 
should be comparable just as they were in chapter 4. The accuracy of the 
three-dimensional semiclassical model in such comparisons remains to be proven 
and no comparisons are made in this study. However, from the results in 
chapter 4, there is good reason to expect that the three-dimensional semiclas- 
sical model, constrained to homonuclear molecules with light collision partners 
and properly interpreted, should give reasonably accurate predictions. 

Finall;/ , we must consider the consequences of the ambiguous closed- 
channel concept on our previous convergence arguments. While the initial dis- 
cussions of convergence in figures 6.11 to 6.13 were keyed to the concept of 
open and closed channels, they serve mainly as a guide to identify the number 
of channels for which convergence may be expected. In that regard, the concept 
is as useful as it would be in an energy- conserving model. Recall, however, 
that the primary aspect of our convergence criterion was to provide companion 
states in all vibrational manifolds to properly handle the transient d 3 mamics 
during the collision. That principle applies to all collision models, 
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regardless of their energy-conserving features or the final rotational-state , 
occupations obtained. 

6.7 Three-Dimensional Inelastic Collisions and Their Relation to Collinear 

Encounters 

The preceding three-dimensional collision model is applied here to evalu- 
ate the validity of collinear models for a variety of molecular types and 
initial conditions. However, examples are constrained to collision partners 
that can be treated accurately by the semiclassical approximation in the 
absence of dynamic coupling between the quantized molecule and the classical 
trajectory. Guided by the results in chapter 4 in that regard, we consider 
only helium atom collisions with H2, N2, and CO. The hydrogen molecule is an 
example in which the rotational eigenenergies are broadly spaced in comparison 
with the vibrational energies and the effects of vibration- rotation coupling 
are expected to be significant. H2 also requires the smallest convergent 
basis set since its homonuclear nature only couples rotational states of the 
same parity (i.e., AJt is always even). N2 also requires only states of the 
same parity, but its close rotational energy spacing places it in a different 
molecular class where numerous rotational states of high angular momentum must 
be included. Finally, CO is similar in rotational structure to N2, but its 
heteronuclear nature couples all rotational states of either parity. Most 
common inert-atom/diatomic-molecule collisions are represented by one of thfed^ ' 
three examples. Excluded cases are those in which reactive atom-exchange 
states participate in the energy- transfer process (as in H2-F collisions) and 
those in which electronic states participate (as in NO). 

The range of initial rotational states of interest in a practical appli- 
cation is first indicated by considering their equilibrium population 
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distributions. Figure 6.15 illustrates some typical distributions for several 
values of the parameter » T/6^, where 0^ is a characteristic rotational tem- 
perature with representative values of 2.8 K for CO and 44 K for H 2 
(0^ * hB^/sk, where s is a symmetry factor). The results in figure 6.15 
suggest that initial rotational states in the range = 0 to 20 are repre- 

sentative of most applications. We shall see that the collision dynamics for 
larger initial values of £ are easily inferred from the predictions for 
£ < 20 . 



Figure 6.15.- Maxwellian rotational state populations. Numbers at the distri- 
bution peaks indicate the corresponding angular momentum quantum numbers. 

Finally, we do not attempt to compare collinear and three dimensional 
predictions directly by obtaining rate coefficients from each model since to 
do so with the three-dimensional model would be too costly. Instead, we shall 
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compare the individual elements contained in the rate coefficient definition 
and their variations with the collision parameters, as discussed in sec- 
tion 6.2. 

6.7.1 H 2 *-He Collisions 

Although vibration-rotation energy transfer for H 2 -atom collisions has 
been studied extensively in the recent literature, ^ ^ ^ ^ we include 
it here as a contrasting example to the behavior of the heavier molecules that 
follow. The interaction potential constants, A and L, for all H 2 "He calcula- 
tions were chosen to resemble the calculated potential of Gordon and Secrest.®^ 

An indication of the coupling strength between vibration-rotation states 
in H 2 is given by the matrix elements shown in figure 6.3(a). As mentioned 
previously, the results show that transitions with A£ < A provide the primary 
path for vibrational energy transfer. The relative levels of rotational 
eigenenergies in each vibrational manifold shown in figure 6.4(a) then suggest 
that transitions from v = 0 to 1 (for example) will occur with increasing 
resonance enhancement from values of £ > 10 in v = 0, and that the state 
v£ = 0,16 is in near resonance with v'£' = 1,14 leaving very little excess 
energy to be traded with translation. Such transitions may be considered as 
the primary contributors to an apparent net V-R mechanism for energy transfer. 

The calculated transition probabilities shown in figures 6.16 to 6.18 
confirm all the foregoing expectations. Figure 6.16(a) shows that, when the 
initial state is v£ = 0,16, transitions occur to v’ = 1 predominantly for 
A£ = 2 downward with large probability and that the occupation of rotational 
states within each vibrational manifold is dispersed, but with less probabil- 
ity, by A£ = 2 transitions in either direction. A sequence of A£ = 2 down- 
ward transitions then populates v’ = 2, etc., at lower £. Conversely, when 
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the initial state has small angular momentum, as shown in figure 6.16(b), 
vibrational states are connected either by single large Ail transitions or by 
successive rotational transitions within v = 0 that precede any vibrational 
change. The vibrational transition probability is correspondingly lower. 

This later example then appears more as a V-T mechanism for energy transfer 
because of the relatively large amount of translational energy required to 
induce a vibrational transition. Intermediate results for other Initial rota- 
tional states are summarized in figure 6.17, where the net vibrational transi- 
tion probabilities are shown for all initial Z below the continuum. (Note 
that the molecular dynamics for initial states near the continuum are not 
accurately treated because of the neglected continuum interaction, but fig- 
ure 6.17 is believed to demonstrate a realistic qualitative behavior.) 

We can see from figure 6.17 that the vibrational transition probability 
is strongly dependent on the initial rotational state and hence a collinear 
model cannot be expected to realistically predict vibrational transition rates 
for The inapplicability of a collinear model to H 2 is further confirmed 

by figure 6.18. For example, curve (a) in figure 6.18 represents the usual 
collinear prediction in which the rotational contributions to eigenfrequencies 
are entirely excluded. The comparative three-dimensional predictions show 
that the probability and threshold energy depend strongly on the initial £ 
and both can be significantly different from the collinear results. Further- 
more, simple corrections to the collinear model, such as the use of vibrational 
frequencies shifted to match the predominant vibration-rotation states, are 
not satisfactory as curves (b) and (c) in figure 6.18 indicate. 

The preceding comparisons were all done for a zero impact parameter, but 
we have shown with equation (6.47) that, although collisions at b = 0 make 
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INITIAL STATE I 

Figure 6.17.- Effect of initial rotational state 
on the net vibrational transition probabili- 
ties for para-H 2 (v = 0)-He collisions, 
v' =1. Collision parameters are the same 
as for figure 6.16. 
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Figure 6.18.- A comparison of net vibrational 

transition probabilities for para-H 2 (v = 0)-He 
collisions from the three-dimensional and 
collinear models; v* = 1. Three-dimensional 
model Includes all states indicated for 
figure 6.16. 
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no contribution to the cross section or rate coefficient, the probabilities at 
b = 0 can be used as a normalizing factor. Cross sections are then propor- 
tional to P(E,0) and the integral of (b/L) P(E,b)/P(E,0) over all b/L. While 

v£->v* v£->v’ 

adequate discouragement for the use of a collinear theory to model H 2 colli- 
sions has already been presented, we display the cross-section integrands for 
later comparisons with heavier molecule results. Figure 6.19 shows the H 2 
cross-section integrand to depend strongly on initial Figure 6.20 shows 

it also to depend on collision energy but with significance only for initial 
rotational states where vibrational energy is transferred predominantly through 
V-R transitions (i.e., as in fig. 6.20(b) for v£ = 0,16). 



b/L 

Figure 6.19.- Effect of initial rotational state on the net vibrational cross- 
section integrand for para-H 2 (v = 0)-He collisions, v’ =1, at 
E/hojg = 1.5. Note that the integral of this parameter over all b/L is 
proportional to the inelastic cross section according to equations (6.46) 
and (6.47). 
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(a) H 2 (vJI = 0,2)-He, v' = 1 (b) H 2 (v£ = 0,16)-He, v’ = 1 

Figure 6.20.- Parameters affecting the net vibrational cross-section integrand 

for para-H 2 (v = 0)-He collisions, v’ = 1. 


From the preceding comparisons, there is no difficulty in concluding that 
the analysis of vibrational energy transfer in or any molecule where near- 
resonant vibration- rotation transitions with small A£ can occur, requires a 
complete three-dimensional treatment including both vibrational and rotational 
motions. Furthermore, even though some initial states with small angular 
momentum are treated by the collinear model with some resemblance to the three- 
dimensional results, extended use of the collinear model to deduce a thermal- 
rate coefficient enveloping a thermal distribution of initial rotation states 
is entirely inappropriate. 


6.7.2 N 2 “He Collisions 

The calculated collision dynamics of N 2 with He serve here as the primary 
example to illustrate the nature of vibration- rotation energy transfer in 
diatomic molecules with closely spaced rotational eigenenergies . The N 2 
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molecule provides a considerable computational convenience because, as with 
all homonuclear molecules, only alternate rotational states are coupled, 
thereby requiring half the basis set demanded by an otherwise similar hetero- 
nuclear molecule. Thus, in this section, we study N^-He collisions in detail 
and later compare the results with more limited calculations for the hetero- 
nuclear CO molecule. 

Unlike the H 2 -He collisions, the interaction potentials for atom colli- 
sions with many-electron molecules like N 2 are relatively unknown. The poten- 
tial parameters, A and L in equation (6.1), are therefore subject to large 
uncertainties and we must ensure that any conclusions made concerning the 
nature of energy- transfer processes are unaffected by the interaction uncer- 
tainties. An indication of the range of uncertainty is obtained by noting the 
range of interaction parameters implied in previous comparisons of experimen- 
tally determined vibrational rate coefficients with their related collinear 
theories.^^»^^»^^ For example, the interaction range, L, is typically found 
to be between 0.02 and 0.03 nm with L = 0.02 nm favored for molecules like 
N 2 * The analysis of molecular beam experiments,®^ which yields vibrationally 
inelastic cross sections directly, suggest similar values for L but also 
produces values of the interaction magnitude A. For molecules like N 2 , 
values are typically near A 1000 eV. However, an alternate means of 
obtaining A is to compare the exponentially repulsive potential model 
(eq. (6.1)) with the repulsive part of Lennard- Jones potentials (eq. (3.22)) 
implied by early viscosity measurements.®’®® Magnitudes as large as A=6 MeV 
are thus obtained. While this later value is not taken seriously nor is our 
potential model realistic enough to warrant much detailed interpretation, both 
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values are used here as limits to demonstrate the effects of A on the calcu- 
lated cross sections. 

The mechanism of vibrational energy transfer in N 2 is significantly dif- 
ferent from that previously shown for H 2 . The reasons may be generally under- 
stood from an examination of the eigenenergies and interaction matrix elements 
of N 2 (Illustrated in figs. 6.21 and 6.22, respectively). A well-known fea- 
ture of most induced transitions is that their probability is enhanced by 
their degree of resonance. For example, figure 6.17 illustrates the case for 
H 2 in which the most resonant transition from 5, = 16 is also the most prob- 
able, causing the principal path for energy transfer to be through V-R transi- 
tions. For N 2 , however, the vibration-rotation eigenenergies shown in 
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Figure 6.21.- Vibration-rotation eigenenergies for N 2 . 
(Note: Not all rotational states are shown.) 
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Figure 6. 22,- Time-independent matrix elements for para-N 2 ; v, v’ =0, 1, 

Open symbols are for L = 0.02 nm, solid symbols are for L = 0.03 nm. 

A similar plot for ortho-N 2 would have no distinguishable differences. 

figure 6.21 allow near-resonant transitions only for large A£, while the 
matrix elements for N 2 (fig. 6.22) suppress large transitions, just as 

they do in H 2 and CO (fig. 6.3). Thus, molecules with a vibration-rotation 
spectrum like N 2 , characterized by a large value for are always con- 

strained to nonresonant V-T transitions with a transition energy approximately 
equal to Consequently, the rotational energy transferred by small A£ 

transitions is always a small contribution to the total, thus rendering the net 
vibrational transition rates insensitive to the initial angular momentum. 

Figure 6.23 demonstrates that insensitivity for both para- and ortho-N 2 * 

The emphasis of nonresonant V-T transitions and the accompanying insen- 
sitivity of vibrational transition rates to the initial-state angular momentum 
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Figure 6.23.- Effect of Initial rotational state on the net vibrational trah- 
sltlon probabilities for para-N 2 (V = 0)-He collisions, v' = 1. Collision 
parameters are b = 0, L = 0.02 nm, E/htOg = 1.1. ' i 

I ' 

are the key factors contributing to a surprising accuracy of the collinear 
collision model for predicting transition probabilities in N^. Vibrational 
transition probabilities from both the three-dimensional and collinear model 
are compared in figure 6.24. Recall that the interaction potential is charac- 
terized only by the range L in the collinear model. Likewise, probabilities 
obtained from the three-dimensional model for zero impaot parameter depend 
only on L and hence are independent of A. The three-dimensional predictions 
in figure 6.24 therefore represent all values of A and the onty difference 
between the two models shown in figure 6.24 is the collision geometry. 

The role of A in the three-dimensional model is shown in figure 6.25(a), 
where the cross-section integrand is shown as a function of the impact param- 
eter. Clearly, A determines the variation of transition probability with 

impact parameter and thus determines the relation between P(b = 0) and the 

v£^v’ 

cross section. In another sense, A determines the equivalent elastic cross 
section (eq. (6.46)), undefined In the collinear model. 
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Figure 6.24.- A comparison of net vibrational transition probabilities for 
para-N 2 (v = 0)-He collisions from the three-dimensional and collinear 
models; v’=l, b-0, L=0.02 nm. Open symbols are three-dimensional 
model results. Shaded symbols are equivalent probabilities for a constant 
Op as defined by equation (6.46) but set equal to the cross section at 
E/hWe = 0.8, thus simulating the procedure typically applied to collinear 
models. 


When the interaction potential parameters are fixed, the equivalent 
elastic cross section (fig. 6.25(b)) increases slightly with collision energy. 
Any variation is contradictory to the assumption made when converting collinear 
probabilities to cross sections by use of a constant hard-sphere cross section, 
Oq (eq. (6.45)). However, as it turns out, the increase in approximately 
equals the increasing difference in P(b=0) between the two collision models 
so that their respective cross sections and rate coefficients are in closer 
agreement than figure 6.24 implies. The shaded symbols in figure 6.24 show 
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(a) Z = 20, E/hwg = 1.1 (b) A = 1000 eV, L = 0.02 nm 

Figure 6.25.- Parameters affecting the net vibrational cross-section integrand 

for para-N^Cv = 0)-He collisions, v* = 1. 


the equivalent position of the three-dimensional results corresponding to a 

e 

constant inelastic cross section obtained using the value of from 

E/ho)^ = 0-8 (i.e., cross sections and rate coefficients from the two models 
would compare graphically as indicated by the collinear model curve and the 
shaded symbols). 

A further point illustrated in figure 6.24 pertains to the required 
vibrational basis set for molecules like N 2 - The collinear model has been 
used as a guide to show that, although the predominant vibrational transitions 
are single-quantum transitions from the ground state to v' = 1, higher ^ ^ r , 
vibrational states participate even at collision energies near threshold. The 
three-dimensional model behaves similarly and it was therefore necessary to 
include four vibrational manifolds from v = 0 to 3 with £ = 0-60 in each 
before acceptable convergence was achieved- On the other hand, a collinear, 
first-order perturbation calculation of the type described in chapter 4 gives 
results nearly identical to the 11-state collinear model over the energy range 
included in figure 6.24. These comparisons suggest that a first-order 
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perturbation treatment of the vibrational motion, in conjunction with the 
"sudden approximation" describing the rotational motion, will be a very useful 
analytical method for dealing with three-dimensional inelastic collisions if a 
workable solution of the resulting integral equation can be found (e.g., recall 
sec. 6.4.1). 

Another interesting feature of the three-dimensional calculations is the 
variation of transition probabilities with impact parameter. As shown in 
figures 6.8 and 6.9, the net vibrational transition probabilities, 
simply decrease monotonically as the impact parameter becomes larger so that 
b = 0 impacts always produce the greatest probability of a vibrational tran- 
sition. However, the same is not true for all individual vibration- rotation 
transitions. The variation of with b is shown in figure 6.26, 

where some vibration- rotation transitions are seen to be more effectively 
induced by nonzero impacts, just as previous classical calculations have 
suggested. 



Figure 6.26.- Vibration- rotation transition probability variations with impact 
parameter for para-N 2 (v = 0)-He collisions at E/hG 0 ^=l.l, A =1000 eV, 

L = 0.02 nm. Collision radius is = 0.3 nm. 
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Finally, having established the applicability of a collinear model for 
predicting the rate of vibrational energy transfer from N 2 initially in the 
ground vibrational state, we can now investigate the accuracy of the collinear 
model in predicting the associated vibrational quantum number dependence of 
vibrational energy- transfer rates* The results validate the conclusions of 
chapter 5. Figures 6.27 and 6.28 compare both collision models for N 2 ini- 
tially in the v = 10 vibrational state and show them to be qualitatively 
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Figure 6.27.- Comparison of net vibrational transition probabilities for para- 

N 2 (v = 10)-He collisions from the three-dimensional and collinear models; 

V* * 11, b = 0, L = 0.02 nm. Shaded symbols are equivalent probabilities 

for a constant a equal to the cross section at E/hm = 0.5. 

o ^ 
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Figure 6.28.- Parameters affecting the net vibrational cross-section integrand 
for para-N 2 (v = 10)-He collisions, v* = 11. Collision parameters are 
A = 1000 eV, L = 0.02 nm, = 0.3 nm. 

similar as in the ground-state predictions. When the initial state is vibra- 
tionally excited, the differences in results are seen (fig. 6.27) to be of sec- 
ondary importance; that is, a larger vibrational basis set is required and the 
error of the collinear model is slightly greater. The three-dimensional model 
predictions must be compared with the five-state collinear model (dashed line) 
in figure 6.27 and, clearly, neither set of calculations contains a sufficient 
number of vibrational states. Thus, convergent three-dimensional calculations 
for excited vibrational states become increasingly impractical, even when the 

effective Hamiltonian approximation is incorporated. However, the variation 
e 

in with collision energy again compensates nicely for the inaccuracy of 

the collinear model as the solid symbols in figure 6.27 indicate. Thus, as 
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before, the colllnear model appears to produce an adequate description of the 
vibrational quantum-number dependence of cross sections and rate coefficients 
for molecules like N 2 * 

6.7.3 CO-He Collisions 

A comparison of figures 6.3(b) and 6.4(b) for CO with figures 6.21 
and 6.22 for N 2 suggests that the structural properties of CO are very similar 
to N 2 . The Intent of this section is therefore to look mainly for effects 
Introduced by the heteronuclear nature of CO associated with the additional 
coupling of even- and odd-parity rotational states. Generally, no significant 
effects were found and the conclusions reached for N 2 appear to apply equally 
well to CO. For example, figure 6.29 illustrates the net vibrational transi- 
tion probability dependence on il. For small the variations of 
are only slightly more pronounced in CO than for N 2 . Figure 6.30 shows that 
the colllnear model is just as applicable and the vibrational quantum-number 



Figure 6.29.- Effect of initial rotational state on the net vibrational tran- 
sition probabilities for C0(v = 0)-He collisions, v' =1, at E/hoig = 1.08, 
b * 0, L ■ 0.02 nm. Basis set includes I =» 0-60 for vibrational states 
V ■ 0-2. 
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Figure 6.30.- Comparisons of the net vibrational transition probabilities for 
CO-He collisions from the three-dimensional and collinear models. Collin- 
ear basis set includes states v = 0-17. Three-dimensional basis set is 
the same as shown in figure 6.28. Collision parameters are b = 0, 

L = 0.02 nm. 

dependence obtained with it is just as reliable. Finally, figure 6.31 demon- 

0 

strates that the equivalent elastic cross section, is not profoundly 

Influenced by the collision conditions, any more than it is for N 2 . Thus, we 
conclude that when the rotational eigenenergies are closely spaced in compari- 
son to the vibrational eigenenergies and, hence, no resonant transitions with 
small Mi are available, vibrational energy transfer will appear as a V-T 
process regardless of the other molecular properties. 
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C0(v=5)-He 


v’ = 4 





Figure 6.31.- Parameters affecting the net vibrational cross-section integrand 
for CO-He collisions; A = 1000 eV, L = 0-02 nm, = 0.3 nm, v’ - v = 1. 


6.8 Classification of Rotational Coupling Effects on Vibrational Energy 

Transfer 

The primary purpose of this chapter is to examine the role of coupled 
rotational motion in a diatomic molecule during the collisional exchange of 
vibrational and translational energy. A corollary to that purpose is the 
identification of conditions for which a collinear collision model will realis- 
tically predict the rate of vibrational energy transfer. We have found that, 
although the analysis of rotational coupling is complex, the nature of its 
influence on the energy- transfer process is conceptually simple. 

The mechanisms of vibrational energy transfer become readily apparent 
when one recognizes that the controlling features of vibration-rotation cou- 
pling are the predominance of coupling between states with small differences 
in angular momentum (i.e., small hZ) and its interplay with the resonance 
enhancement of transition probabilities. Thus, while we deal with collision 
energies that have many rotational states occupied in each vibrational mani- 
fold, their occupation occurs primarily through sequential rotational 
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transitions of small A£ during the encounter. Likewise, vibrational mani- 
folds are connected primarily by vibration-rotation transitions with small 
hZ* Hence, the predominance of small Ail transitions is a common feature of 
all inelastic mechanisms within the molecule. Furthermore, the range of Ail 
for which rotational states are closely coupled is independent of the degree 
of resonance associated with a transition and it is insensitive both to the 
interaction potential parameters and to the inertial properties of the mole- 
cule. Thus, the effective range of Ail (which we shall refer to as ^^max* 

where 0 < AJl < A£^^^) is similar for all the diatomic molecules and interac- 

i 

tidn potentials considered here and it is presumably similar for all others as 
weil. We have shown (fig. 6-3) that, typically, ^ 4 to 6. 

With small Ail transitions as a common characteristic, the different 
mechanisms of vibrational energy transfer separate into three natural classes. 
The molecules belonging to each class are identified by their inertial proper- 
ties, as specified first and foremost by the ratio of fundamental vibrational 
and rotational frequencies, and, second, by the proximity of the 

initial rotational state to a resonant companion. Given these two identifiers, 
we can then anticipate the qualitative nature of vibrational energy transfer 
for any diatomic molecule that does not involve electronic motion or reactive 
atom-exchange in the process. 

The first class pertains to all molecules in which ^^max> 

regardless of the initial rotational state. Our examples were CO and 
When the frequency ratio is very large, near-resonant vibration- rotation 
transitions of small AJl do not exist an 3 rwhere in the practical range of 
rotational states and any resonant enhancement of large AJl transitions is 
suppressed by the lack of coupling between such states. Consequently, the 
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energy- transfer process is dominated by nonresonant small M transitions in 
which the energy exchange is primarily between vibrational and translational 
modes. The initial-state angular momentum then has little influence on the 
rate of energy transfer and the process is described as a V-T mechanism. Mole- 
cules in this class are justifiably treated by a collision model in which the 
rotational contribution is either averaged or omitted. We found the collinear 
model to be surprisingly accurate for this class of molecules. 

The second and third classes of energy transfer pertain to molecules in 
which comparable with ^^max’ example was H2. Since 

is not large, near- resonant vibration-rotation transitions with small M are 
available. The proximity of the initial rotational state to rotational states 
capable of near-resonant vibration-rotation transitions then determines the 
class of energy transfer in which the molecule belongs. For example, the 
second class may be chosen as those molecules with initial rotational states 
remote enough from the near-resonant transitions for their angular momentum to 
be different by an amount greater than ^^xnax" Since the rate of energy 
transfer through the near-resonant transitions is rapid but proportional to 
the occupation of the resonant states, at least some of the resonant states 
must first become occupied before the near-resonant transitions can serve as 
an effective energy- transfer path. However, the restriction to small A£ 
transitions requires a sequence of induced rotational transitions to first 
take place within the vibrational manifold containing the initial state. Such 
a multistep process for successive small rotational energy changes during the 
collision is collectively as inefficient as a single-step nonresonant vibra- 
tional transition directly from the initial state. As a result, the energy- 
transfer process will not favor either path and the process must be labeled 
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as either a V-T mechanism in the extreme case or as a V-R-T mechanism in which 
near- resonant V-R and nonresonant V-T transitions compete. 

The final class also pertains to molecules with ^ ^^max> with 

an initial state within A^^iax ^ near-resonant vibration-rotation transi- 
tion with small A£. Near- resonant vibration- rotation transitions then 
immediately dominate the energy- transfer process. Pure V-R mechanisms of this 
type characteristically transfer vibrational energy at rates far exceeding the 
previous two classes. 

Clearly, these latter two classes involve the rotational motion of the 
molecule in a significant manner and a collision model omitting the rotational 
coupling would not distinguish their separate characteristics. However, we 
have shown that, while the collinear model is not applicable to such molecules, 
we can at least identify those molecules for which it may be applied by using 
simple identifiers. 
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CHAPTER 7 


A REVIEW AND SOME CONSIDERATIONS FOR FUTURE STUDY 

7.1 Review of the Newfound Aspects of Vibrational Energy Transfer 

The primary emphasis of this study has been directed toward the factors 
that influence the collisional exchange rates of vibrational and translational 
energy from excited vibrational states of diatomic molecules. In particular, 
emphasis has centered on the dependence of energy- transfer rates on initial- 
state vibrational quantum number and on the role of coupled rotational transi- 
tions in the energy-transfer process. As a consequence of the emphasis on 
excited vibrational states, two fundamental aspects were included in the col- 
lision model that are not often considered, namely, the anharmonicity of the 
molecular vibrations and the coupled interaction of multiple vibrational 
states. However, an overriding limitation to the realism of the collision 
model is the uncertainty of the interaction potential between collision 
partners. Thus, the conclusions of this study pertain mainly to the qualita- 
tive nature of vibrational energy transfer with no attempt made to predict 
absolute rates. 

A point made early in this study was that the analysis of a macroscopic 
nonequilibrium process is most conveniently carried out using simple analytic 
formulas to generate the necessary vibrational rate coefficients. One of the 
primary objectives of this study, therefore, was to evaluate the several ana- 
lytic approximations in popular use, as they apply to excited state transi- 
tions. To that end, a semiclassical description of the collision dynamics 
was adopted because, from previous comparisons, it showed the greatest promise 
as a theoretical framework leading to accurate analytical solutions. However, 
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a semiclassical formulation is itself an approximate description of the colli- 
sion dynamics. The first investigative step was therefore to determine the 
limitations of the semiclassical approximation when applied to a multistate 
anharmonic oscillator initially in an excited state. 

A comparative evaluation of the semiclassical approximation was conducted 
(ch. 4) based on a collinear collision model entirely equivalent to a fully 
quantum-mechanical formulation appearing in the literature. Transition 
probability predictions were compared for a wide range of anharmonic molecular 
types, initial states, and collision parameters. The comparisons also 
included heteronuclear molecules. Generally, they were unlike previous com- 
parisons in the literature which are typically confined to homonuclear 
harmonic oscillators in the ground vibrational state. The results illustrated 
some notable and previously unrealized effects of oscillator anharmonicity on 
the semiclassical approximation. For example, when computing the motion of 
the incident particle, the usual semiclassical procedure is to consider the 
oscillator nuclei as stationary relative to the molecular mass center. 

However, the effects of oscillator compression and recoil are amplified when 
the oscillator is anharmonic, and the agreement between semiclassical and 
quantal theories is significantly degraded unless the time-dependent average 
positions of the oscillator nuclei are introduced into the classical path 
determination. In the absence of such coupling, errors in the semiclassical 
approximation are largest when the incident particle mass is comparable to or 
larger than the mass of either molecular nucleus. When the molecule is 
heteronuclear, the semiclassical errors can become so large that anomalous 
resonances appear in the transition probability predictions. The anomalous 
resonances are strictly a consequence of the anharmonic coupling and do not 
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occur when the oscillator Is harmonic. Coupling of the enharmonic oscillator 
motion with the classical path removes the resonances and brings the semi- 
classical model into acceptable agreement with the exact quantal predictions. 
These results place new limitations on the use of analytic solutions based 
on a semiclassical approximation. For example, the analytic models do not 
include the effects of oscillator compression on the classical path and they 
should not be applied to heteronuclear anharmonic oscillators like the 
hydrogen-halides. For the same reason, the analytic models are also inaccu- 
rate when applied to the collision of light homonuclear oscillators with a 
heavy incident particle. This latter result is contrary to earlier arguments 
based on a concept that the semiclassical approximation should be most accu- 
rate for heavy incident particles because their wave packets are more local- 
ized and hence their motion corresponds more closely with a classical descrip- 
tion. Finally, the limitations of the semiclassical approximation observed 
for oscillators Initially in the ground vibrational state were found to be no 
more restrictive for oscillators in an excited state. Thus a semiclassical 
collision model should be adequate for studies of the dependence of fransition 
rates on the initial vibrational quantum number. 

With the semiclassical approximation validated and its limitations 
understood, the factors influencing the dependence of vibrational rate coeffi- 
cients on the Initial state quantum number were investigated next (ch. 5). 
Again, a collinear model was used, this time because it corresponds to the 
collision geometry adopted in all the analytic solutions to be tested. Com- 
parisons of the several analytic solutions available from the literature with 
multistate numerical solutions for anharmonic oscillators in excited states 
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showed that the most accurate analytical description in applications where 

the accuracy is important is a first-order perturbation treatment of anhar- 

monic oscillators. Conversely, the approximation in greatest popular use was 

found to give the poorest results for highly excited states. The influence 

of multiple-quantum transitions on the vibrational relaxation process from 

highly excited states was also examined. Two- and three-quantum transitions 

from highly excited states were found to be generally unimportant at kinetic 

temperatures less than the characteristic vibrational temperature of the 

oscillator (defined as hco /k) . Since vibrational relaxation is usually 

e 

superceded by other kinetic mechanisms at higher temperatures, the usual 
assumption that single-quantiom transitions prevail is adequate for the range 
of initial vibrational states typically considered. 

The remaining question pertains to the role of coupled rotational transi- 
tions in the transfer of vibrational energy and their impact on the previous 
conclusions obtained with a collinear collision model, A three-dimensional 
collision model was developed (ch. 6) that allows an arbitrary number of 
coupled rotational states to be included in the arbitrary set of vibrational 
manifolds. However, the collision calculations only confirm what is apparent 
(in retrospect) from the matrix elements associated with all diatomic 
vibration-rotation states: namely, that the two controlling factors of rota- 
tional coupling are (a) a restriction to vibration-rotation transitions with 
small changes in angular momentum and (b) the interaction of that restriction 
with the rate enhancement given to near-resonant transitions. Based on these 
general features, one can classify the mechanisms of vibrational energy 
transfer between a diatomic molecule and a structureless particle into three 
distinct types: vibration-translation (V-T) , vibration-rotation-translation 
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(V-R-T), and vibration-rotation (V-R) . The molecules belonging to each type 
are easily identified, first and foremost by their ratio of fundamental vibra- 
tional and rotational frequencies, and, second, by the proximity of 

their initial rotational state to a near-resonant transition invoking a small 
change in angular momentum. A result particularly important to the analyst of 
macroscopic kinetics is the finding that molecules belonging to the class 
dominated by V-T transitions (i.e., those where 1) accurately 

treated by a collinear collision model that may be reduced to yield analytic 
solutions. 

As a consequence of this study, we reach the broad conclusion that the 
collisional exchange rates of vibrational and translational energy can be 
accurately estimated for diatomic molecules in excited vibrational states 
using a simple analytic semiclassical model if the following conditions are 
met: 

(a) The fundamental vibrational frequency of the molecule is larger than 
its rotational frequency by several orders of magnitude. 

(b) The molecule is homonuclear or only slightly heteronuclear . 

(c) The incident particle mass is less than either nuclear mass of the 
molecule. 

(d) The Interaction potential is accurately modeled in the region of 
closest approach. 

(e) A theoretical model is used that includes oscillator anharmonicity 
in its primary formulation. The first-order perturbation treatment of Morse 
oscillators appears to be the most satisfactory choice. 

While the preceding conditions are numerous, they only exclude light or 
heteronuclear molecules, like H 2 and the hydrogen-halides, or collision 
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partners heavier than Argon. Hence the limitations of a well-chosen analytic 
model are not severe. 

7.2 Considerations for Future Study 

Long before this study was conceived, many Investigators were aware of 
the need for more detailed descriptions of the Interaction potentials between 
simple diatomic molecules and atoms or Ions. The computational and experi- 
mental determination of such potentials continues to be an activity of fore- 
most importance if quantitatively accurate predictions of vibrational and 
rotational energy-transfer rates are to be achieved. Clearly, many other 
related collision phenomena such atom exchange reaction rates. Ion and atom 
recombination rates, and colllslonal radiative line broadening also await the 
same potentials. However, the theoretical and experimental methods for deter- 
mining Interaction potentials are usually somewhat remote from the physics of 
Inelastic collisions discussed here. We only acknowledge their Importance 
to future studies of this type. In fact, with exact quantum-mechanical calcu- 
lations of Inelastic collisions now effectively a routine numerical exercise, 
much of the new work on vibrational and rotational energy transfer Is based 
on the availability of Improved potentials. 

Nevertheless, there are new practical applications, particularly those 
associated with lasers, that require further study Into several untouched 
aspects of vibrational energy transfer. For example. If we limit our Interest 
just to the exchange of vibrational and translational energy and exclude the 
multitude of other vibrational exchange mechanisms such as vibrational energy 
transfer between two oscillators or the interactions between vibrational and 
electronic states, we are still left with the following considerations: 
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(a) The absence of satisfactory experimental determinations of the depen- 
dence of vibrational rate coefficients on initial-state quantum number was 
indicated in chapter 5. Clearly, such measurements will be difficult, but 
recent improved techniques for selective excitation of upper states using 
txmable lasers and multiphoton absorption offer possibilities for new 
approaches . 

(b) The use of a semiclassical approximation and the deletion of long- 
range forces from the interaction have made the collision models of this 
study inappropriate for predicting low-temperature rate coefficients. Yet we 
have shown that the rate predictions for transitions from highly excited 
vibrational states are most sensitive to the collision parameters at low 
temperatures. Several infrared lasers of great practical importance operate 
at such conditions. Thus, rate predictions using a fully quantum-mechanical 
model and a more complete description of the long-range interactions would be 
extremely useful for the analysis of such lasers. The same collision models 
would also advance the study of heteronuclear molecules at all conditions 
since, for those molecules, the semiclassical approximation is generally 
inappropriate . 

(c) Molecules like H 2 are shown to transfer vibrational energy with high 
probability from rotational states near the continuum. Clearly, more realis- 
tic predictions of such energy-transfer rates should include interactions with 
the continuum. Furthermore, a collision model including continuum states 
would allow further study into the nature of vibration-dissociation coupling 
and also into the reverse process, namely, three-body recombination to excited 
vibrational states. 
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The foregoing considerations have come to mind during the course of this 
study as a result of the particular topics Investigated. However, recent 
Innovative techniques using the selective excitation of vibrational states 
In such applications as laser Isotope separation, photoenhancement of chem- 
ical reactions, and fluorescence enhancement have brought Importance to many 
other aspects of vibrational energy transfer not considered here. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 94035, March 4, 1976 
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APPENDIX A 


NOTATION 


The following catalog of symbols Includes only those used repeatedly. 
All symbols are defined locally In the text. Equations and figures cited 
locate explicit definitions and usage. 


a 



B 

e 

c 

n 

D 

D 

o 

E 

E 


E 


V 


E 


T 


h 

k 

k» 

Ir Ir 

m,n* nmi 


Morse intramolecular potential range, equation (3.3) 
trajectory coefficient for nonzero impacts, equation (6.36) 
interaction potential magnitude, equation (3.18) 
impact parameter, figure 6.1 

Morse oscillator wave*function constants, equation (3.10) 
molecular rotational frequency constant 
wave-function amplitude in basis state n 

Lennard-Jones interaction potential well depth, equation (3.22) 
Morse intramolecular potential well depth, equation (3.3) 
relative kinetic energy in a cent er-of -mass reference frame 
average relative kinetic energy 
oscillator energy in eigenstate v 

total energy of a colliding system in a center-of-mass reference 
frame, equation (4.15) 
h/2TT, where h is the Planck constant 
modified^ spherical Bessel function, equation (3.31) 

Boltzmann's constant 

oscillator anharmonicity constant, equation (3.9) 

rate coefficient for transitions from quantum state m to n 
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L 

i‘<y) 


angular-momentum orbital quantum numbers 
' ' Interaction potential range parameter, equation (3.18) 
Laquerre polynomial, equation (C.2) 
angular-momentum projection’ quantum numbers 


S 

”l 

j(J) 

•^I'm'lm 

N 

V 

Pj(y) 




r 

r 

e 



(r) 

V 


(J) 

w 


t 


m'-m 

mass of nucleous 1 

rotational matrix element coupling term, equation (6.22) 
oscillator radial wave-function normalization factor, 
equation (3.6) 

Legendre polynomial In the variable y 
transition probability between states m and n 
final transition probability from state v'l'm' to v)lm 
net vibration-rotation transition probability averaged over 
all Initial m states and summed over all m' stated, 
equation (6.42) 

net vibrational transition probability averaged over all 
Initial m states and summed over all V and m* states 
in manifold v' , equation (6.43) 
oscillator Internuclear separation distance 
oscillator equilibrium Internuclear separation distanee 

elementary radial overlap integral, equation (6.19) 
oscillator radial wave function for eigenstate v, 
eqtiation (3.6) 
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radial matrix element, equation (6.16) 

time measured from the instant of closest approach 



T 


u 


U 


v,v' 

V,V 

V 


V 

nm 




V® 

v’i'vil 

x,x^ 


X 

c 


X 



kinetic temperature 

average relative collision speed in a center-of-mass 
reference frame 

interaction-potential spherically symmetric term, 
equation (6.9) 
vibrational quantum number 

interaction potential between colliding nuclei 
spherically averaged interaction potential, equation (6.30) 
Intramolecular potential 

time-dependent Interation matrix element (overlap Integral) 
for states m and n 

time-independent matrix element for states v'i'm’' and 
vim, equation (6.26) 

''effectlve-Hamlltonlan'' matrix element, equation (6.55) 
separation distance between colliding nuclei, figure 3.4 
"hard sphere" collision radius 

mass-center separation distance between collision partners, 
figure 3.4 

mass-center separation distance at closest approach 
anharmonlc second-order frequency coefficient, equation (3.1) 
spherical harmonic function 

osclllator-lnternuclear separation parameter, equation (3.7) 


Y 

6 

LI 


molecular mass ratio, figure 3.4 
molecular orientation angle, figure 3.4 
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r(y) 

r(a,y) 

e,4> 


y 


y 


O 


P 


o 


a' 

o 


a 


e 


o 


a 

mm 


T 


T 

C 




0) 

m 


0) 

e 


0) 

inn 


Q 


gamma function of argument y 
incomplete gamma function 

spherical polar coordinate angles, figure 6.1 
reduced mass of the collision pair 
reduced mass of the molecule 

Lennard-Jones zero potential radius, figure 3.3 
”hard-sphere” constant cross sections, equation (5.4) 
equivalent hard-sphere cross section computed from the 
three-dimensional collision model, equation (6.47) 
total cross section for transitions from state m to n 
vibrational relaxation time constant, equation (1.1) 
mean collision-interaction time 
oscillator steady-state wave function 
perturbed oscillator time-dependent wave function 
circular frequency of eigenstate m 
fundamental oscillator frequency, equation (3.1) 

0 ) - 0 ) 

m n 

trajectory azimuthal angle, figure 6.1 



APPENDIX B 


GENERAL FORMULATION OF THE SEMICLASSICAL COLLISION THEORY 

The semiclassical procedure developed here is assembled from the contents 
of typical textbooks describing classical and quantum mechanics • For example, 
the classical equations of motion are derived from first principles in H, 
Goldstein, Classioal MedhanioSy Addison-Wesley (1950), chapter 3, and the 
quantum-mechanical methods are discussed in E. Merzbacher, Quantum MeahaniaSy 
John Wiley (1970), chapter 18. Both aspects are included here to provide a 
unified description of the complete theory and to identify explicitly the 
notation and assumptions associated with the collision model. 

The formulation to follow is based on a center-of-mass reference frame 
in which and y denote the reduced mass of the target molecule and the 

complete collision system, respectively. The incident particle is limited 
here to a structureless point mass whose motion is pictured classically. 

The target molecule is capable of intranuclear motion and its dynamic response 
to a time-dependent disturbance induced by the incident particle is described 
quantum mechanically. The motions of both collision partners are coupled 
through an interaction potential that depends on the relative separation of 
all nuclei in the system. For these purposes, the potential is represented 
here only by an arbitrary function V(q,R) when q specifies the molecular 
coordinates in configuration space and R locates the incident particle 
position relative to the molecular mass center. The remaining discussion may 
then be divided into a section describing the classical motion of the incident 
particle and a section detailing the quantum-mechanical formulation for the 
molecular dynamics. 
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B.l Classical Trajectory- 

In most semlclasslcal collision theories, the classical path Is deter- 
mined from just a spherically symmetric average of the Interaction potential 
centered on the target mass center. This approach reduces the encounter to a 
simple two-body central-force problem and, more Importantly, It avoids most 
of the difficulties arising otherwise from a need to define the molecular 
coordinate q, contained In the Interaction potential. In classical terms. 

The spherically symmetric average potential, V, may be obtained by quantum 
mechanically averaging the potential over all molecular coordinates In a 
manner suggested by 

V(R) = <j|v(q,R)|n> 


where molecular states |j> and |n> may be Initial states or some combina- 
tion of Initial and final states. In many cases, the method of averaging has 
little Influence on the final results. 

Given a central potential, the trajectory remains In a single plane 
described by two coordinates as shown In figure B.l. The potential may then 
be denoted as V(R) = V(x) and the trajectory Is conveniently described 
by a Lagranglan development of the equations of motion. As a result, the 
conservation of total energy, E, leads to 



+ V(x) 


(Bl) 


and the conservation of angular momentum, L, requires that 


L = 


-2 dfl 

dF 


(B2) 
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Figure B.I.- Classical path in a center-of-mass reference frame 
for a two-body, central force interaction. 

The initial conditions are defined by the initial speed, u^, and the impact 
parameter, b. At t = -<», we then have E = (l/2)iju? and L = ybu^. However, 
before introducing the initial conditions, we must recognize that no account 
has been taken of the energy or angular momentum traded inelastically with 
the target. This inconsistency is the origin of the lack of conservation in 
a semiclassical theory. While, in principle, a further approximation could 
be invented for keeping the system conservative, the usual method has been 
to adopt an equally approximate approach in which E and L are simply 
interpreted as "effective" constants of the motion averaged over the trajec- 
tory. In this formulation, we shall consider the impact parameters b and b*, 
shown in figure B.l to be identical but acknowledge that u^ ^ in a fully 
conservative system. We then define an average speed, u = u(u^,u£) and average 
energy E = 1/2 yu^, where the method of averaging is determined by that giving 
the best results. With these interpretations, the equations of motion that 
determine the time dependence of the coordinate R = R (x,fi) in the inter- 
action potential are obtained from equation (Bl) and (B2) as 
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(B3) 



(B4) 


B.2 Quantum-Mechanical Molecular Motion 

The molecular motion is driven by the complete interaction potential, 
V(q,R). To emphasize the time dependence, however, we use the equivalent 
notation V(q,R) = V(q,t). The Hamiltonian describing the molecular motion 
is then 

+ V(^,t) (B5) 

where 3f^(q) is the stationary-state Hamiltonian containing the intra- 
molecular potential according to 

3C(q) = - ^ + V (q) 

o 2vi^ q o 

This definition of 3f(q,t), in which the Incident particle motion is only 
implied by the time dependence of V(q,t), is the essence of the semlclasslcal 
approximation. Otherwise, ^f(q,t) would contain a momentum operator related 
to the incident particle. 

In the Schrodinger picture, the equation of motion is 

-ih = 3C(q,t)T(^,t) (B6) 

For t ±~, the molecule is undisturbed and giving a stationary- 

state solution of the form 

'i'o(q.b) c(0) e ” i|)^(q) + f dk (B7) 

n 
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where the probability amplitudes, and are constant in time. 

The first term in equation (B7) represents bound states of with quantum 

numbers n and satisfying 


JC t/; = ho) 

On n n 

The second term accounts for any continuum states allowed by with 

energy and momentum h^. While continuum states may exist, we justify 

their neglect by arguing that they will never participate in the dynamics 
of the molecule for the conditions of interest here. Consequently, we shall 
always choose total energies (internal plus kinetic) well below the level 
where any continuum states are energetically accessible. With that 
stipulation, the bound-state eigenfunctions, provide a natural and 

complete basis set in which to expand the solution to equation (B6) . Thus, 
we can write the time-dependent wave function as 

-id) t 


’I'Cq.t) = (t) e ^ (q) 

n n 


(B8) 


The probability amplitudes, c^(t) , are analogous to time 

dependent. The probability of occupation in state n at any time during 
the encounter is = |c^(t)|^. Since all that we desire are the 

occupation probabilities, a description of the texrms provides an 

adequate solution to the problem. Equation (B6) may be transformed into a 
set of equations — one for each ^ By the usual procedure of substituting 

equation (B8) and utilizing the orthogonality of the eigenfunctions, 

. The result is 

dc.(t) i(o) - 0 ) )t 

ih — ^ — = 2J ® ^ <j|V(q,t)|n> (B9) 

n 
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where the bracket notation refers to 

<j|v(q,t)ln> = J’^'j V(q,t) dq (BIO) 

integrated over all q space. 

To solve equation (B9 ) , the molecule Is considered to be Initially In a 
pure eigenstate |l), thus creating the Initial condition 

for all n. The Initial phase of c^(-“) is unimportant and is chosen 

arbitrarily since we are Interested only In Ic (t)l^. The final state of 

n 

the molecule at t -H» then determines the transition probabilities 
resulting from the collision. For transitions to state |j>, the probability 
is 

^i-^j “ 

At this point, note the role of various terms In equation (B9) and how 
they contribute to the transition probability. The matrix element, 

^3 I V(q,t) I n> , is a coupling factor that connects states |j> and |n>. It 
contains the primary quantal properties of the transition and it introduces 
the appropriate selection rules, if any exist. However, while selection rules 
will control individual matrix elements, the coupling of more than two states 
In the molecular dynamics can allow alternate routes for the molecule to 
reach a selected final state. Thus, energetic collisions Involving numerous 
Intermediate states will not always display the selection properties appear- 
ing in low-energy collisions where only two states participate in a transi- 
tion. 


194 



Resonance in a transition plays no direct role in determining the laatrix 
element properties. Thus, the matrix-element properties will prevail regard- 
less of the degree of resonance. However, resonance will have an additional 
effect on the transition probability by way of the phase term in equation (B9). 
Clearly, a resonant transition (in which Wj =* will not be degraded by 
phase interference during the collision and will achieve the maximum proba- 
bility determined by the matrix element. Conversely, as [wj - o)^l 
increases, phase interference can add an oscillatory structure to the final' 
transition probability that varies with the collisional parameters. 

The practical aspects of solving equation (B9) are made simpler if the 
time-dependent aspect of the matrix elements defined by equation (BIO) may be 
factored according to 


<j|V(q,t)|n> = U(t)<j|V(q)|n> 


(BID 


= U(t)V. 

jn 

Then the time- independent elements which contain all the quantum- 

mechanical selection properties and often require considerable numerical labor, 
can be computed in advance of the time-dependent solution. The function U(t) 
is obtained from the classical trajectory and applies to the entire set of 
equations (B9) for all quantum states. 

Another practical aspect in solving equation (B9) is to adopt the 
so-called "close-coupling” approximation in which not all eigenstates in the 
complete set are Included. Guided by trial solutions, only those states 
contributing to the dynamics of selected states of Interest are retained in 
the coupled set of equations (B9) . Usually, many states that are energetically 
inaccessible from a classical point of view may be neglected, although 
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experience has shown that inaccessible states with eigenenergies close to a 
state of interest contribute to the transient d 3 mamics during the collision 
even though they are unoccupied afterward. 

Finally, a necessary criterion that will be satisfied at any time, if 
the numerical solution of equations (B9) proceeds accurately, is the closure 
relation conserving probability, that is, 

T}|c(t)|2=l (B12) 

n ^ 

Hence, equation (B12) may be used as one test for regulating the numerical 
step size, although, in practice, it seldom becomes a limiting factor, 

B.3 First-Order Perturbation Theory 

When the kinetic energy is very low, the occupation of all states other 
than the initial state remain very small. Thus, if [i> denotes the initial 
state, we can assume |c^(t)l^ ~ 1 and |c^(t)|^ « 1 for all times and all 
n ^ i. To obtain the equivalent of a first-order perturbation theory in 
these circumstances, we must make the further stipulation that only two 
states interact, with one of them the initial state. This suggests that the 
eigenenergy of the second state is remote from all others and lies adjacent 
to the initial-state eigenenergy. If, for simplification of notation, we 
denote 

<n|v(^,t)|i> E V^^(t) (B13) 

equation (B9) can then be written: 



dc. (t) 



■Lh 

dt 

= V^^(t) c.(t) 

(B14) 


dc (t) 

•£(u -u.)t 



n'‘ ■' 
“(Tt 

= V'. . (t) c. (t) + V’ . (t) c (t) e " ^ 

11 1 ni n 

(B15) 
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The neglect of in equation (B14) constitutes the perturbation approximation 
and leads to a violation of closure as given by equation (B12) . The set of 
equations (B14) and (B15) are therefore not exactly equivalent to a two-state 
description using equations (B9). 

A solution of equation (B14) is 

c^(t) = exp 

thus fixing the occupation |c^(t)p = 1 
to suggest the form of by writing 


J V^^(T)dTj (B16) 

for all t- Equation (B16) is used 


where |cn(t)|^ = |bn(t)p. Thus, we may solve for b^Ct). Equation (B15) , 
in terms of bj^(t), takes on the convenient form 

. dbn(t) . r i 1 

^h = V* .(t) exp k(w - 0 ).)t+|- / (V - v:.)dx (B18) 

dt ^ ^ n ± h ^ nn ii I 

By defining a phase frequency F . = o) - cj. + (V’ - V*. ,)/h the transition 

probability is then obtained in simple form as 


i->n 


(B19) 


Note that since ^ from equation (B16), probability is not conserved 

and the solution of equation (B19) will be accurate only if 
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APPENDIX C 


RADIAL MATRIX ELEMENT INTEGRAL 

The radial matrix element expressed in equation (6.18) as a 

series expansion containing the Integrals where, according to equa- 

tion (6.19), 

OD 

The integral must be evaluated for a sufficient range of n to reach 

convergence in therefore seek an analytic solution to equa- 

tion (6.19) for arbitrary n that will allow rapid calculations of the 
numerous matrix elements required for a typical basis set of vibrational 
states. 

The first step in evaluating equation (6.19) is to express it in 
explicit algebraic terms. To that end, we recall that (r) ^^d ^^(r) 
are Morse oscillator radial wave functions described in chapter 3. Mbrse^^ 
shows that, in the absence of rotational centrifugal forces, the wave 
function may be written: 


where 


«^(r) 


N e 

V 


-Z/ 2 ^v/ ^ 


l;v(,) 


(Cl) 


a(r-r») 


b = k' - 2v - 1 
V 

1/2 

N^ = [ab^r(v + l)/r(k* - v)] 
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and V Is an Integer denoting the vibrational eigenstate quantum number, 
b 

The function L^^^(z) is a Laguerre polynomial defined as'^ 

in«0 

where 

(v) _ (-1)^ r(k* - v) 

\i m! (v - m) ! F(k' - 2v + m) 

and r(y) Is a gamma function. By comparing the elgenenergles associated 
with with the spectroscopic term expression: 


E /h = w (v + 1/2) - u X (v + 1/2)^ 

V e e e 

Morse shows that the parameter k' Is then a measure of the oscillator 
anharmonlclty and related to the spectroscopic parameters by 


k' ® 0) /o3 X (C3) 

e e e 

Substitution of the radial wave-function expression, equation (Cl), Into 
equation (6.19) leads to a transformed Integral over the variable z 
according to 


,(n) » ^ fe- 

w' \L/ w* (aL)^ 


z^ Lv^(z) L^y'(z) 


It' 


dz 


(C4) 


where 6 , is a Kronecker delta and the new terms are defined by 

w 

N , = N N , /a 

w V V 


^ V(k* - 2v - l)(k' - 2v' - pvlJT 
r(k’ - v) r(k' - v') 


(C5) 


and 


X = k’ - 2 - (v + v') 


(C6) 


200 



How invoking equation (C2) and noting that the product of two finite poly- 
nomials may be written as 


where 

equation (C4) becomes 


V V* « 

L (z)L , (z) = 2] V 
S=0 


(C7) 


j=0 


"j ^8-j 


W 


* 


-z X+s 
e 2 


Zn 




n 


dz (C8) 


Equation (C8) is an explicit algebraic equivalent of equation (6.19) for 
a Morse oscillator, and the reduced integral to be solved is readily apparent. 
Several exact analytic solutions of equation (C8) have been obtained in the 
past, but only for specific values of n. For example, Herman and Schuler^® 
found a solution for n * 1 that may be written as 

= -^ 6 - r(k* - v)r(k* - v) - v~ /pQN 

L w* aL v! (v* - v) (k* - v - v* - 1) f T(k* - v) ^ ' 

where v* > v. Heaps and Herzberg^^ extended the solution to n = 2 and 
indicated a procedure for obtaining solutions with larger n. However, the 
formulation for n « 2 is extensive and the implied formulation for larger 
n appears impractical for the required calculation of numerous n terms. 
Generalization of the solutions for arbitrary n also appears impractical. 
Hence, further considerations of exact analytic solutions to equation (C8) 
were abandoned. The iterative numerical procedure of Cashion^^^ was also 
rejected for similar reasons. Instead, we seek an approximate analytic 
solution based on the observation that the integrand in equation (C8) is a 
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I 


Ipcallzed function, uCpnflned. to a narrow, range of z when X Is large, . 
Since the value of., X.. Is. dominated by the .parameter k' and anharmpnlclby 
Is a second-prder. feature of all dlatpmlc molecules, k' Is alwa^^s large .. 
compared to unity. Hence an approximation, based on X, >> 1 will b@ 
applicable. 

The generality of the approximation to be made Is demonstrated In another 
sense, by, writing, the Integral In equation (C8) In the generalized form: 

(CIO) 


1(a) 




z a . . 
z u(z) 


dz 


where m(z) is an arbitrary function. The Integration., of equation (CIO), 
by parts proceeds according to 


1(a) 


00 

■I 


/ ^ do(z) , 
u(z) dz 


-p 


«(z)u(z)| u(z) dz 

uZ 


(Cll) 


'o ■'0 


where 


u(z) 


f -z a 
Je z 


dz 


(C12) 


and the Integral In equation (C12) Is Indefinite. To evaluate equation (C12) , 
repetitive Integration by parts leads to the series solution 

u(z) ■= -e"^z“jl + -| + + . . . j (C13) 

which may be recognized as the asymptotic expansion for the Incomplete 
gamma function^ 


u(z) = -r(o + l,z) 


(C14) 
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Equation (C14) is an exact relation, but It does not simplify the solution of 
equation (Cll). We now seek an approximation to u(z) that Is Integrable In 
equation (Cll). The nature of the approximation Is Indicated In figure C.l. 
Note In figure C.l(a) that the derivative 

. e"*z® 
dz 


Is very localized for large o. It has maximum at z <• a and a half-width 
at a half-height of dz/z ■ 'll Inlla.. Thus, the range of the Integrand In 
equation (Cll) becomes narrower as a Increases. Noting these features 
of du(z)/dz and the fact that u(») •• 0 from equation (C14) , the function 
u(z) approaches a step function as sketched in figure C.l(b). Thus, we can 
approximate u(z) by Introducing the Heaviside step function: 

f 0 f or y < 0 ^ 


H(y) 



1 for y > 0 


(C15) 



OL “Z 

Figure C.l.- Properties of u(z) and Its derivative du(z)/dz * z e 
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and writing 


u(z) « -o(0)H(a - z) (C16) 

where 

u(0) = -r(a + 1,0) 

= -r(a + 1) 


from the definition of the gamma function. Hence, we have, for a » 1, 

u(z) = -r(a + l)H(o - z) (C17) 

Substituting equation (C17) into the expression for 1(a) by equation (Cll) 
and invoking the step-function properties then leads to 

1(a) = -r(a + 1) j^-u(O) - /“^^dzl 


or, finally, 


1(a) = -r(a + l)u(a) 


(C18) 


Equation (C18) is a generalized approximate solution of equation (CIO) 
for arbitrary u(z) with only the stipulation that a >> 1, It may be 
applied to the solution of equation (C8) for A » 1 with the result: 

v+v' 

+ s + 1) 


^ y D r(A 

w \^/ vv (aL)^ s 

s— 0 




n 


(C19) 


From the definition of X given by equation (C6) , we see that the require- 
ment, X >> 1, is met when 

k' - (v + v') » 3 (C20) 

Since k' is always much larger than 3, the accuracy of equation (C19) will 
depend primarily on the sum v+v' and will decrease as the sum Increases. 
Correspondingly, the accuracy of equation (C19) is nearly Independent of n 
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since its derivation was done for an arbitrary function u(z), in this case, 
equated to [Zn(kVz)]’^, Thus, we can use the exact solution for n = 1 
given by equation (C9) to evaluate the accuracy of equation (C19) . 

An instructive first step in assessing the accuracy of equation (C19) 
is to examine its variation with k’ . We do so most easily by choosing the 
simplest and most frequently applied case of v' = 1 and v = 0. The exact 
solution for n = 1 is then 


[exact] 

while equation (C19) gives 

^01^ [approximate] 


-Vk* - 3 1 

aL k' - 2 



The relative error is simply 


Error = (k' - 2)Zn|-j-r 


k’ - 2] 


- 1 


(C21) 


Sample error values are tabulated below for some diatomic molecules covering 


a broad range of k' . 


Species 

k' 

aLrQp [exact] 

Error 

H 

37.25 

-0.1660 

0.0144 

2 




CO 

161.22 

-.0790 

.0032 


163.23 

-.0785 

.0031 

Br, 

302.05 

-.0576 

.0017 


Comparisons for H 2 (the worst case) at higher v and v’ (smaller X) are 
shown in figure C.2. As the figure shows, equation (C19) is acceptably 
accurate for single-quantum transitions (l.e. , v*-v = ±1) in for all 
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0 4 8 12 16 20 

V 

Figure C.2.- Radial integrals and their approximation for H„; 

k» = 37.25, aL = 0.509. 

initial eigenenergies at least up to half the dissociation energy of 
implied by the Morse potential. For other molecules with larger k', a much 
larger fraction of the total number of vibrational states will be treated 
accurately by equation (C19) . 

The final consequences of using equation (C19) to approximate 
are shown in figure C.3 where, again using H 2 as the worst example, we 
compare the radial matrix element » computed using equation (C19) , with 

numerically exact values obtained from a Gauss-Laguerre quadrature solution^ 
of equation (6.19). This comparison then encompasses values of ^ 

wide range of n. Recall that small J values correspond to small changes 
in angular momentum since |£*-£| < J < £ + ii; we showed in chapter 6 thac 
small contribute most to the energy-transfer process. As figure C.3 


206 



20 r 



Figure C.3.- Exact and approximate radial matrix elements for H,; 

k' - 37.25, aL = 0.509. 

Illustrates, equation (C19) Is exceptionally accurate for small J where 
Is the largest and most effective. Furthermore, the accuracy Is not 
strongly degraded even for large v where r^^ Is poorly approximated. 

Finally, In applying equation (C19) In a nvimerlcal calculation with 
large k' eind Increasing v+v', we find that the practical limit to its use 
Is not due to the error of the approximation but rather to the loss of 
numerical precision. More specifically, when k' Is large, as for CO or N^, 
the approximation Is basically very accurate to vfv* < 150, thus Including 
single-quantum transitions from Initial states to v •• 75. However, the terms 
In the summation of equation (C19) alternate In sign and the numerical range 
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between the largest term and the final value of the sum can exceed the largest 
number (with all digits significant) possible in most computers (l.e., as 
vW* Increases, we require decreasing differences between increasing numbers). 
For example, using a CDC-7600 computer with 28 digits in double precision, 
meaningful values of r^ from equation (C19) for CO (k' = 161.22) were 

obtained only up to v < 12 before all significant digits were lost. The 
numerical quadrature solutions were developed to obtain matrix elements for 
larger v. However, the orthogonality properties of r ^^ , are retained in 
the approximation and they may be used to at least monitor the numerical 
precision when using equation (C19) . To do so, one simply calculates 



or the equivalent surviving terms, 

v+v' 

Dr(X+s+l) =0 (C22) 

S^Q ® 

concurrently with r^ ^ ^ . The precision with which equation (C22) is 
satisfied is then a measure of the precision obtained with equation (C19) . 
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APPENDIX D 


SPHERICAL MATRIX ELEMENT INTEGRAL 

A general definition of the spherical matrix elements Is given by 
equation (6.15) as 

’^rm'Jlm“ <^'®’|Pj(cos 6)|£m> (6.15) 

where 6 is defined by figure 6.1 and Pj is a Legendre spherical poly- 
nomial of order J. In section 6. 1.2.2, the spherical harmonic wave func- 
tions of a rigid rotor are then introduced and the matrix elements take the 
specific integral form give by equation (6.21) as 


^ T ^ f TT ^ . 

-V2JTT J[ 

where angles 6 and <|) are polar angles (also defined in fig. 6.1). The 
purpose of this appendix is to derive an analytic solution of equation (6.21) 
using the properties of spherical harmonic functions commonly applied in 
angular momentum theories .104,105 

We begin by noting that equation (6.21) is similar to the integral of 
three spherical harmonics for which the solution is known to be (Edmonds, 
p. 63): 

-2TT -TT 

n, '^0 m '^0 n, (0.<(>)sin 6 d0 d<f. 

/z, Z z\/z Z„ ZA 

-yl(2Z,+l)(2Z^^l)(2Z, + l) Q oAmJm2in3) 

The bracket symbols are Wigner 3-j symbols. Equation (6.21) may be made to 

correspond to equation (Dl) by converting Y*, r(9><(>) to its complex conju- 

Xj m 

gate using the spherical harmonic property (Rose, 1*1^ p. 241): 
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'‘I’m- ‘ 

« nd by equating Y (6,0) to a new spherical harmonic In terms of 6 and 4>. 

o ^ 0 

To accomplish the latter, we recall the addition theorem (Rose,^^^ p. 60): 

ffl 

where the angles In equation (D3) are related by 

cos 6 - cos 6j^ cos 0^ + sin 0j sin 6^ cos(<J>j - <t>2) (D4) 

However, from figure 6.1, 6 Is related to the polar angles In equation (6.21) 
by 

cos 6 =* sin 0 cos(0 - (J)) (D5) 


80 that, by assigning 

0j »ir/2 , (J>j 
62“® » <t'2 * "f* 

the addition theorem may be written 

Vo'*’®* ■ <■>« 

m 

With the aid of equations (D2) and (D5), the integral In equation (6.21) may 
now be correlated with equation (Dl) , giving the solution, by Inspection, as 


,(J) 

^£*m*lm 


' ^2J + 1 


( 21 ' + 1)(21 + 1 ) 


0 0 0/ Z/ ^Jm(2*^)\-m' mm) 


m 

Equation (D7) Is further simplified by noting that the 3-j symbol coupling 
the projection states m* ,m and m Is nonzero only If 
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in “ in — m 


(D8) 


Hence the summation over m reduces to a single term and the matrix element 
becomes 


Jl’m'Am ^ 


-»"Vm ''3.„.-.(f.s)(o 0 oX-l- uLl) 


(B9) 


Equation (D9) is the desired analytic matrix element expression, but 

it is in a symbolic notation that requires further reduction to obtain an 

algebraic equation suitable for calculation. An algebraic form will also 

allow the classical time-dependent terms introduced by f^(t) to be isolated. 

To achieve an algebraic formula, we first evaluate the function 

Y_-(tt/ 2,^2) . A comparison of the Rodrigues formula for an associated Legendre 
Jm 

polynomial, P™(x) , with the definition of shows that the two are 

related by (Edmonds,^®** p. 24) 




+ 1 (A - x im3 

o.M - C-1) ^ ( - + „) ij * 


(DIO) 


Thus, we have 


r* o^ ( + 1 (J - m)!y/2 ffi 

[^7^ JT+ m) ! J ^ 


(t) 


(Dll) 


where (from Abramowitz and Stegun,^^ p. 334) 


- oin 

Pj(0) = — cos 


Jj + m + l^ 

ri» - 

L J * i) 


(D12) 


Equation (D12) may be reduced to simpler algebraic terms by noting that 
since J and Hi are Integers, the cosine term has the property: 
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cos 


-^(j - i) 


|(-1) ^ for J ±m even 
0 for J ±m odd 

Hence, nonzero matrix elements are obtained only for values of J ±m even. 
With that stipulation, the gamma functions in equation (D12) are also 
reducible as follows: Define a parameter z so that 


(D13) 


r ^J + m + ij ^ + 1/2) 


Then z = (J + m)/2 is always an integer because J + m is an even integer. 
With integer z, the duplication formula for gamma functions gives^^ 

/JT (2z) ! 


r(z + 1/2) = 




zl 


or 


/j + m + l\ /r ( J + m) ! 


(D14) 


Furthermore, for Jim even, the other gamma function in equation (D12) becomes 

. i) . (£^)! 

SO that 

+ in + 


/tT 


(J + 5) ! 


Hence, nonzero values of Y* (iT/2,n) are given for J ±m even by 

Jm 

J+m I ; :: — 

V* o\ = 2 j 2J + 1 /(J + m) 1 (J - m) ! -imQ(t) 

^JmU’ / ^ * 4 tt + mj i ^ - m^ | ® 


(D15) 


(D16) 


With equation (D16) , the time-dependent terms are easily isolated in the matrix 
element formula by defining a new matrix with constant elements given by 
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^(J) 


(-i>^ : Tn - w a2.' . 1H2, . i/ ’ " “Y‘' " 

2J(i^)i(i^)! Vo 0 oA-.’ s 


m; 


(D17) 


Again, the stipulations on that J ±m is even and m == - m. 

X, m 

Equation (D9) may then be written 


m(J) = -i(m'-m)jj(t) 

A'm'Am 


(D18) 


thus obtaining the results given by equations (6.22) and (6.23). 

Equation (D17) has been maintained in terms of the symbolic 3-j 

coefficients to simplify the notation, but they may be evaluated by the well- 

known formulas to follow^®**’ 

» J £ \ 

Q Q Q j is nonzero only for L 5 Z' + J + S, even. With 

that stipulation. 



(- 1 ) 


: J 


L/2 (L - 2i!,')l(L - 2J)!(L - 21) ! 


(L + 1) ! 




(L/2) ! 

(L/2 - i!.')!a/2 - J)!(L/2 - il) ! 

For even values of Jl' + J + £., we can also equate 



(2£' + 1)’^/^C(£J£' :m S m’ ) 


(D19) 


where C()lJJl' :m m m' ) is a Clebsch-Gordan coefficient defined by Rose,^°^ 
(p. 39), with a convenient algebraic expression. The result then leads to 
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/ ' ^ r (1* + A - J)ig' - t + J)ia + J- t»)lU* + m')Ki» - 

[,_j L (t + j + »' + l)!(t + m)I(l - m)l(J + 5>!(J - 1)1 J 

■>•»- v)ig - m -«- v)l ■ 

via* - t + J - v)!<r + m' - v)(v + £ - J - m')l 

V . : 

idiere v ranges over all Integer values giving nonnegative factorial argu- 
ments. Since the indiae constraints giving nonzero values of the vector- 
coupling coefficients (3-j symbols) are 


m' I < i , |m| < J , |m| < i 


i' - £| < J < i' + i 


(D21) 


m » m' - m 


the summation limits in equation (D20) are those that define the narrowest 
range of v within 


Mlnlinum viO,J-i + m’ 

Maximum viJ + ^f'+A,l'+m' 

Equations (D17) to (D20) are sufficient to calculate for all 

J,i*,m’,i,m combinations satisfying equation (D21) and 

i' + J + i even 
J ±m even 

All other matrix elements are zero. 
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APPENDIX E 


ANALYTIC TRAJECTORY EQUATIONS FOR NONZERO IMPACT PARAMETER 


We have defined a trajectory function In chapter 6 by equation (6.33) as 


U(b,t) = 


Vii 




-(x-5Eo)/L 


(6.33) 


where the constants,, are specified by the spherically averaged interaction 
potential ' 


V = A V 


11 


(El) 


Then, to solve the set of coupled dynamical equations describing the collision, 
U(b,t) must be determined explicitly in terms of time t. One approach would 
be to numerically integrate the classical trajectory equation given in appen- 
dix B as / _\2 


(f) 


(E2) 


but, as we point out in chapter 6, there are considerable advantages provided 
by an approximate analytic description of U(b,t) that benefit both the 
remaining numerical analysis and future analytic descriptions of the collision 
dynamics. In this appendix, we therefore develop an analytic form of U(b,t) 
by following the work of Hansen and Pearson. 

An indication of the functional form of U(b,t) is obtained by noting that 

.... 

equation (E2) , with the potential given by equation (El), may be solved 
exactly for b = 0. The result is 

[,-C*-Xo)/l] . sech2(|£) (E3) 


where x^ denotes the distance of closest approach. Furthermore, at closest 
approach for b = 0, all the Initial kinetic energy is converted to potential 
energy so that V = E, leading to 
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(E4) 


In view of equation (E3) , the primary nature of solutions to equation (E2) 
for nonzero impact parameters should be represented by 

exp|-[x(b,t) - XQ(b)]/L| = sech2|a^(b,t) (E5) 

where a^(b,t) is a slowly varying function of both b and t that can be 
approximated by a low-order expansion. 

By solving equation (E2) numerically, we can obtain the exact values of 
a|j(b,t) required to satisfy equation (E5). Some sample results are tabulated 
below for a representative collision energy and for small- and large-impact 
parameters. 


TABLE El.- EXACT VALUES OF a^j(b,t) FOR E/V^^ = 10"‘^ 


b/L 

ut/L 

V/E 

ajj(b,t) 

1.2 

0 

0.99 

0.994 


6 

.01 

.995 

8 

0 

.38 

.70 


6 

.015 

.74 

16 

0 

.0010 

.35 


6 

.0004 

.38 


Note that when the change in a|^(b,t) with t becomes noticeable, the inter- 
action potential is extremely small. We can therefore approximate a|^(b,t) 
without introducing significant error in the collision dynamics by assuming 
a|j(b, t)^ a^^(b,0) . Then, expanding both sides of equation (E5) about t = 0, 
we obtain, to first order, 

1 - [x(b,t) - x^Cb)]/! = 1 - (a|^ut/2L)^ (E6) 
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similarly 


S(b.t) - i„(b) + ^ (E7) 

From equation (E2), 



and 

-OWat)t,o ;|,2l 
* 

Since * -V/L from equation (El), equation (E6) finally leads to 

a,^(b,0) |l - [j^y[l - 2L/Xo(b)]|l/2 (E8) 

Note that x^Cb) must still be computed by iteratively solving equation (E2) 
for t ■ 0. However, as we show in figure 6.6, x^Cb) is closely approximated 
by equation (E4) when the Impact parameter is small. Similarly, when the 
impact parameter is large, the trajectory path is nearly a straight line and 
x^Cb) b. In that case, equation (E8) becomes 

a^(»,0) = /2L/b (E9) 

All three cases are compared in figure E.l, where we see that a completely 

1.0 
.8 
.6 

Ob (b.o) 

.4 

,2 

0 4 8 12 16 20 

b/L 

Figure E.l.- Trajectory coefficients for E/AV^^ = 10“**. 




217 



analytical approach would be to compute = 0) from equation (E4) for 

Increasing b until It gave values of a^^ equal to those from equation (E9). 
Beyond that point, a better approximation Is to assume x^ ^ b. In that case, 
U(b,t) ^ 0 according to equation (6.33) and the Interaction may be considered 
negligibly small. 
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APPENDIX F 


SYMMETRIES OF VIBRATION-ROTATION MATRIX ELEMENTS 
AND THEIR COMPACT COMPUTER STORAGE 

Experience with the collinear collision model made clear the facts that: 
the basis set required for a three-dimensional vibration-rotation model would 
be large, it would vary in size with the initial conditions, and it should be 
minimized for computing economy. A basic criterion of the computational 
scheme was therefore to permit an arbitrary basis set of vibrational and rota- 
tional eigenstates to be specified as part of the input information. Conse- 
quently, the calculation requires a large and variable number of time- 
independent matrix elements, V . . Since the matrix elements are con- 

stant in time, the obvious procedure was to compute them in advance of the 
dynamical solution and store them in the computer memory. 

When specifying the basis set in problems of this nature, the available 
size of accessible memory can be as severe a limitation as the computing time 
required. In large time-sharing systems, the operating cost is affected by 
both factors, while, in smaller systems, adequate memory volume is often not 
available. Hence the programmer's task becomes one of minimizing the memory 
volume that must be allocated to accommodate matrix elements of the largest 
basis set of interest. Since the allocation must usually be done in advance 
of any input Information, the storage scheme must also be optimized in advance 

Each matrix element is identified by six quantum numbers for which the 
simplest storage scheme would be a six-dimensional array. However, advance 
memory allocation for such an array would be extremely wasteful because each 
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dimension would have to be set to the largest value of Interest, 
suppose we choose the random basis set: 



For exanqple. 


(FI) 


where i is an index identifying the state, and v^ and are the vlbra*- 
tlonal and angular momentum quantum nijmbers specifying the state. If 
represents the largest values to be considered, then there are Vq + 1, + 1 

possible values of v, Z and 21^ + 1 values of m. To accommodate the 
example, let “5,3. Then the array would be dimensioned 

V(v' £'m' vim) “ V(6,4,7,6,4,7) , thus allocating 28,224 memory elements, while 
only 256 are filled by the sample basis set above. 

On second thought, a more efficient storage scheme is based on the index i 
Identifying each state and m^^, the projection quantum number. Two small 1-D 
(one-dimensional) arrays, Vj^(i) and £^^(1), may be established to give the 
v,£ quantum numbers when needed and the matrix elements are stored in a square 
2-D array, = V(j’,j), where j is conq>uted from 


j = 2 + 1) - £j^ + m^ (F2) 

k=i 

Such a matrix element array for the example basis set is Illustrated schemati- 
cally in figure F.l. The array dimensions, allocated in advance, are now 
required to be only as large as the total number of differential equations 
that can be solved in a reasonable computing time. Thus, the storage alloca- 
tion and the computing time limits are kept compatible. As .an example, suppose 
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Figure F,l.- A representation of the matrix elements correspond- 

ing to the sample basis set, equation (FI). Each group, i or j , repre- 
sents a pair of quantum numbers, and The matrix element is then 

identified by relating V^i£i^n*v£m array V(i,ra^, j ,mj) . The 

shaded squares are locations containing zeros. The filled circles • 
denote primary elements, unrelated by symmetry. The open squares and 
open circles c are elements related to the primary elements by the 
symmetry equations (F9) and Fll) . The open circles are additional ele- 
ments included in the index equation, equation (F17) , and stored with the 
primary elements in memory. 


that the maximum computing time limits the total number of states to 16. That 
limit would then encompass the example basis set and require a matrix element 
array with 256 elements. 

While the preceding storage scheme is a notable improvement, it is still 
extremely wasteful. Closer examination of the sample matrix element array in 
figure F.l reveals that, as a result of the constraints on rotational state 
coupling given by equations (6.24) and (6.25), almost half of the elements are 
zeros; In addition, approximately three quarters of the nonzero values are 
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numbers with opposing signs but duplicate magnitudes, as a consequence of the 
symmetries of there is an obvious opportunity to further 

reduce the storage requirements by storing only the elements of unique and 
nonzero magnitude. Such a scheme requires additional computing of indices to 
^relate an arbitrary set of quantum numbers, vH'm*vJ?,m, to the storage indice 
containing the appropriate matrix element. But, for large basis sets, that 
additional small effort is offset by a reduction in the number of elements 
that must be computed and by the savings in storage. We may use, as an 
illustrative example, the largest basis set represented in this study. It was 
for a heteronuclear molecule (CO) containing two vibrational manifolds with 
rotational states from iJ, = 0 to 10 in each one. Thus, 242 states were 
included and a array of 58,564 elements would have been required. 

However, by storing only the unique and nonzero matrix elements, the storage 
requirement was reduced to 14,883 elements. While the storage requirement was 
still large, the difference decided between possible and impossible storage 
allocation. The remaining paragraphs of this appendix are therefore devoted 
to a study of the symmetry properties of ^Y»£»^nfv£m select the 

unique elements and a derivation of the index equations for locating the 
matrix elements in a reduced storage scheme. 


F.l Symmetries of 

The symmetry properties of V , . , , . are revealed by the terms def in- 

V A/ m V A/m 

ing it. According to equation (6.26), 


^v’«,'m*vilm 


SL'+SL 

Z 


^ (2J + 1)eW> 


£m 


(F3) 
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where 


^(J) _ ^ - V ( J+m) / 2 


rAJ + m)!(J-m)!(2jt» + l)(2£ + l)~ [ A* ^ ^ 

I JVooo/UaJ 


(F4) 

m = m' - m (F5) 

R^! = <v’|f(r,J)|v> (F6) . 

and f(r,J) is an algebraic function of r and J given in equation (6.16). 

We first note that the vibrational states are freely interchangeable 
because the vibrational wave functions are purely real and not operated upon 
by the algebraic function f(r,J). Thus, equation (F6) may be rewritten to 
give 


= <v|f(r,J)|v»> = (F7) 

for all J. 

Next, the summation limits in equation (F3) are seen to be unaffected by 
an interchange of i and Hence, all the remaining symmetry properties of 

V I n determined entirely by the rotational coupling term, n • 

Furthermore, having generalized the vibrational symmetry, only three possible 
elementary symmetry operations remain: (1) an exchange symmetry between rota- 
tional states to relate „ and ^91, ,, (2) a sign reversal of m only, 

and (3) simultaneous sign reversals of m and m’ . All other operations would 
correspond to combined applications of the above. In the following, we deal 
with exchange and sign-reversal symmetries separately. 
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F.1.1 Exchange Symmetry 

The relation of \ « and i is easily shown in a general 

£’m’£m ^ ^ 

fashion by starting with its definition in symbolic notation. From chapter 6, 
Zm equivalent to 




-m)n 


(F8) 


where the combination of terms containing the time-dependent variables, S and 
Q, render the result constant in time. (Note that the bracket notation in 
eq. (F8) implies the integration over all configuration space.) Again, 
the operator Pj(cos 6) is algebraic and hence not operable on the wave func- 
tions. Thus, equation (F8) is unchanged when rewritten as 




But equation (F4) shows that the rotational coupling terms are always real so 




ZmZ m 


Z m £m 


(J) 

"£m£ 'm' 


Correspondingly , 


^v' £’m' v£m ^v£mv^ £'m' 


(F9) 


F.1.2 Projection-State Sign- Reversal Symmetry 

As indicated previously, sign reversal may be implemented in two ways. 

The first, a sign reversal of m alone may be immediately dismissed as an 
unsymmetric operation by noting that it would induce a change in the magnitude 
of m via equation (F5) and thereby lead to different magnitudes for 

n • The second case is a simultaneous sign reversal of 
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both m and m* . Only the sign of m Is then reversed and equation (F4) 
leads to 


(FIO) 


The corresponding 3-j symbol symmetry Is given by Edmonds (p.47) as 

/£l *2 * 3 \ _ ^2 ^3 \ 

\“l ®2 “3/ \~“l ““2 "“ 3 / 

but 1' + J + 1 (l.e. » Ij + I 2 ^3 Tn the above equation) must be even to 

( 1* J H\ 

Q Q 0/ equation (F4). Thus, the ratio of 
3-j S 3 nnbols In equation (FIO) Is always unity and we are led to the final 
result : 


V 


(-1)"^“V . 


(Fll) 


'v'l'm’vlm ^ V,£’,-m’,v,l,-m 

With the symmetries given by equations (F9) and (Fll) , one choice of pri- 
mary matrix elements Is Illustrated In figure F.l by the filled circles. In 
the notation of figure F.l, V(l,m^: j ,±m^) are Included for each 1 and all 
m^^ from -1^ to zero. All j and im^ are Included that fall to the right 
of the diagonal, with the exception of those related to preceding m^ In the 
same row by symmetry. All other matrix elements (Indicated by open spaces and 
open circles In fig. F.l) are then obtained by the symmetry relations, 
equations (F9) and (Fll). 


F.2 Primary Matrix Element Storage 

The remaining task Is now to devise a scheme of Indexing the primary 
matrix elements so that they may be stored and retrieved using the Identifier 
set (1, m^, j, method chosen here Is to Index them sequentially from 

left to right In figure F.l, starting with the top row and continuing, row by 
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row, toward the bottom. The primary elements may then be stored in a minimum 

memory volume by computing the index P(i,m, j>m.) and locating them in a 1-D 

J 

array, V(P). Similarly, the matrix elements are retrieved during the dynami- 
cal solution by again computing P(i,mj,j,m^) and applying the symmetry equa- 
tions (F9) and (Fll) when necessary. Hence, we need only define P(i,m^,j,mj) 
explicitly in terms of the identifiers to complete the storage scheme. 

Before developing the index equation, we first note that the formulation 
will be somewhat simplified if we slightly relax the requirement that aVL 
matrix elements related by. symmetry be excluded from the primary set. Very 
little redundancy is introduced by reinstating the few excluded matrix ele- 
ments to the right of the diagonal in figure F.l in rows where < 0. 

Such elements are indicated in figure F.l by open circles. With those ele- 
ments included, we compute P(i,m^,j,mj) by first defining the following 
component terms: 

An operator is required to be identified with each (i,m^) row and with 
the properties 

I I for il. + m. even 

0 for ii. + m. odd 
1 1 

Then the total number of primary and symmetric nonzero elements in row (i,m^) 
is 


3 ?. 


where I is the total number of states in the basis set. Similarly, 


the number of nonzero symmetric elements to the left of the diagonal in row 



- l)/2 


(F12) 


1-1 


3 i»®£ 1 1 i»®£ 


so that the total number of primary elements in row (l,m^) is 

N(i..,) . - D/2 

The number of rows preceding row (i,m^) is 

± 

nd.rnJ - 23 <21. + 1) + m - 1 - 1 

^ k-1 ^ 


(F13) 


(F14) 


so that the total number of primary elements in rows preceding row (l,m^) is 


then 


n(i,mj^) 

^ ] N(r,my) 
r=l 


(F15) 


Now choosing a specific element in row (i,m^), the number of nonzero primary 
and symmetric elements preceding V(l>m^,j,m^) is 

j-1 

Z (l. ^4 m ) <*'4 ^4 m " l)/2 

K.= l ^ ^’®i ^ J i 

while those symmetric elements to the left of the diagonal are again given by 
equation (F12) but rewritten as 

1-1 

2 

k= 

The difference in the two terms above is then the number of primary elements 


23 <^4 ^4 m ^ <^4 ®4 ^4 m “ l )/2 


preceding V(i,m^,j,mj) in row (l,m^), given by 

3-1 


<*k 


>/2 


(F16) 
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The Index of j »nij) Is now a combination of terms (F15) and (F16) with 

the result 

j-1 n(^l) 

» 1 + (Ij - + mj - m^)/2 + ^ NCr.m^) 

(F17) 

Equation (F17) requires that j > !• The identifiers may be exchanged to read 
P(j,mj,l»m^) if 1 > j. Note that the summations In equation (F17) may also 
be reduced to more efficient forms for computer calculation. 
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